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EXPERIMENTAL  INVESTIGATION  OF  THE  TURBULENT  LARGE  SCALE 
TEMPORAL  FLOW  IN  THE  WING-BODY  JUNCTION 
by 

Edwin  P.  Rood,  Jr. 

An  experimental  Investigation  of  the  fluid  dynamic  flow  In 
the  wing-body  junction  was  made  to  reveal  the  existence  of  large 
scale,  time-dependent  structures.  These  temporal  features  are 
discussed  relative  to  the  three  major  parts  of  the  mean  flow:  the 
onset  boundary  layer,  the  wing-body  junction  flow  characterized  by 
the  horseshoe  root  vortex,  and  the  wing-body  junction  wake  flow. 
Distinct  flow  structures,  distinguished  by  bandwidths,  were 
observed  in  all  parts  of  the  flow.  One  type  of  structure  was  due 
to  the  distortion  of  the  existing  structure  of  the  onset  boundary 
layer  by  the  wing.  A  second  and  a  third  structure  were  newly 
created  in  the  flow  by  the  presence  of  the  wing.  These  unique 
structures  along  the  body  were  boundary  layer  phenomena,  clearly 
distinguishable  from  ordinary  Strouhal-type  shedding  at  the  wing 
trailing  edge. 

The  measurements  distinguishing  the  structures  were  made 
using  applications  of  two-point  spectral  correlation  analyses  of 
the  flow  velocities.  The  application  of  this  established  procedure 
was  made  possible  by  modern  high  speed  data  techniques.  An  orderly 
framework  of  definitions  and  discussion  is  presented  around  which 
additional  research  can  be  performed  to  discover  the  details  of  the 
flow  structure. 
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CHAPTER  1 


INTRODUCTION 


1.1  Objective  of  the  Research 

The  objective  of  this  research  was  to  measure  the  large 
scale  temporal  structure  of  the  turbulent  boundary  layer  flow  about 
a  wing-body  junction  which  gives  rise  to  a  horseshoe  root  vortex. 
This  flow  is  a  complex  interaction  between  the  onset  boundary 
layer,  the  vortex  legs,  and  the  wing  wake.  Measurements  were  made 
to  confirm  the  existence  of  quasi-periodic  structure  and  to  deter¬ 
mine  some  of  the  governing  physics  of  the  temporal  features  of  the 
vortex.  These  temporal  features  are  not  evident  from  time-averaged 
measurements  of  the  flow  and  have  not  been  previously  studied. 
Consequently,  the  knowledge  gained  from  this  study  is  fundamental, 
and  should  prove  useful  for  guiding  further  research  applicable  to 
specific  problems  in  machinery  or  propulsor  design  such  as  unsteady 
(intermittent)  separation,  and  rotor  vibration  problems. 

In  the  context  of  this  investigation,  the  following  defini¬ 
tions  are  given  to  the  terms  large  scale,  low  frequency,  quasi- 
periodic  end  coherent .  The  meaning  of  large  scale  is  that  the 
organized  motion  of  the  structure  encompasses  a  spatial  extent 
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large  relative  to  the  wing  thickness  or  the  boundary  layer  thick¬ 
ness.  Low  frequency  Implies  a  temporal  periodicity  at  least  an 
order  of  magnitude  larger  than  the  ratio  of  the  wing  thickness  to 
the  flow  velocity.  The  temporal  structure  of  the  flow  Is  quasl- 
perlodic,  meaning  that  the  motion  Is  cyclic  for  short  time  dura¬ 
tions,  but  Is  uncorrelated  over  long  time  durations.  There  Is  an 
element  of  randomness  present  in  the  velocity  field.  The  flow 
is  coherent  within  the  meaning  that  organized  motion,  or  structure, 
superposed  on  the  mean  fow  is  spatially  and  temporally  confined. 
"Spatially  confined"  means  that  at  every  instant  in  time  the  defi¬ 
nite  phase  relationship  between  the  particle  velocities  in  the 
structure  is  over  a  subregion  of  the  entire  fluid  domain.  "Temp¬ 
orally  confined"  means  there  is  a  definite  time  scale  characterizing 
the  autocorrelation  of  the  structure. 

There  is  no  information  in  the  open  literature  directly 
addressing  the  large  scale,  low  frequency  flow  about  a  wing-body 
junction.  This  study  is  an  attempt  to  characterize  the  problem, 
which  is  best  approached  through  a  thorough  understanding  of  the 
complex  mean  flow  patterns  in  the  wing-body  junction.  The  mean 
flow  may  suggest  physical  mechanisms  (and  the  locations  where  those 
mechanisms  might  be  expected  to  be  observed)  that  govern  the  tem¬ 
poral  structure.  It  was  concluded  from  previous  studies  of  the 
mean  flow  that  unsteady  motion  could  be  expected  in  the  vicinity  of 
the  primary  three-dimensional  separation  at  the  wing  leading  edge, 
and  in  the  separated  flow  at  the  wing  trailing  edge. 

Preliminary  measurements  using  this  hypothesis  as  a  guide,  and 
a  subsequent  reevaluation,  have  led  to  the  premise  that  primary 
attachment,  and  possibly  secondary  separation,  at  the  leading  edge 
region  may  be  responsible  for  the  initiation  of  conditionally  stable 
temporal  structure.  Furthermore,  it  is  shown  as  a  result  of  this 
investigation  that  the  trailing  edge  pressure  condition  on  the  wing 
is  an  important  factor  in  determining  the  organization,  growth  and 
spatial  dominance  of  these  and  additional  instabilities  generated 
near  the  wing  trailing  edge.  Evidence  is  also  presented  that  the 
large  scale  time  dependent  structure  emitted  from  the  wing-body  junc- 
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tion  Is  in  part  related  to  the  large-scale  structure  in  the  onset 
boundary  layer. 

1.2  Experimental  Method 

The  wing-body  junction  flow  was  established  by  mounting  a 
streamlined  wing-like  protuberance  onto  a  flat  wind  tunnel  wall 
(the  body).  The  set-up  is  shown  in  Figure  1-1.  The  wing  was  a 
constant  two-dimensional  section  with  a  semi -elliptic  nose  profile, 
a  section  of  parallel  middle-body,  and  a  circular  arc  tail  profile. 
The  onset  wall  boundary  layer  was  generated  along  the  entire  inlet 
length  of  the  wind  tunnel,  was  thick  relative  to  the  wing,  and  was 
fully  turbulent  with  a  momentum  thickness  Reynolds  number  of  at 
least  5500.  The  detailed  description  of  the  experiment  can  be 
found  in  Chapter  3. 

A  pitot  tube  was  used  to  measure  the  mean  velocity  in  the 
onset  boundary  layer  and  an  x-array  hot  film  probe  was  used  to 
measure  the  velocity  components  at  selected  positions  around  the 
wing.  The  surface  static  pressure  was  measured  on  the  body  (the 
wall)  in  the  vicinity  of  the  nose  of  the  wing.  The  large  scale 
temporal  flow  was  measured  using  two-point  velocity  correlations 
obtained  with  single  element  hot  film  probes.  Specifically,  the 
auto-spectral  densities,  the  coherence,  and  the  relative  phase  were 
experimentally  obtained  and  used  to  determine  the  space-time  extent 
of  the  temporal  flow.  Oil  dot  and  oil  film  techniques  were 
employed  to  assess  mean  features  of  the  flow  to  provide  clues  to 
explain  the  physical  sources  for  the  temporal  flow. 


CHAPTER  2 


PHYSICAL  FEATURES  OF  THE  MEAN  FLOW 


2.1  Introduction  to  the  Global  Features 

The  flow  around  a  wing-body  junction  is  complex.  It  con¬ 
sists  of  a  turbulent  boundary  layer-pressure  gradient-viscous  wake 
interaction  that  Is  pressure  driven  in  some  restricted  regions. 
Reynolds  stress  driven  in  other  regions,  and  viscous  dominated  in 
most  other  regions.  The  essential  feature  of  the  interaction,  the 
horseshoe  root  vortex,  was  shown  in  the  schematic  of  Figure  1-1. 

The  figure  depicts  a  streamlined  protuberance,  the  wing, 
mounted  normal  to  a  surface,  the  body.  The  body  is  sufficiently 
large  that  a  turbulent  boundary  layer  exists  upstream  from  the 
wing.  The  wing  extends  out  from  the  body  a  distance  far  enough 
that  there  is  two-dimensional  potential  flow  around  a  considerable 
portion  of  the  wing  between  the  boundary  layer  and  the  wing  tip. 
This  specification  differentiates  the  wing  from  a  simple  roughness 
element  submerged  within  the  boundary  layer.  The  wing  is  stream¬ 
lined  so  that  there  are  no  flow  separations  attributed  to  break¬ 
away  at  a  body  corner.  More  precisely,  there  is  no  "two-dimension 


wj>  ^  .  ■  . 


■7.V  V  k.  ■ 


6 


al"  separation,  with  separation  line  parallel  to  the  wing  mid-chord 
line,  other  than  that  which  may  occur  very  near  the  trailing  edge 
as  the  wake  is  formed. 

In  the  simplest  conception  the  vorticity  present  in  the 
onset  boundary  layer  is  wrapped  around  the  leading  edge  of  the 
wing  to  form  the  horseshoe  root  vortex.  The  fluid  approaching  the 
wing  encounters  an  adverse  pressure  gradient  coincident  with  the 
presence  of  the  wing,  and  is  slowed.  Meanwhile,  fluid  off  to  the 
side  of  the  wing  passes  by  the  wing  without  experiencing  much 
effect.  Consequently,  fluid  elements  in  one  part  of  a  vortical 
line  element  are  decelerated  relative  to  fluid  in  other  parts.  The 
subsequent  wrapping  of  the  vorticity  line  element  around  the 
leading  edge  results  in  the  appearance  of  a  region  of  concentrated 
vorticity  known  as  the  root,  or  horseshoe,  vortex. 

The  effect  of  this  skewing  of  the  transverse  vorticity  is 
the  generation  of  streamwlse  vorticity  and  a  three-dimensional  flow 
velocity.  Whereas  the  flow  upstream  from  the  wing  has  one  com¬ 
ponent  of  velocity,  which  is  directed  in  the  streamwlse  direction, 
the  flow  in  the  root  vortex  has  components  in  a  plane  normal  to  the 
streamwlse  direction  due  to  the  existence  of  the  streamwlse  vor¬ 
ticity.  These  components  comprise  the  so-called  secondary  flow. 

The  flow  does  not  recover  to  its  original  state  downstream 
from  the  wing.  The  root  vortices  persist.  Furthermore,  there  are 
complex  interactions  between  the  vortices,  the  separated  flow  near 
the  wing  trailing  edge,  and  the  viscous  wake  of  the  wing. 

Throughout  all  of  this  disruption  of  the  flow,  the  root  vortices 
remain  inside  the  boundary  layer  from  which  they  were  generated; 
that  is,  they  do  not  lift  up  and  away  from  the  body.  Nevertheless, 
the  customary  boundary  layer  approximations  useful  in  flow  calcula¬ 
tions  are  invalid;  the  flow  is  described  by  elliptic  differential 
equations. 

The  complete  detailed  description  of  the  flow  is  quite 
complicated.  The  flow  will  now  be  explained  in  more  depth. 

2.2  Detailed  Description  of  the  Flow 

In  order  to  discuss  the  many  features  of  the  flow  about  the 


wing-body,  it  is  convenient  to  divide  the  flow  into  various  regions 
according  to  the  governing  dynamics  of  the  physical  features.  The 
various  regions  are  shown  in  Figure  2-1.  In  addition  to  the  single 
spiral  corresponding  to  the  primary  root  vortex,  there  are  addi¬ 
tional  spirals  corresponding  to  two  secondary  vortices.  These 
spiraling  flows  are  associated  with  three-dimensional  flow  separa¬ 
tion,  which  further  adds  to  the  complexity  of  the  flow.  The  high¬ 
lights  of  the  flow  regions  are  presented  below. 

To  a  considerable  degree,  the  description  is  conjecture, 
consistent  with  observations  and  measurements.  The  reason  for  this 
is  that  the  detailed  flow,  especially  in  the  separated  region,  has 
not  been  sufficiently  measured.  There  is  further  discussion  of 
this  observation  in  Section  2.3. 

The  discussion  below  is  strictly  applicable  only  to  tur¬ 
bulent  onset  boundary  layers.  The  flow  for  laminar  boundary  layers 
is  less  ambiguous  because  it  can  be  studied  with  low  speed  flow 
visualization  such  as  the  smoke  wire. 

2.2.1  Classification  of  the  Flow 

The  turbulent  boundary  layer  -  wake  interaction  flow  around 
a  wing  attached  to  a  body  consists  of  three  major  parts,  which  are 
further  divided  into  seven  regions .  Part  One  of  the  flow  is  con¬ 
cerned  with  the  onset  turbulent  boundary  layer  and  the  three- 
dimensional  separation  at  the  leading  edge.  It  is  in  this  part 
that  the  primary  root  vortex,  and  possibly  secondary  vortices,  are 
established.  The  flow  leaves  the  leading  edge  region,  and  the 
associated  strong  pressure  gradients ,  and  enters  Part  Two.  The 
exact  division  between  Part  One  and  P8rt  Two  is  arbitrary,  but  it 
is  clear  that  there  is  no  strong  pressure  gradient  in  Part  TVo, 
although  there  may  be  a  small  gradient  as  the  trailing  edge  is 
approached.  The  horseshoe  root  vortex  experiences  viscous  dif¬ 
fusion  and  the  effects  of  counterrotation  effected  by  anisotropy  of 
the  Reynolds  stresses.  The  flow  is  slowly  varying  in  the  stream- 
wise  direction;  the  flow  is  in  some  sense  in  equilibrium.  The  flow 
experiences  a  rising  pressure  gradient  as  it  nears  the  trailing 
edge  of  the  wing  and  enters  Part  Three.  The  exact  division  between 
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Part  Two  and  Part  Three  of  the  flow  is  also  arbitrary;  however,  the 
secondary  flow  in  Part  Three  is  rapidly  changing  in  the  streamwise 
direction.  Characteristics  of  the  flow  are  a  downward  movement 
along  the  wing  towards  the  body  and/or  a  movement  away  from  the 
wing  characteristic  of  imminent  flow  separation.  A  new  vortex 
springs  from  the  wing-body  junction  to  form  a  v-shaped  wedge  behind 
the  wing.  Inside  of  this  wedge  the  flow  is  viscous  wake. 

As  shown  in  Figure  2-1,  the  three  parts  of  the  flow  are 
denoted  the  Leading  Edge  Dominated  Flow,  the  Wing-Body  Junction 
Flow,  and  the  Wing-Body  Junction  Wake  Flow.  These  designations  are 
a  mix  of  terms  gleaned  from  the  literature;  there  is  no  set  ter¬ 
minology  that  has  evolved  yet.  In  fact,  the  problem  at  hand  is 
known  variously  as  the  wing-body  junction,  boundary  layer  -  wake 
interaction,  blade-end  wall,  and  numerous  other  names. 

2.2.2  Discussion  of  the  Flow  Regions 

REGION  I:  Zero  pressure  gradient  turbulent  boundary  layer  flow  far 
upstream  from  the  wing 

The  onset  boundary  layer  in  the  ideal  situation  is  two- 
dimensional  and  fully  developed  in  that  the  mean  velocity  profile 
follows  the  nominal  1/7  -  power  law,  and  the  momentum  thickness 
Reynolds  number  is  of  the  order  5000  or  larger.  A  spectral  analy¬ 
sis  of  the  flow  shows  no  peaks  or  "bumps"  on  the  otherwise  smooth 
distribution.  The  velocities  are  uncorrelated  in  the  horizontal 
plane  for  nominal  distances  beyond  one  boundary  layer  thickness  for 
the  frequencies  characterizing  the  temporal  flow.  There  may  be 
measurable  correlation  for  very  low  frequencies,  but  these  are 
clearly  separable  by  spectral  analysis  techniques  from  the  temporal 
features  of  the  flow  to  be  Investigated. 

REGION  II:  Pressure-driven  three-dimensional  boundary  layer  flow 
Immediately  upsteam  from  the  wing-generated  flow 
separation 

The  onset  turbulent  boundary  layer  encounters  the  adverse 
pressure  gradient  established  by  the  presence  of  the  wing.  The 
three-dimensional  pressure  profile  exhibits  a  ridge  extending 
upstream  with  diminishing  height,  and  sloping  off  on  either 
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side  of  the  ridge.  The  pressure  ridge  extends  upstream  a  distance 
of  order  10  wing  thicknesses  for  a  6:1  semi-elliptic  (major  axis) 
nose.  As  the  nose  is  made  more  slender  (the  ratio  n:l  made 
larger),  or  the  boundary  layer  made  thinner  (6/t  decreased  in 
value)  the  pressure  ridge  length  upstream  diminishes. 

The  pressure  gradient  has  two  effects  on  the  velocity 
field.  First,  the  pressure  gradient  produces  a  change  in  the  velo¬ 
city  distribution  that  decreases  the  streamwise  velocity  component 
and  increases  the  vertical  velocity  component.  The  effect  on  the 
vorticity  is  the  appearance  of  "true"  vortical  motion  as  6v/Sx 
obtains  non-zero  values  at  the  expense  of  Su/Sy  with  the  constraint 
that  the  vorticity  (Sv/6x  -  Su/6y)  remains  constant.  Second,  the  vor¬ 
ticity  vector  is  skewed  by  the  transverse  velocity  gradient  to 
generate  streamwise  vorticity  from  transverse  vorticity. 

REGION  III:  Pressure  generated  flow  separation  and  secondary  flow 
in  leading  edge  region 

This  region  includes  the  actual  three-dimensional  separa¬ 
tion  of  the  flow  and  the  concurrent  establishment  of  concentrated 
vorticity  in  the  so-called  horseshoe  root  vortex  which  wraps  around 
the  leading  edge  with  its  legs  trailing  downstream.  Several 
descriptions  of  the  separated  flow  appear  in  the  literature.  A 
reasonable  depiction  is  shown  in  Figure  2-2. 

» 

The  existence  of  the  secondary  attachment  A  ,  shown 
in  the  figure,  is  a  matter  for  debate.  It  may  be  Important  for 
explaining  the  physical  origin  of  temporal  flow  associated  with  the 
separated  region.  This  is  because  the  concentration  of  attachment 
fluid  in  a  small  area  is  one  clue  suggesting  the  existence  of  large 
scale  time  dependent  flow  structures.  In  this  small  area  there  are 
large  velocities  and  velocity  gradients.  Therefore  small  movement 
of  the  attachment  zone  would  produce  relatively  large  changes  in 
the  velocity  at  a  given  point  in  the  field. 

In  any  case,  there  is  no  doubt  that  a  primary  vortex  is 
established  in  the  leading  edge  region.  The  vortex  diameter  is  of 
the  order  of  the  boundary  layer  thickness  and  is  located  within  the 
region  on  the  order  of  the  boundary  layer  thickness  ahead  of  the  wing 
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leading  edge  for  wing  thicknesses  greater  than  the  boundary  layer 
thickness. 

REGION  IV:  The  developed  secondary  flow  downstream  from  the  leading 
edge  (Reynolds  stress  driven) 

In  this  region  the  pressure  gradients  in  the  streamwlse  direc¬ 
tion  are  relaxed.  The  vortex  diffuses,  and  at  the  same  time  rela¬ 
tively  small  secondary  flow  opposing  the  root  vortex  is  generated 
by  the  anistropy  of  the  Reynolds  stress  tensor.  (This  flow  is 
secondary  flow  of  the  second  kind,  or  stress-induced  secondary 
flow;  it  is  differentiated  by  Prandtl  (1926)  from  secondary  flow  of 
the  first  kind,  or  skew-induced  secondary  flow,  and  by  Young  (1977) 
from  secondary  flow  of  the  third  kind,  or  vortex-shedding  secondary 
flow. ) 

REGION  V:  The  separated  flow  at  the  trailing  edge 

Adverse  pressure  gradients  reappear  as  the  flow  encounters 
the  trailing  edge  region.  There  is  a  curving  flow  along  the  wing 
towards  the  body  and  the  secondary  vortex  in  the  junction  breaks 
away  from  the  wing,  although  it  remains  close  to  the  body.  A 
viscous  wake  is  generated  at  the  trailing  edge  of  the  wing. 

Immediately  outside,  or  at  the  edge  of  the  wake,  another  secondary 
vortex  is  established  as  a  result  of  the  downward  flow  along  the 
wing.  This  secondary  vortex  is  counterrotating  relative  to  the 
secondary  vortex  originated  at  the  wing  leading  edge. 

At  the  trailing  edge,  and  downstream,  the  primary  vortex 
legs  (one  on  each  side  of  the  wing)  are  free  to  mutually  interact. 
Because  they  are  large  and  transport  significant  amounts  of  stream- 
wise  momentum  across  the  flow  it  is  expected  that  they  contribute 
significantly  to  separation  features  in  the  wing  trailing  edge 
region. 

REGION  VI:  The  vortical  flow  in  the  wing  wake-boundary  layer 
interaction  flow 

The  four  counterrotating  secondary  vortices  (two  on  each  side 
of  the  wing  center  plane)  diverge  as  the  wake  in  the  boundary  layer 
expands  downstream.  The  primary  vortex  legs  continue  downstream, 
crossing  over  (on  top  of)  the  secondary  vortices.  The  mutually 


Interacting  vortices  and  the  separated  flow  from  along  the  wing 
span  provide  the  basic  features  of  the  flow.  The  temporal  flow 
Initiated  in  the  wing-body  junction  is  convected  into  this  compli¬ 
cated  flow. 

2.3  Review  of  the  Literature 

A  literature  review  revealed  that  there  is  a  considerable 
amount  of  experimental  data  applicable  only  to  the  mean  boundary 
layer.  Generally,  the  studies  can  be  classified  according  to  one 
of  the  three  flow  regimes  discussed  previously.  Studies  of  the 
flow  in  Part  One,  the  leading  edge  dominated  flow,  pertain  mostly 
to  flow  visualization  studies  using  oil  film  techniques  with  sparse 
supplementing  pressure  and  velocity  measurements.  Studies  of  the 
flow  in  Part  TWo  involve  extensive  velocity  measurements  to  examine 
the  diffusion  of  the  root  vortices  and  the  generation  of  secondary 
flow  of  the  second  kind  (Reynolds  stress  driven).  Studies  of  the 
flow  in  Part  Three  concentrate  on  oil  film  flow  visualizations  to 
assess  the  effects  of  the  vortices  on  trailing  edge  separation  and 
the  wake  growth,  and  on  velocity  measurements  to  evaluate  root 
vortex  growth. 

There  is  no  documentation  or  suggestion  of  any  studies  of 
the  large  scale  flow  structure  in  the  wing-body  interaction  flow 
for  a  turbulent  boundary  layer.  The  studies  of  turbulence  are  not¬ 
withstanding  because  they  are  restricted  to  measurements  with  zero 
time  and  space  separation  although  by  definition  "turbulence"  means 
time  dependent.  On  the  other  hand,  the  review  of  the  literature 
for  the  mean  flow  for  the  wing-body  junction  reveals  gross  features 
which  provide  clues  revealing  the  possibility  and  origin  of 
coherent  flow  structures.  It  is  expected  that  flow  regions  encom¬ 
passing  large  changes  in  the  mean  flow  quantltltes  would  be  asso¬ 
ciated  with  the  onset  of  large  scale  temporal  flow.  These  regions 
are  associated  with  flow  separations. 

The  applicable  literature  is  directed  toward  two  areas  of 
related  interest:  the  aerodynamic  wing-body  junction  and  the  tur¬ 
bomachinery  blade-end  wall  interaction.  The  aerodynamic  relevance 
is  for  the  local  drag  resulting  from  the  secondary  flow  and  for  the 
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downstream  influence  of  the  trailing  vortices  on  local  separation. 

The  turbomachinery  relevance  adds  the  effects  of  the  vortices  on 
heat  transfer  and  the  impact  of  the  vortices  on  the  performance  of 
downstream  blades.  There  is  no  single  reference  that  thoroughly 
addresses  all  three  Parts  of  the  flow.  However,  there  are  two 
references  that  stand  out  from  the  others  because  of  their  compre¬ 
hensiveness  although  they  are  lacking  detail  in  many  respects. 

These  primary  references  are  the  doctoral  dissertation  by  Shabaka 
(1979)  and  a  report  by  Sepri  (1973). 

Shabaka  investigated  the  flow  for  Parts  One  and  TWo,  with 
special  emphasis  on  the  relative  Importance  of  turbulence  terms. 

His  work  was  intended  for  input  to  the  development  of  computer 
prediction  methods,  although  he  provides  no  computations.  Consequently, 
he  measured  the  spatial  extent  of  the  Interaction  flow  and  ascer¬ 
tained  the  general  invalidity  of  law  of  the  wall  and  simple  tur¬ 
bulence  models  for  the  wing-body  junction.  His  measurements,  for  a 
6:1  semi-elliptic  nose  on  a  wing  of  very  large  chord  length.  Include 
wall  static  pressures,  mean  three-component  velocities,  and  single- 
point  double  and  triple  velocity  correlations  relevant  to  Reynolds 
stress  and  kinetic  energy  evaluations. 

The  significance  of  Shabaka* s  results  to  the  present  work 
is  the  information  about  the  physical  size  of  the  root  vortices. 

This  Information  was  required  to  choose  appropriate  instrumentation 
relative  to  the  physical  size  of  the  flow.  The  work  of  Shabaka  is 
noticeably  lacking  in  definition  of  the  precise  flow  in  the  vici¬ 
nity  of  the  leading  edge  separation.  However,  such  experimentally 
obtained  data  is  also  absent  from  all  other  references,  a  fact 
which  reflects  the  difficulty  Inherent  in  making  the  required 
measurements. 

Sepri  provides  the  most  complete  flow  visualization  analy¬ 
sis  published,  although  it  consists  solely  of  oil  film  studies. 

His  elaborate  investigation  of  an  NACA  0012  wing  mounted  on  a  flat 
plate  Includes  all  three  parts  of  the  flow,  and  In  particular 
demonstrates  the  effects  of  wing  Incidence.  The  notable  feature 
revealed  by  his  studies  is  the  flow  line  cf  attachment  In  the 
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leading  edge  Region  III  between  the  primary  separation  and  the 
wing.  The  distinguishing  feature  in  this  respect  for  the  wing  at 
incidence  Is  the  unmlstakeable  divergence  of  the  oil  film  streaks 
on  both  sides  of  that  attachment  line.  In  the  limit  as  the  inci¬ 
dence  angle  approaches  zero  degrees  the  attachment  line  moves 
toward,  and  appears  to  lie  very  close  to,  the  line  marking  the 
physical  intersection  of  the  wing  and  the  body. 

The  oil  film  studies  are  enhanced  by  three-component  mean 
velocity  measurements  in  the  trailing  edge  region  and  downstream 
of  the  wing.  Notable  again  is  the  effect  of  wing  incidence  on  the 
relative  sizes  of  the  root  vorticies;  it  is  observed  that  the  vortex 
on  the  high  pressure  side  of  the  wing  spatially  grows,  and  the  vor¬ 
tex  on  the  low  pressure  side  spatially  contracts,  with  the  result 
that  the  vortex  legs  remain  unequal  in  size  downstream  from  the 
wing. 

Seprl  made  no  flow  velocity  measurements  in  the  vicinity  of  the 
three  dimensional  separation  at  the  leading  edge.  Seprl  provides 
a  hypothetical  description  of  the  flow  in  the  separated  region, 
including  four  counterrotating  vortices,  two  separation  and  two 
attachment  lines,  and  an  interior  saddle  point  (where  flow 
streamlines  both  converge  and  diverge).  Although  the  description 
is  reasonable,  and  obeys  topological  rules,  it  is  not  substantiated 
because  flow  field  measurements  were  not  obtained. 

The  above  two  references  are  recommended  primary  sources  for 
information  about  the  mean  flow  for  the  wing-body  junction  flow 
because  they  are  the  most  comprehensive.  Among  the  important  refer¬ 
ences  supplementing  the  above  are  three  doctoral  disserations.  Oguz 
(1981)  examines  the  effect  of  boundary  layer  thickness  on  a  wing 
similar  to  that  of  Shahaka.  Baker  (1978)  examines  the  root  vortex 
about  a  circular  cylinder.  The  applicability  of  the  results  to  the 
present  study  is  restricted  by  the  bluff  trailing  edge.  Pal  (1981) 
examines  wake  development  downstream  from  the  wing-plate  junction. 
All  of  these  dissertations  are  deficient  in  detail  sufficient  to 
precisely  describe  the  three-dimensional  separations  both  upstream 
from  the  leading  edge  and  in  the  trailing  edge  region  of  the  wing. 


An  important  source  of  experimental  data  for  downstream 
velocity  profiles  is  a  report  by  Love  (1963).  It  is  the  earliest 
known  effort  to  systematically  examine  the  features  of  the  hor¬ 
seshoe  vortex  for  relatively  blunt  wing  profiles.  Love  measured 
transverse  velocity  profiles  across  the  flow  and  parallel  to  the 
wall  at  downstream  stations  behind  circular  cylinders,  circular 
arc  profiles,  and  elliptic  profiles.  Among  the  conclusions  are  two 
Import  ant  observations: 

1)  the  shape  of  the  wing  has  a  marked  effect  on  the  general 
nature  of  the  flow  downstream,  while  the  size  or  the  flow  velocity 
have  little  effect. 

2)  the  downstream  flow  is  strongly  dependent  on  the  point  of 
separation  of  the  boundary  layer  on  the  wing. 

These  conclusions  mean  that  vortex  development  at  the  leading 
edge  of  the  wing  is  conditioned  by  trailing  edge  factors. 

Love  is  also  to  be  credited  with  making  the  distinction  bet¬ 
ween  (a)  downward  flow  at  the  leading  edge  attributable  to  the 
total  pressure  gradient,  and  (b)  downward  flow  attributable  to  the 
horseshoe  vortex.  Although  Love's  explanation  is  somewhat  dif¬ 
ficult  to  follow,  the  distinction  between  the  two  is  apparently  the 
concentration  of  vorticity  in  the  upstream  region  near  the  body  in 
contrast  to  the  distribution  of  vorticity  far  from  the  body  but 
inside  the  boundary  layer  that  Impinges  on  the  leading  edge  region. 
Both  effects  produce  downward  flow  components. 

There  are  additional  references  that  provide  supplemental 
information.  Chu  and  Young  (1974)  continued  the  work  of  Sepri, 
especially  with  respect  to  the  effect  of  wing  incidence  on  the 
downstream  velocity  profiles.  Barber  (1978)  provides  insight  into 
the  individual  effects  of  wing  thickness  and  boundary  layer 
thickness.  He  specifically  discusses  the  possible  interaction  bet¬ 
ween  the  root  vortex,  trailing  edge  separation,  and  the  downstream 
velocity  profiles.  The  data  are  sparce,  and  serve  to  support  the 


conjecture  rather  than  to  experimentally  define  the  phenomena. 

Dechow  and  Felsch  (1977)  and  Pierce  and  McAllister  (1980)  have  made 
measurements  of  the  flow  upstream  from  the  leading  edge  in  Region 
II.  The  measurements  were  not  made  in  the  actual  vortical  flow 
inasmuch  as  the  Intent  of  the  measurements  was  to  verify  the  appli¬ 
cability  of  stress  models  for  predictive  techniques.  In  this  sense 
the  wing-body  configuration  was  a  device  to  produce  a  pressure- 
driven  three-dimensional  boundary  layer.  Another  reference  pro¬ 
viding  supplementary  information  with  respect  to  three  dimensional 
flow  separation  is  the  AGARDograph  by  Peake  and  Tobak  (1980).  The 
presumption  that  separated  flow  obeys  topological  rules  is  further 
investigated  in  papers  by  Tobak  and  Peake  (1979),  (1981). 

The  only  documented  attempt  to  directly  associate  flow 
field  measurements  with  oil  film  flow  visualizations  at  the  leading 
edge  Region  II  was  made  by  East  and  Hoxey  (1969).  The  data  are 
incomplete  although  an  apparent  relationship  is  revealed  between 
features  of  the  oil  film  and  the  location  of  the  minimum  static 
pressure  in  the  flow  field. 

The  flow  in  the  wing-body  junction  has  escaped  numerical 
study.  An  accurate  computation  would  involve  a  solution  of  the 
complete  Navler-Stokes  equation.  This  procedure  is  precluded  by 
computer  memory  and  time  limitations.  The  assessment  by  Sung 
(1983)  provides  a  critical  evaluation  of  the  problem. 

As  stated  earlier,  there  are  no  reported  measurements  in 
the  open  literature  of  time  dependent  flow  structure  in  the 
wing-body  junction  for  the  turbulent  boundary  layer.  On  the 
other  hand,  cyclic  flow  structure  has  been  observed  by  Peake,  et  al 
(1965)  and  by  Schwind  (1962)  for  the  laminar  boundary  layer.  Those 
authors  specifically  warn,  however,  that  the  observed  differences 
in  the  oil  film  patterns  and  the  known  differences  in  dynamics  bet¬ 
ween  laminar  and  turbulent  boundary  layers  might  lead  to  the  expec¬ 
tation  that  the  turbulent  flow  is  not  cyclic. 

One  physical  source  for  time  dependent  large  scale  flow  struc¬ 
ture  is  the  periodic  oscillation  of  the  three  dimensional  separa¬ 
tion.  A  series  of  journal  articles  written  by  Simpson,  et  al 


(1981),  and  (1983)  describes  comprehensive  experimental  results 
for  unsteady  two-dimensional  separation.  It  is  notable  that  ordi¬ 
nary  separation  produced  by  an  adverse  pressure  gradient  is  a  dyna¬ 
mic  time-dependent  phenomena  although  such  techniques  as  oil  film 
flow  visualization  imply  that  the  separation  is  fixed  in  location. 
These  references  are  listed  here  because  of  the  potential  applica¬ 
bility  to  the  three-dimensional  flow  about  the  wing-body  junction. 

Additional  references  not  discussed  above  are  listed  in  the 
bibliography  to  complete  the  list  of  references  on  the  wing-body 
turbulent  boundary  layer  Interaction. 

2.4  Rational  for  the  Temporal  Measurement 

Examination  of  the  mean  velocity  features  prompted  conjecture 
that  the  flow  in  the  wing-body  junction  contains  large  scale 
temporal  motion.  The  fact  that  the  flow  is  separated  was  suf¬ 
ficient  reason  to  perform  an  investigative  study.  In  the  absence 
of  geometry  fixing  separation,  such  as  a  sharp  edge,  the  flow  is 
generally  unsteady  for  several  simple  flows.  Consider  as  examples 
the  separating  flows  on  a  cylinder  and  on  a  flat  plate  subjected  to 
an  adverse  pressure  gradient. 

The  vortex  shedding  phenomenon  associated  with  a  cylinder  is 
well  known.  The  point  of  separation  for  the  flow  on  the  cylinder 
surface  is  a  periodic  function  of  time  as  the  wake  rolls  up  to  form 
and  shed  the  vortices  Schlicting  (1979,  p.28).  Peake  and  Tobak 
(1980)  summarize  data  showing  that  the  surface  location  for  separa¬ 
tion  for  inclined  axisymmetric  bodies  is  a  function  of  time.  It  is 
concluded  that  separation  is  associated  with  the  formation  of  tem¬ 
poral  flow  structure,  at  least  for  the  case  of  wake  formation. 

This  situation  is  certainly  applicable  to  the  trailing  edge  flow 
(Region  V).  However,  the  presence  of  body  curvature  in  the  above 
examples  is  not  necessary  to  produce  the  temporal  flow. 

Simpson,  et  al  (1981)  culminated  a  decade  of  research  with  the 
Important  observation  that  a  flat  plate  two-dimensional  separating 
flow  is  characterized  by  a  frequency  associated  with  the  presence  of 
coherent  flow  structures.  The  significant  difference  between  the 


flat  plate  temporal  flow  and  the  cylinder  turbulent  flow  is 
apparently  that  the  velocity  fluctuations  of  the  fluid  in  the  tem¬ 
poral  structure  are  small  for  the  flat  plate  and  large  for  the 
cylinder.  Although  the  separating  flow  for  the  wing-body  junction 
Is  three-dimensional,  the  flow  reversal  at  the  origin  of  the  root 
vortex  is  similar  to  the  two-dimensional  case  at  least  in  the  plane 
of  symmetry.  Furthermore,  it  would  be  expected  that  the  flow  in 
the  trailing  edge  region  of  the  wing,  and  inside  the  boundary 
layer,  would  experience  unsteady  separation  phenomena  that  would  be 
clearly  distinguishable  from  the  usual  trailing  edge  vortex 
shedding  that  occurs  outside  the  boundary  layer. 

The  conclusion  from  the  above  arguments  was  that  a  search  for 
turbulent  large  scale  temporal  flow  associated  with  the  wing-body 
junction  would  be  fruitful. 


CHAPTER  3 

DESCRIPTION  OF  THE  EXPERIMENT 


3.1  The  Low  Turbulence  Wind  Tunnel 

The  experiment  was  conducted  In  the  Low  Turbulence  Wind 
Tunnel  at  the  David  Taylor  Naval  Ship  R&D  Center.  The  general 
arrangement  of  the  tunnel  is  shown  in  Figure  3-1  from  Scottron  and 
Shaffer  (1965).  For  this  investigation  the  test  section  walls  were 
positioned  straight  and  in  parallel  although  they  are  deformable  as 
shown  in  the  figure.  The  center  of  the  test  section  is  charac¬ 
terized  by  a  turbulence  level  of  less  than  0.2  percent,  both  as 
reported  by  Scottron  and  Shaffer  and  as  verified  in  the  present 
study.  The  test  section  of  the  tunnel  is  4.4  m  long,  1.2  m  high, 
and  0.6  m  wide.  The  removable  test  section  top  contains 
transparent  panels  which  were  replaced  with  slotted  rectangular 
pieces  through  which  the  measurement  probe  traverse  assembly 
protruded;  see  the  photograph  in  Figure  3-2. 

The  maximum  tunnel  speed  is  about  40  m/s;  however,  the 
ambient  noise  level  outside  the  tunnel  at  that  speed  presents  a 
hardship  to  test  personnel.  Although  the  preliminary  exploratory 
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Figure  3-1:  General  Arrangement  of  Wind  Tunnel 
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invesdgatlon  was  for  a  speed  of  38.1  m/s,  the  speed  was  therefore 
reduced  to  30.5  m/s  for  the  detailed  measurements  without  an 
observed  change  in  the  flow  features  but  with  a  substantial  reduc¬ 
tion  in  noise  level  outside  the  tunnel. 

Preliminary  investigation  of  the  wind  tunnel  boundary  layer 
indicated  a  flow  structure  of  large  spatial  extent  with  a  6  Hz 
characteristic  frequency  at  a  tunnel  speed  of  38.1  m/s.  That  flow 
structure  was  removed  by  placing  a  6.4  mm  x  2. A  m  long  cylindrical 
rod  on  the  wall  of  the  large  entrance  section  on  the  downstream 
side  of  the  screens  as  shown  in  Figure  3-3.  The  rod  presumably  sta¬ 
bilized  otherwise  transient  conditions  produced  by  the  low  Reynolds 
number  inlet  flow. 

The  dynamic  pressure  determined  from  measurements  with  a 
pitot-static  tube  placed  in  the  tunnel  30.5  m/s  freestream  was  low 
frequency  spectrally  analyzed,  and  was  determined  to  have  a  maxi¬ 
mum  period  of  30  s.  The  spectral  output  is  in  Figure  3-4.  The 
integrated  pressure  spectrum  is  equivalent  to  a  velocity  amplitude 
fluctuation  of  2  percent  at  a  tunnel  free  stream  speed  of  38  m/s. 

The  relatively  long  tunnel  speed  fluctuation  period  established  the 
minimum  data  acquisition  time  for  mean  flow  measurements. 

3.2  The  Wing 

The  investigated  wing  profiles  had  semi-elliptic  noses, 
parallel-sided  mid-sections,  and  circular  arc  tails.  In  all  there 
were  two  nose  shapes,  three  lengths  of  parallel-sided  mid-section 
(including  none  at  all),  and  one  tall  shape.  Not  all  combinations 
of  geometry  were  investigated;  those  that  were  are  in  Figure  3-5. 

The  nose  shape  is  designated  by  the  ratio  a/b  of  the 
ellipse  semi-axes,  where  "a"  is  the  length  of  the  axis  in  the 
chordwise  direction  of  the  wing.  The  ratios  a/b  -  1.5  and  6.0  were 
examined  in  this  study,  although  wot  in  a  systematic  manner.  The 
a/b  -  6.0  nose  was  used  for  the  preliminary  studies  to  assess  the 
feasibility  of  studying  large  scale  time  dependent  flow  structure 
in  the  wing-body  junction  and  to  test  the  proposed  measurement 
techniques.  The  detailed  investigation  that  ensued  used  the 
a/b  -1.5  nose.  That  blunter  nose  produced  a  larger  root  vortex 
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3-5(a):  Profile  for  Prellalnsry  Investigation 


Figure  3-5:  Investigated  Wing  Profiles 
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which  lead  Co  spatially  more  extensive,  hence  more  easily  measured, 
flow  structure.  As  an  added  bonus,  it  turned  out  that  the  observed 
frequency  bands  associated  with  various  flow  structures  were  more 
easily  distinguished  for  the  a/b  -  1.5  nose  than  for  the  a/b  *  6.0 
nose. 

The  a/b  ■  6.0  semi -elliptic  nose  profile  for  the  preliminary 
investigation  was  directly  joined  to  the  28-degree  circular  arc 
tail;  see  Figure  3-5(a).  On  the  other  hand,  the  a/b  ■  1.5  semi- 
elliptic  nose  profile  had  a  length  of  parallel-sided  mid-section 
between  the  nose  and  the  28  degree  circular  arc  tall;  see  Figures 
3-5(b)  and  3-5(c).  The  total  length  of  parallel-sided  mid-section 
was  short  (2.25  wing  thicknesses)  for  the  finite  chord  length  wing 
and  long  (10.0  wing  thicknesses)  for  the  "infinite"  chord  length 
wing.  In  all  there  were  three  profiles  examined,  which  are 
designated  the  a/b  ■  6.0  finite  chord  length  wing,  the  a/b  ■  1.5 
finite  chord  length  wing,  and  the  a/b  -  1.5  Infinite  chord  length 
wing. 

Other  variations  in  the  wing  profile  design  geometries  were 
held  to  a  minimum  to  isolate  the  effects  of  nose  radius  and  wing 
chord  length.  -All  three  profiles  were  the  same  thickness  (50.8 
mm).  The  profiles  for  the  a/b  «  6.0  and  the  a/b  *  1.5  finite  chord 
length  wings  were  equal  in  chord  length  (254  nm),  with  a  chord 
length  to  wing  thickness  ratio  of  5.  The  profile  for  the  infinite 
chord  length  a/b  ■  1.5  wing  had  a  chord  length  to  wing  thickness 
ratio  of  17. 

The  wings  were  of  constant  profile  in  the  spanwise 
direction  and  extended  from  the  wall  a  distance  of  229  mm,  or  4.5 
times  the  nominal  boundary  layer  thickness.  The  wing  tips  were  cut 
square;  i.e.,  there  was  no  rounding  or  contouring  of  the  tip.  The 
loci  of  wing  leading  and  trailing  edge  points  were  along  straight 
lines  normal  to  the  wall.  That  is,  the  wing  was  a  right  cylinder 
normal  to  the  wall. 

The  wing  nose  and  tail  pieces  were  manufactured  from  wood 
and  were  finished  with  a  hard  sealant.  The  extended  parallel  mid¬ 
section  was  made  of  two  plastic  sheets  sandwiching  a  block  of  wood. 
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The  various  pieces  were  configured  to  dovetail  together  at  assembly 
in  the  wind  tunnel* 

The  wings  were  attached  to  the  wall  using  a  threaded  rod 
protruding  normal  to  the  wall  and  through  spanwise  holes  in  the 
wings.  A  nut  on  the  rod  was  used  to  tightly  press  the  wing  against 
the  wall.  The  spanwise  hole  was  located  through  the  nose  pieces 
for  the  two  finite  chord  length  wings.  A  clamp,  shown  in  Figure 
3-6,  was  in  turn  attached  to  the  threaded  rod  with  a  second  nut  to 
press  the  tall  piece  against  the  wall.  The  short  parallel  mid- 
section  for  the  a/b  -  1.5  profile  was  a  part  of  the  nose  piece  and 
therefore  did  not  require  clamping.  For  the  larger  infinite  chord 
length  wing  both  the  nose  and  the  tail  were  secured  by  threaded 
rod,  while  the  extended  parallel  mid-section  piece  was  clamped  at 
each  end  to  the  tips  of  the  nose  and  tail  pieces.  The  clamping 
arrangement  is  shown  in  Figure  3-7. 

Boundary  layer  trips  were  Installed  on  each  side  of  the 
noses  of  the  wings  a  short  distance  back  from  the  leading  edge. 

The  trips  were  0.5  ram  diameter  wire  taped  to  the  wing  at  spanwise 
intervals  to  hold  the  wires  down,  as  shown  in  Figure  3-8.  The  pur¬ 
pose  of  the  trips  was  two-fold:  first,  to  stimulate  turbulent  flow 
over  the  wing,  and,  second,  to  spatially  fix  transition  to  reduce 
or  eliminate  "contamination"  of  the  flow  by  unsteady  effects. 

Without  the  trips,  the  location  of  the  transition  would  have  been 
expected  to  fluctuate  and  produce  time  dependent  flow  structures 
that  perhaps  could  not  be  readily  distinguished  from  the  other  unsteady 
effects  of  Interest. 

3.3  The  Wall 

The  wall  on  which  the  wing  was  mounted  was  a  large  flat  plate 
of  clear  plexiglas.  This  plate  was  installed  parallel  to  the  wind 
tunnel  side  walls,  rather  than  to  either  the  top  or  the  bottom 
wall,  to  minimize  "end  wall  effects"  on  the  flow  around  the  wing.  As 
described  above,  the  tunnel  side  walls  were  122  cm  in  height 
whereas  the  tunnel  top  and  bottom  walls  were  only  61  cm  in  width. 

The  plexiglas  plate,  which  effectively  replaced  the  wind 
tunnel  side  wall,  extended  from  the  tunnel  top  to  the  tunnel  bot- 


Figure  3-6:  Photograph  of  Compression  Clamp  Holding  Tall  Piece 
for  Finite  Chord  Length  Wing 
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Figure  3-7:  Photograph  of  Compression  Clamps  Holding  Extended 

Parallel  Mid-Section  Piece  for  Infinite  Chord  Length 


ton.  The  wing  was  mounted  normal  to  the  plexiglas  wall  and 
protruded  horizontally  Into  the  tunnel  flow,  with  equal  distances 
between  the  wing  and  the  tunnel  top,  and  the  wing  and  the  tunnel 
bottom.  The  wing  and  wall  configuration  Is  shown  In  Figure  3-9. 

The  plexiglas  wall  was  Installed  with  a  clearance  between 
It  and  the  tunnel  wall  to  permit  Installation  of  pressure  tap 
tubing.  As  shown  In  Figure  3-9,  an  upstream  ramp  was  used  to 
literally  lift  the  tunnel  boundary  layer  from  the  tunnel  wall  and 
onto  the  plexiglas  wall.  This  procedure  established  a  relatively 
thick  onset  boundary  layer  for  the  wing,  and  eliminated  problems  of 
flow  alignment  at  the  plate  leading  edge.  A  turning  block  of 
radius  61  cm  was  used  to  turn  the  flow  from  the  ramp  and  onto  the 
plexiglas  wall.  Oil  dot  and  tuft  flow  visualization  studies  showed 
no  evidence  of  flow  separation  on  the  wall  downstream  from  the 
block.  Furthermore,  laterally  displaced  two-point  velocity  corre¬ 
lations  did  not  show  non-zero  values  indicative  of  flow  separation. 

The  plexiglas  wall  terminated  downstream  at  a  support  block 
which  presented  a  backward  facing  step  transition  for  the  flow. 

The  step  transition  was  Immediately  upstream  from  the  divergent 
flow  outlet  for  the  test  section.  Certainly  there  was  Unsteady 
flow  downstream  from  the  block  and  into  the  divergent  outlet,  but 
Its  presence  was  not  evident  In  the  upstream  flow  from  measure¬ 
ments  In  this  Investigation. 

Plexiglass  was  chosen  for  the  wall  material  to  allow  obser¬ 
vation  of  the  numerous  pressure  tap  tubes  to  verify  their  proper 
connections  when  the  tunnel  wall  was  installed  in  place.  The 
plexiglas  material  also  made  the  wall  easier  to  handle  than  a 
heavier  metal  wall. 

The  wing  selected  for  study  was  initially  aligned  geometri¬ 
cally  to  the  assumed  wind  tunnel  flow  direction  by  locating  It 
along  a  scribe  line  on  the  plexiglas  wall.  The  scribe  line  was 
located  from  the  locus  of  points  equidistant  from  the  tunnel  floor. 
Changes  in  the  wing  incidence  of  12  deg.  from  the  assumed  zero 
degree  alignment  produced  relatively  small  measurable  effects  for 
the  time  dependent  flow  near  the  wing.  It  was  concluded  that  pre- 


else  alignment  of  the  wing  was  not  necessary  to  quantify  the  flow 
near  the  wing  for  the  purpose  of  this  investigation.  However,  the 
preliminary  time  dependent  flow  measurements  several  wing  chord 
lengths  downstream  from  the  wing  were  observed  to  be  sensitive  to 
small  changes  in  incidence  angle.  For  those  measurements  the  zero 
incidence  alignment  of  the  wing  was  verified  by  obtaining  symmetric 
measured  mean  wake  profiles.  These  observations  of  incidence  sen-* 
sltivlty  mean  that  a  more  precise  method  of  determining  wing  angle 
of  incidence  may  have  to  be  developed  for  precision  wake  measure¬ 
ments. 

3.4  The  Coordinate  System 

The  coordinate  system  for  the  experiment  is  shown  in  Figure 
3-10.  The  coordinate  directions  are  defined  similarly  to  the  con¬ 
ventions  used  by  Shabaka  (1979)  and  by  Oguz  (1981)  in  previous  stu¬ 
dies  of  the  mean  wing-body  junction  flow.  The  coordinate  system 
origin  is  located  on  the  wall  at  the  intersection  with  the  wing 
leading  edge.  The  coordinate  system  is  identified  with  customary 
boundary  layer  directions  in  Figure  10(a).  The  positive  direction 
of  x  is  in  the  longitudinal,  or  streamwlse,  direction,  and  the 
positive  direction  of  the  coordinate  y  is  across  the  boundary  layer 
and  normal  to  the  wall;  l.e.,  in  the  transverse  direction.  The 
coordinate  z  is  defined  in  the  direction  across  the  boundary  layer, 
perpendicular  to  the  mean  flow  direction  and  parallel  to  the  wall; 
i.e.,  in  the  lateral  direction. 

As  in  the  previous  studies,  the  coordinate  z  in  the  nota¬ 
tion  P^  (x,y,z)  specifying  the  location  of  the  point  PA  is  defined 
relative  to  the  wing  surface.  Hence  otherwise  possible  positions 
P^  inside  the  wing  are  excluded  from  the  coordinate  range  of  P^ 
(x,y,z)  where  x,y,z  can  assume  all  values  from  -»  to  +». 

3.5  Probe  Positioning  Devices 

The  temporal  flow  was  measured  by  correlating  flow  velocity 
fluctuations  at  two  probes,  one  of  which  was  fixed  in  position 
while  the  second  was  traveled  to  sequential  positions. 

The  traverse  for  the  traveling  probe  consisted  of  two  devi¬ 
ces  to  provide  independent  probe  positioning  in  the  y  and  z  direc- 


3-10(a):  Boundary  Layer  Dlrectlona 
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Figure  3-10:  The  Coordinate  System 
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tions  from  outside  the  tunnel.  The  first  device  was  a  massive  mill 
table  mounted  outside  the  tunnel  and  supported  by  the  twin  I-beams 
forming  the  tunnel  top  framework;  see  Figure  3-11.  The  mill  table 
provided  continuous  positioning  in  the  y-direction  with  a  scale 
resolution  of  0.0001  in  (0.00254  mm).  The  second  device,  also 
shown  in  Figure  3-11,  was  a  reciprocating  ball  guided  plate  with 
scale  resolution  of  1/64  in  (0.40  mm).  The  guided  plate  provided 
positioning  in  the  z-direction,  and  was  rigidly  bolted  to  the  mill 
table. 

Positioning  of  the  traveling  probe  in  the  longitudinal 
direction  was  accomplished  by  adjustment  of  the  probe  holder 
assembly,  described  below,  or  by  moving  the  entire  mill  table  along 
the  tunnel  I-beams. 

A  piece  of  angle  arm,  guided  by  a  notched  plate  attached  to 
the  guided  plate  support,  extended  from  the  traverse  toward  the 
tunnel  roof.  Circular  section  rod  was  mounted  into  a  block 
attached  to  the  end  of  the  angle  arm,  and  extended  through  a  sxot 
in  one  of  the  plexiglas  viewing  windows  of  the  tunnel  and  into  the 
test  section;  see  Figure  3-12. 

The  slot  in  the  viewing  window  was  covered  nearly  entirely 
by  tape  except  where  the  traverse  arm  passed  through.  This  proce¬ 
dure  minimized  tunnel  cross  flows  produced  by  the  pressure  dif¬ 
ferential  between  the  tunnel  and  the  atmosphere. 

A  probe  holder,  shown  in  Figure  3-13,  was  attached  to  the  rod 
and  held  the  hot  film  probe.  The  probe  holder  assembly  allowed  con¬ 
siderable  adjustment  flexibility  in  probe  positioning  through  the 
use  of  set  screws  and  clamps  to  hold  parts  in  position. 

The  fixed  probe  for  the  two-point  velocity  correlation 
measurements  was  held  in  position  by  a  similar  assembly  rigidly 
attached  to  the  tunnel  floor;  see  Figure  3-14.  The  position  of  the 
fixed  block  holding  the  circular  section  rod  was  changed  by  bolting 
it  at  ^rearranged  positions  on  a  flat  plate  secured  to  the  tunnel 
floor. 

The  shedding  frequencies  of  the  probe  support  apparatus 
were  much  larger  than  the  measured  frequencies  of  interest  and  were 
observed  to  be  predicted  by  the  customary  formula  fd/U-0.,2. 
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The  accuracy  of  the  traveling  probe  position  is  estimated 
to  have  been  ±1.3  am  in  in  both  the  x  and  y  directions  attributed 
to  flexibility  in  the  cantilevered  arms  between  the  traverse  mecha¬ 
nisms  and  the  probe.  This  error  bound  was  determined  by  comparison 
of  visual  observations  of  the  probe  with  the  tunnel  at.  zero  speed 
and  at  test  speed.  The  stiffness  in  the  z-direction  was  larger; 
the  error  is  estimated  to  have  been  ±0.8  mm.  The  accuracy  of  the 
fixed  probe,  attached  to  the  tunnel  floor,  is  estimated  to  have 
been  ±0.8  mn  in  all  directions. 

3.6  Instrumentation 

3.6.1  Velocity  Measurement  System 

The  time  dependent  velocities  and  the  mean  values  of  the 
three  components  of  velocity  in  the  wing-body  Junction  flow  were 
measured  with  hot  film  anemometry.  An  arbitrarily  chosen  pitot- 
static  tube  was  used  for  the  primary  velocity  measurement,  and  was 
also  used  to  measure  the  tunnel  free  stream  speed. 

A  schematic  of  the  velocity  measurement  system  is  3hown  in 
Figure  3-15*  The  two  hot  films  used  for  velocity  measurements  were 
maintained  approximately  200°C  above  the  ambient  temperature  using 
temperature  compensation  probes  for  control  resistances.  The  anemo 
meter  bridge  voltages  were  low  pass  filtered  and  d.c.  suppressed 
before  being  input  to  a  Nicolet  660B  Dual  Channel  FFT  Annalyzer  for 
spectral  correlation  analysis.  The  power  spectral  density, 
coherence,  and  relative  phase  between  the  signals  was  output  from 
the  Nicolet  660B  to  a  Tektronix  plotter  to  produce  hard  copies  of 
the  processed  data. 

3. 6. 1.1  Pitot-Static  Tube  Measurement  of  Free  Stream  Velocity 

The  anemometer  bridge  voltages  were  calibrated  against  the 
wind  tunnel  free  stream  speed  as  determined  by  a  reference  pitot- 
static  tube  positioned  in  the  free  stream.  The  reference  pitot- 
static  tube  was  located  at  PA  (-20,  36,  36  cm)  for  all  the 
measurements.  A  photograph  of  the  tube  is  shown  in  Figure  3-16. 
The  pressures  at  the  pitot-static  tube  were  measured  with  the 
pressure  measurement  system  described  in  Section  3.6.2.  The  tube 
was  commercially  manufactured  and  is  shown  in  Figure  3-17.  The 
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Figure  3-15:  Schematic  of  the  Velocity  Measurement  System 
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Figure  3-16:  Photograph  of  Pitot-Static  Tube 
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nose  was  rounded,  and  the  stem  was  3  ran  circular  tubing.  The  pitot 
pressure  orifice  was  13  stem  diameters  forward  from  the  vertical 
support,  and  the  static  pressure  orifices  were  U  diameters  aft  of 
the  nose  and  9  stem  diameters  forward  from  the  vertical  support. 

A  determination  of  the  tunnel  air  speed  was  raa.de  by- 
subtracting  the  reference  static  pressure  from  the  pitot  pressure, 
via  the  differential  transducer,  to  obtain  the  tunnel  dynamic 
pressure  at  the  tube.  The  customary  conversion  from  pressure  to 
velocity  was  made  after  a  determination  of  air  density  from  a  chart 
using  measurements  of  the  tunnel  air  temperature,  pressure,  and 
relative  humidity.  Those  quantities  were  measured  with  the  stan¬ 
dard  tunnel  instrumentation,  and  are  believed  to  have  been  of  suf¬ 
ficient  accuracy  to  determine  the  density  to  within  0.5  per  cent. 
Although  there  were  no  readily  available  standards  against  which 
the  instrumentation  could  be  compared,  they  were  checked  against  a 
second  set  of  instrumentation  and  found  to  be  in  agreement  within 
the  reading  resolution  (±0.25*F,  *0.511  R.H.,  ±0.005  in.  Hg).  Al¬ 
though  a  common  bias  error  may  have  been  present  for  both  instru¬ 
mentation  sets,  the  environmental  conditions  did  not  change  sub¬ 
stantially  throughout  the  test  so  that  relative  differences  in 
tunnel  speeds  resulting  from  determination  of  the  density  are 
negligible. 

3. 6. 1.2  Numerical  Processing  of  Hot  Film  Output 

The  anemometer  bridge  voltages  were  low  pass  filtered  at 
10K  Hz  to  remove  unwanted  signals  such  as  the  megahertz  range  noise 
from  the  anemometer  feedback  circuit.  The  filtered  signals  were 
input  to  the  Nicolet  660B,  which  was  used  to  obtain  discrete 
numerical  values  of  the  time  dependent  analog  bridge  voltage  out¬ 
put. 

To  maximize  the  A/D  resolution  of  the  Nicolet  660B,  the 
input  signal  was  d.c.  suppressed  to  limit  the  voltage  input  to 
values  between  -1  volt  and  +1  volt  for  the  range  of  velocities 
under  measurement.  With  this  procedure  the  span  for  the  anemometer 
bridge  voltages  was  resolved  into  l»096  parts  producing  an  A/D  reso¬ 
lution  of  0.25  n w.  Examination  of  the  instantaneous  digitized  time 
history  of  the  voltages  produced  by  the  time  dependent  flow  in  the 
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wing-body  Junction  showed  that  the  resolution  interval  was  at  least 
an  order  of  magnitude  smaller  than  the  observed  fluctuations  in  the 
wing-body  Junction  flow. 

The  Nicolet  660B  was  operated  in  the  SIG  ENH  mode  to  otain 
numerical  averages  of  the  input  voltages.  The  voltages  were  digi¬ 
tized  at  a  rate  of  two  samples /second  for  at  least  30  seconds,  the 
tunnel  fluctuation  period.  These  numerical  averages  were  related 
to  the  flow  speed  by  calibration. 

The  hot  films  were  calibrated  by  placing  them  in  the 
assumed  free  stream  flow  about  the  model  at  a  position  in  the  same 
y-z  plane  as  the  measurements  were  to  be  mde.  The  tunnel  speed  was 
varied  and  determined  from  the  reference  pitot-static  probe  located 
at  (-20,  36,  36  cm).  It  is  estimated  that  the  "free  stream"  at 
the  probe  calibration  position  may  have  been  2  percent  larger  than 
at  the  reference  pitot-static  tube  position  because  of  blockage 
effects . 

For  low  frequency  resolution  the  filtered  and  d.c. 
suppressed  anemometer  analog  outputs  were  d.c.  coupled  to  the 
Nicolet  660B.  The  Nicolet  660B  was  operated  in  the  DUEL  IK  mode 
for  the  correlation  and  spectral  analysis  of  the  two-point  measure¬ 
ments.  To  accurately  perform  the  analysis  over  the  frequency  range 
of  interest  (usually  0.5  to  200  Hz),  an  analog  filter  internal  to 
the  Nicolet  660B  was  used  to  provide  at  least  70dB  attenuation  of 
frequencies  that  would  cause  aliasing.  The  filtered  analog  signals 
were  digitized  to  provide  numerical  time  histories  of  the  data. 

The  digitized  time  histories  were  numerically  processed  using  Fast 
Fourier  Transform  techniques  to  produce  spectral  correlation  analy¬ 
ses  of  the  data.  The  graphical  output  from  the  analyses  was 
displayed  nearly  real  time  on  a  CRT  screen,  and  was  hard  copied 
using  a  Tektronix  plotter. 

Both  the  frequency  resolution  and  the  number  of  ensemble 
time  histories  of  data  were  programmable  input  to  the  Nicolet  660B. 
The  resolution  was  chosen  to  be  either  0.25  Hz  or  0.50  Hz,  and  the 
number  of  time  history  ensembles  was  usually  chosen  to  be  400.  The 
choice  for  the  number  of  ensembles  was  governed  ty  statistical 


uncertainty  considerations,  Which  are  discussed  later  in  Section 
3.7.  The  number  of  ensembles  was  occasionally  less  than  400  in 
circumstances  where  the  coherence  of  the  signals  was  observed  to 
be  negligible  for  less  than  400  ensembles.  In  either  case,  the 
data  were  analyzed  with  30  percent  redundancy  to  correct  for  power 
leakage  at  the  spectral  window.  This  procedure  produced  a  sta¬ 
tistically  independent  number  of  ensembles  equal  to  1/2  the  number 
of  analytical  time  history  ensembles. 

The  selection  of  the  frequency  resolution  automatically 
determined  the  sample  rate,  and  hence  the  maximum  resolvable  fre¬ 
quency,  in  the  operation  of  the  Nicolet  660B.  The  maximum  fre¬ 
quency  resolved  was  either  100  or  200  H2  depending  on  the  flow 
conditions. 

It  is  noteworthy  for  additional  experimental  work  that  400 
ensembles  of  data  analyzable  to  100  Hz  in  bandwidth  Increments 
of  0.25  Hz  required  thirteen  minutes  for  acquisition.  A  typical 
correlation  contour  with  the  fixed  probe  on  one  side  of  the  wing 
and  the  traverse  probe  on  the  other  side  required  40  measurements 
for  a  total  test  time  of  10  hours.  It  is  apparent  that  even  a 
small  matrix  of  test  conditions  and  measurement  plane  locations  can 
require  an  impractical  total  test  time  in  some  circumstances. 

The  unfiltered  hot  film  analog  signals  were  also  input  to 
a.c.  coupled  oscillocopes  to  provide  convenient  visual  display  of 
the  signal.  This  practice  provided  assurance  that  the  output  from 
the  anemometers  remained  stable;  the  turbulence  signal  had  a  par¬ 
ticular  “look**  to  it  that  disappeared  whenever  certain  instrumen¬ 
tation  problems  occurred. 

3. 6. 1.3  Three-Component  Mean  Velocity 

A  TSI  Model  1247  Cross  Flow  "X"  Probe  hot  film  array  was 
used  to  measure  the  three  components  of  the  mean  velocity  in  the 
wing-body  junction  flow.  The  hot  film  sensor  diameter  and  length 
were  0.025  mm  and  0.51mm,  respectively.  The  distance  between  the 
sensors  was  also  0.51mm,  and  therefore  the  measurement  volume  was 
approximately  (0.5  mm)^.  The  x-array  film  was  calibrated  by  moving 
the  film  to  the  tunnel  free  stream  and  measuring  the  analog  output 
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for  each  film  as  a  function  of  tunnel  air  speed.  This  procedure 
automatically  accounted  for  the  inclinations  of  the  films  to  the 
flow  and  the  resultant  longitudinal  cooling  effect.  Because  the 
measured  flow  direction  relative  to  the  probe  did  not  deviate 
substantially  from  the  calibration  flow  direction,  the  calibrations 
were  sufficiently  accurate.  The  error  in  the  longitudinal  velocity 
was  negligible,  and  the  error  in  the  cross-velocity  was  approxima¬ 
tely  10  percent  low  in  magnitude. 

The  x-array  hot  film  probe  was  rotated  to  two  successive 
positions  to  make  measurements  of  the  three  velocity  components. 

The  error  in  this  procedure  was  only  a  few  percent  because  the 
longitudinal  velocity  was  much  larger  than  either  the  lateral  or 
the  tranverse  velocity.  In  all,  four  velocities  were  measured  at 
each  position  in  the  flow:  a  nearly  transverse  component,  a  nearly 
lateral  component,  and  a  redundent  pair  of  longitudinal  component 
measurements.  Measurement  of  the  actual  probe  orientation  and  the 
angle  between  the  film  elements  permitted  trigonometric  resolution 
of  the  velocities  into  the  appropriate  coordinate  system  com¬ 
ponents. 

The  angles  between  the  films  was  determined  using  the 
method  reported  by  Shabaka  (1979),  which  essentially  assumed  a 
linearized  form  for  the  hot  film  output  for  velocity  directions 
within  about  10  degrees  of  the  mean  direction.  A  TSI  Model  1125 
Calibrator,  modified  with  fixtures  to  produce  the  required  probe 
orientation,  was  used  to  provide  empirical  input  to  the  calcula¬ 
tion.  The  measurement  set-up  is  shown  in  Figure  3-18. 

The  probe  orientation  to  the  flow  was  determined  by  placing 
it  in  the  tunnel  free  stream  and  visually  aligning  it  parallel  to 
the  (assumed)  mean  flow  direction  nearly  as  possible.  The  lateral 
and  transverse  velocities  were  defined  to  be  zero  at  this  calibra¬ 
tion  position.  The  x-array  hot  film  was  then  calibrated  for 
voltage  output  as  a  function  of  tunnel  air  speed. 

The  probe  was  then  placed  in  the  appropriate  measurement 
position.  The  measured  velocity  directions  were  therefore  relative 
to  the  assumed  velocity  vector  at  the  location  where  the  calibara- 


/VySfv.s 


46 


tion  «as  performed.  The  calibration  and  measurement  positions  were 
alvays  in  the  same  plane  (x  *  constant)  normal  to  the  free  stream 

flow. 

3.6.1. U  Two-Point  Velocity  Correlations 

TSI  Model  1262  Miniature  Probes  with  Sensor  Upstream  hot 
films  were  used  to  measure  the  two-point  velocity  correlations. 

The  cylindrical  sensor  element  was  0.023mm  in  diameter  and  0.31mm 
long.  The  single  element  films  were  aligned  parallel  to  the  wall 
to  be  mostly  sensitive  to  fluctuations  in  the  x-y  plane.  Because 
the  flow  was  mostly  in  the  longitudinal  direction,  the  sensitivity 
of  the  films  was  much  larger  for  the  longitudinal  fluctuations 
(x-direction)  than  for  the  tranverse  (y-direction)  fluctuations. 

3*6. 1.5  Temperature  Compensation  for  the  Hot  Films 

The  hot  films  were  temperature  compensated  with  TSI  Model 
1310  Temperature  Compensating  Probes  with  TSI  Model  1325  Trimming 
Network  resistors.  To  select  values  for  the  two  adjustable  compen¬ 
sation  resistors,  both  the  film  and  the  compensation  probe  were 
placed  In  a  refrigerator  and  cooled  to  near  0°C.  From  measurements 
of  the  film  and  temperature  compensation  probe  resistances  at  a 
temperature  of  approximately  25°C  (room  temperature)  and  at  the 
temperature  near  0°C,  values  of  the  parallel  and  series  resistances 
In  the  temperature  compensation  probe  were  set  to  suppress  the  hot 
film  response  to  ambient  temperature  variations.  A  laboratory-type 
mercury  thermometer  was  used  to  measure  the  temperature  at  which 
the  film  and  temperature  compensation  probe  resistances  were 
measured.  Operation  of  the  hot  film  included  meaurements  of  the 
tunnel  temperature.  This  measurement  was  made  by  placing  a  third 
TSI  Model  1310  Temperature  Compensating  Probe  in  the  flow,  and 
using  it  with  the  constant  current  bridge  of  a  TSI  Model  lOS^A  ane¬ 
mometer  to  produce  an  output  voltage  proportional  to  temperature. 

The  output  voltage  from  this  probe  was  measured  with  a  Fluke  8050A 
Digital  Multimeter  and  recorded  along  with  the  hot  film  velocity 
outputs  to  provide  information  for  temperature  corrections  in  the 
event  that  large  temperature  changes  were  encountered;  none  were. 

The  temperature  measurement  probe  and  the  two  temperature  com- 
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pensating  probes  used  to  control  the  hot  film  temperatures  were 
placed  downstream  from  the  wing  and  to  one  side  of  the  wake.  The 
photograph  In  Figures  3-19  shows  the  probes  installed  in  the  tunnel 

3.o.2  Pressure  Measurement  System 

The  pressure  measurement  system  was  used  to  acquire  the 
mean  wall  pressures  and  the  reference  pitot-static  tube  pressures 
from  which  the  free  stream  velocity  was  determined.  A  schematic  of 
the  pressure  measurement  system  is  shown  in  Figure  3-20.  Two  scan! 
valves  were  used  to  route  the  desired  pressure  to  be  measured  to  a 
Barocel  Model  372D  10  Torr  pressure  transducer,  which  was  electron¬ 
ically  conditioned  by  an  electronic  manometer.  The  output  from 
the  electronic  manometer  was  in  both  digital  and  analog  form.  The 
manufacturers  claimed  accuracy  for  the  Barocel  system  is  an  indi¬ 
cated  pressure  within  0.7%  of  the  applied  pressure,  or  0.001  Torr, 
whichever  is  larger.  The  accuracy  of  the  signal  conditioning  was 
periodically  verified  by  placing  the  Internal  calibration  load  on 
the  system  electronics.  There  was  no  pressure  calibration  of  the 
system  other  than  to  compare  its  output  against  the  response  of  an 
inclined-tube  manometer  to  ascertain  that  there  was  no  gross  inac¬ 
curacy  (relative  to  the  manufacturer's  specifications). 

3. 6. 2.1  Reference  Pressure 

The  reference  pressure  for  the  differential  pressure  trans¬ 
ducer  for  all  pressure  measurements  was  the  tunnel  static  pressure 
measured  by  the  reference  pitot-static  tube. 

3. 6. 2. 2  Wall  Pressure  Measurements 

Mean  static  wall  pressure  measurements  were  made  using  care 
fully  constructed  pressure  taps  in  the  wall.  Special  care  was 
taken  to  produce  clean,  square,  and  flush  taps.  The  pressure  taps 
consisted  of  3/64  in.  l.D  (5/64  in.  O.D.)  steel  tubing  mounted 
through  holes  drilled  in  the  plexlglas  wall.  The  holes  were 
drilled  on  a  numerically  controlled  drilling  machine  to  assure  uni¬ 
form  holes  normal  to  the  plate  and  without  burrs.  However,  in 
manufacture  the  plate  rippled  where  the  holes  were  closely  spaced. 
The  surface  distortion  was  not  measured,  and  was  not  visually  noti¬ 
ceable.  However,  the  uneveness  of  the  surface  could  be  felt  by 
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Schematic  of  the  Pressure  Measurement  System 
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running  a  hand  over  it.  The  surface  distortion  is  not  believed  to 
have  affected  the  measurements  in  the  experiment. 

The  pressure  tubes  were  manufactured  by  hand.  Each  tube 
was  cut  to  the  desired  length  after  first  running  a  solid  wire 
through  it  to  minimize  crimping  at  the  cut.  The  end  to  be  inserted 
into  the  drilled  wall  holes  was  carefully  ground  to  make  the 
finished  end  flat,  free  of  burns,  and  square  to  the  tube.  The 
finished  tubes  were  cemented  into  the  wall  holes.  Tygon  1/16-in. 

ID  tubing  was  used  to  route  the  pressures  to  the  Scanivalve  Model 
48J4GM  pressure  switch  and  manifold.  This  scanivalve  switch  was 
used  to  sequentially  select  the  pressures  to  be  measured;  see 
Figure  3-20. 

A  second  scanivalve  switch  was  used  to  alternately  route 
the  pre-selected  wall  pressure  and  the  reference  pitot  pressure  to 
the  Barocel  pressure  transducer.  The  pressure  measurement  output 
from  the  electronic  monometer  consisted  of  both  a  digitized  display 
and  an  analog  signal.  The  digital  display  output  was  updated  so 
rapidly  that  it  was  difficult  to  read  fluctuating  voltages  by  eye. 
Therefore,  the  analog  output  was  input  to  the  Nicolet  660B  and 
dlgltlally  averaged  to  produce  the  actual  measurement  data.  The 
Nicloet  660B  output  was  calibrated  against  the  electronic  manometer 
digital  output  display  under  conditions  when  the  pressures  were 
steady. 

The  resolution  of  the  Nicolet  660B  was  set  at  0.024  mv.,  or 
0.0024  percent  of  the  range  of  the  pressure  transducer.  This  reso¬ 
lution  was  two  orders  of  magnitude  smaller  than  the  voltage 
corresponding  to  the  smallest  measured  pressures. 

The  analog  signal  was  digitized  at  one  sample/two  seconds 
for  32  seconds.  This  record  length  corresponded  to  the  nominal 
period  of  the  small  amplitude  tunnel  oscillations;  see  Section  3.1. 

3.6. 2. 3  Free  Stream  Velocity  for  the  Hot  Film  Measurements 

A  different  arrangement  was  provided  for  the  measurement  of 
the  output  voltage  from  the  electronic  manometer  for  measurement  of 
the  free  stream  velocity  for  the  hot  film  measurements;  see  Figure 
3-15.  For  those  measurements  the  pressures  on  the  pitot-static 


Cube  were  measured  and  converted  to  free  stream  velocities. 

Because  the  Nicolet  660B  was  dedicated  to  processing  hot  film  data, 
an  averaging  DVM  was  used  to  obtain  the  analog  output  from  the 
electronic  manometer.  The  averaging  period  was  10  seconds. 

3.7  Error  Analysis 

The  dominant  errors  in  this  investigation  are  attributable 
to  two  sources:  1)  probe  positioning  error  and  2)  statistical 
uncertainty  in  the  spectral  correlation  analysis.  Explanations  of 
the  terms  for  which  error  bounds  are  discussed  below  are  found  in 
Chapters  5  and  6.  However,  the  error  analysis  is  presented  in  this 
section  as  part  of  the  description  of  the  experiment. 

The  probe  positioning  error  given  in  Section  3.5  was 
(Ax,  Ay,  Az)  «  (±  1.3,  ±1.3,  ±0.8  nn)  for  the  traverse  probe  and 
(±0.8,  ±0.8,  ±0.8  ran)  for  the  fixed  probe.  As  the  consequence,  the 
error  in  the  positions  of  the  correlation  contours  are  (±2.1,  ±2.1, 
±1.6  mm),  or  (A  x/t,  A  y/6,  A  z/t)  ■  (±0.04,  ±0.04,  ±0.03). 

Because  the  errors  were  determined  from  estimates  of  traverse  rigi¬ 
dity,  the  errors  are  bias-type  in  that  whole  sections  of  the  corre¬ 
lation  contours  may  be  in  error  in  the  same  direction. 

A  review  of  the  measurement  results  shows  that  shifts  in 
the  contours  within  that  range  do  not  produce  dramatic  changes  in 
the  data;  do  not  change  any  of  the  conclusions;  and  do  not  require 
caveats  to  qualify  the  data. 

The  statistical  uncertainty  in  the  data  is  attributed 
almost  entirely  to  precision  error  limited  by  the  number  of 
ensembles  N  of  time  history  data.  The  reference  for  the  formulas 
used  for  the  error  analysis  is  Bendat  and  Plersol  (1980). 

The  relative  uncertainty  in  the  power  spectral  densities  is 
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where  -  200  in  the  present  investigation,  in  which  N  •  400 
ensembles  of  data  were  processed  with  50  percent  redundancy. 

The  relative  uncertainty  for  the  quantity  |  |  ■  (y^)^ 


(see  Section  6.1)  Is  derived  from  the  relative  uncertainty  In  the 
coherence : 
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and  Is  presented  in  the  following  table: 


Interval  in  Which  True  Value  |Sy|  Lies 

Measured  Value 
for  IsyI 

2/3  Probability 

9/10  Probabilit 

0.2 

0.14  -  0.24 

0.04  -  0.28 

0.3 

0.26  -  0.36 

0.20  -  0.40 

0.4 

0.36  -  0.44 

0.30  -  0.48 

The  interpretation  of  the  table  is  that  there  is  a  2/3  (or 
9/10)  probability  that  the  true  value  of  |  |  falls  within  the 

interval  of  values  shown  in  the  appropriate  column  opposite  the 
particular  value  for  the  measured  |  |  . 

The  uncertainty  in  the  relative  phase  is  given  in  the 
following  table  using  the  formula: 


The  table  shows  that  the  true  phase  is  indicated  by  the 
measured  phase  if  it  is  known  in  advance  that  the  phase  must  be 
either  "in-phase"  or  "out-of -phase,"  even  for  small  values  of  |  | 

Hence  observations  of  the  measured  phase  were  a  more  sensitive 
indication  of  spatial  extent  of  the  flow  structure  than  was  the 
actual  value  for  |  Sy  |  . 


CHAPTER  4 


MEAN  FLOW  MEASUREMENT  RESULTS 


4.1  The  Onset  Boundary  Layer 

Mean  values  of  the  velocity  and  of  the  turbulence  intensity 
were  measured  In  the  onset  boundary  layer.  The  measurements  were 
compared  with  classical  flat  plate  data  to  demonstrate  that  the 
onset  boundary  layer  was  similar  and  can  be  considered  a  classical 
flat  plate  turbulent  boundary  layer.  The  mean  velocities  were 
also  measured  to  provide  boundary  value  input  to  numerical  com¬ 
putations  to  be  made  in  the  future. 

4.1.1  Velocity  Profiles 

Total  head  pressure  profile  measurements  were  obtained  at 
the  four  locations  shown  in  Figure  4-1  to  determine  velocity  profi¬ 
les  for  the  onset  boundary  layer.  In  the  figure,  the  position 
P(0t0,0)  is  where  the  wing  leading  edge  intersected  the  wall  when 
the  wing  was  installed;  it  was  not  installed  for  these  measure¬ 
ments.  The  profiles  were  evaluated  to  determine  the  similarity  of 
the  onset  flow  to  textbook  published  data  for  developed  flat  plate 
turbulent  boundary  layers.  The  measured  profiles  were  used  in 
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integrations  to  obtain  values  for  the  momentum  thickness  6  and  the 
displacement  thickness  6*,  and  which  are  tabulated  in  Table  4-1. 

The  values  for  9  were  used  to  calculate  the  momentum 
thickness  Reynolds  number,  Rq  ■  UQ0/v,  which  was  compared  to  the 
minimum  value,  Rq  ■  5000,  considered  necessary  by  Rao  (1971)  for  a 
developed  turbulent  boundary  layer.  As  shown  in  the  table,  the 
higher  speed  UQ  ~  37m/s  produced  a  very  large  value  for  the  Reynolds 
number  indicative  of  a  fully  developed  turbulent  boundary  layer. 

The  data  for  UQ  ~  18m/ s  produces  a  Reynolds  number  below  5000,  yet 
the  flow  parameters  indicate  a  nearly  turbulent  flow.  It  is 
concluded  that  the  higher  speed  turbulent  flow  was  well  beyond  the 
onset  of  full  development.  The  shape  parameters  H  were  computed 
and  found  comparable  to  values  expected  for  smooth  wall  turbulent 
boundary  layers  (H  <  1.5).  Those  computed  values  for  H  are  also 
tabulated  in  the  table. 

The  total  head  profiles  were  converted  to  velocity  profiles 
by  assuming  constant  static  pressure  across  the  boundary  layer. 

The  velocity  profiles  were  matched  to  an  accepted  logarithmic  law 
of  the  wall  velocity  distribution  to  determine  the  shear  stress  at 
the  wall.  The  matching  is  shown  in  Figures  4-2(a)  through  4-9(a). 
The  obtained  shear  stress  was  then  input  to  logarithmic  velocity 
defect  law  calculations  for  comparison  with  published  data.  Those 
results  are  shown  in  Figures  4-2(b)  through  4-9 (b). 

The  purpose  of  the  comparison  was  to  establish  consistency 
with  previously  published  data;  the  defect  law  calculations  depend 
on  the  logarithmic  fit  in  the  above  calculations  but  the  two  methods 
are  empirically  related  in  the  classical  literature.  It  is 
observed  from  the  figures  that  the  compared  data  are  generally  con¬ 
sistent  in  that  they  are  near  the  empirical  line  representing  the 
classical  relationship.  The  error  bias  is  above  the  line  for  UQ  - 
37  m/s  and  below  the  line  for  U0  *  18  m/s.  Some  Improvement  in  the 
consistency  between  law  of  the  wall  and  the  velocity  defect  matches 
could  be  expected  by  an  adjustment  in  the  coefficients  for  the 
"universal''  profiles.  However,  the  purpose  of  the  present  measure¬ 
ment  was  to  demonstrate  a  reasonable  "equilibrium"  turbulent  boun- 
•dary  layer,  and  to  provide  credence  that  there  was  nothing  peculiar 


•v 


b)  velocity  defect  lev 


X 


y  /  * 


Velocity  Profile  Matches  for  P(-2,5,y/<,1.5) 


64 


in  Che  onset  flow  Co  produce  Che  measured  temporal  flow  In  the 
wing-body  junction. 

In  a  second  check  of  consistency  in  the  measured  data, 
shear  stresses  obtained  from  the  computed  H  by  direct  comparison 
with  previously  published  data  were  compared  with  the  values  that 
were  obtained  for  law  of  the  wall  matching.  Both  values  are  shown 
in  the  table.  The  values  obtained  from  the  law  of  the  wall  profile 
matching  were  up  to  25  percent  less  than  those  obtained  directly 
using  the  calculated  value  for  H. 

A  second  shape  parameter,  G,  was  therefore  calculated  using 
White  (1974)  to  evaluate  the  accuracies  of  the  two  values  for  the 
shear  stress.  The  calculated  values  for  G  for  the  present  measure¬ 
ments  are  given  in  Table  4-1.  For  both  methods  of  obtaining  the 
shear  stress,  the  value  for  G  (>6.5)  indicates  an  adverse  pressure 
gradient  in  the  tunnel.  The  smaller  value  for  the  shear  stress 
leads  to  a  large  calculated  value  for  G,  indicating  a  substantial 
adverse  pressure  gradient  in  the  tunnel,  which  is  not  reasonable. 
Although  the  larger  value  for  the  shear  stress  also  Indicates  an 
adverse  pressure  gradient,  it  is  small  with  an  equivalent  change  in 
velocity  equal  to  1  percent /foot. 

The  inconsistency  in  the  two  methods  for  obtaining  the 
shear  stresses  is  observed  relative  to  previously  obtained  flat 
plate  data,  and  does  not  necessarily  Indicate  that  there  is 
something  strange  about  the  onset  boundary  layer.  An  explanation 
may  be  found  in  an  error  in  the  measurement  distance  y  from  the 
wall,  or  in  the  geometric  conditions  producing  the  boundary  layer. 
It  may  be  that  the  converging  tunnel  inlet  and  the  ramp  leading  to 
the  plate  on  which  the  wing  was  mounted  had  effects  on  the  tur¬ 
bulent  boundary  layer  not  present  for  a  true  flat  plate  condition. 

The  data  may  be  made  wholly  consistent  if  either  the  shear 
stress,  calculated  from  the  law  of  the  wall  profile,  is  increased 
or  the  value  for  H,  calculated  from  the  velocity  profile,  is 
decreased.  However,  it  is  not  likely  that  the  calculated  value  for 
H  is  in  substantial  error. 

A  third,  separate  calculation  of  the  shear  stress  was  made 
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from  the  gradient  of  the  momentum  (loss)  thickness,  Cf  *  2  d0/dx, 
which  also  indicated  that  the  shear  stress  obtained  from  the  law 
of  the  wall  profile  is  too  small.  If  this  is  the  case  then  an 
explanation  may  be  that  the  measured  value  for  y  is  larger  than  the 
effective  value  for  y.  Because  the  law  of  the  wall  matching  was 
made  for  y  <_  0.2  in.,  a  small  change  in  the  value  for  y  could  lead 
to  substanital  relative  changes  in  the  value  for  u*. 

The  conclusion  from  the  above  computations  and  comparisons 
is  that  the  turbulent  boundary  layer  for  UQ  -  37  m/s  is  fully 
developed  since  it  nearly  fulfills  the  expected  characteristics 
for  flat  plates.  Specifically,  the  momentum  thickness  Reynolds 
number  is  large  (>5000)  and  the  data  can  be  described  by  both  the 
law  of  the  wall  and  the  velocity  defect  law.  The  shape  factor  H 
has  values  expected  for  a  developed  flat  plate  turbulent  boundary 
layer. 

4.1.2  Turbulence  Intensity 

The  turbulence  intensity  was  obtained  from  the  low  pass 
filtered  data  from  the  x-array  hot  film;  the  analog  signal  was 
filtered  at  200  Hz.  The  turbulence  intensity  profile  is  shown  in 
Figure  4-10.  Also  shown  in  the  figure  is  a  classical  flat  plate 
result  from  Hlnze  (1975).  The  lower  intensity  of  the  present 
measurements  is  attributed  to  the  filtering  of  the  analog  signal. 

It  is  believed  that  the  measurements  are  otherwise  comparable. 

4.2  Flow  in  the  Wing-Body  Junction 

The  mean  wall  pressure  and  flow  velocity  were  measured  at 
selected  positions  in  the  wing-body  junction  flow  of  the  a/b  ■  1.5 
infinite  chord  length  wing.  The  measurements  were  made  to  assess 
the  overall  characteristics  of  the  mean  flow.  Specifically,  the 
wall  pressures  were  measured  to  determine  the  length  and  breadth  of 
the  upstream  pressure  ridge  produced  by  the  stagnation  on  the  wing 
leading  edge,  and  to  verify  the  presence  of  low  pressures  to  the 
side  of  the  wing  produced  by  flow  accelerations.  In  addition,  all 
three  components  of  the  mean  velocity  were  obtained  in  the  plane 
(xQ/t,  y/t,  z/t)  ■  (2.00,  y/t,  z/t),  where  the  secondary  flow 


"horseshoe  vortex  was  well  developed,  to  demonstrate  the  strong 
flow  distortion  that  resulted  from  the  vortical  transport. 

Flow  visualization  studies  using  oil  dot  and  oil  film  tech¬ 
niques  were  also  conducted  to  assess  the  type  and  degree  of  three- 
dimensional  separations  present  both  in  the  wing-body  junction  flow 
and  at  the  tip  of  the  wing. 

The  results  from  the  pressure,  velocity,  and  flow  visuali¬ 
zation  measurements  were  compared  and  found  to  be  mutually  sup¬ 
porting. 

4.2.1  Wall  Pressures 

The  measured  mean  wall  pressures  are  presented  as 
pressure  coefficients 


where  p  is  the  static  pressure  at  the  position  (x/t,  0.0,  z/t)  and 
pt  and  ps  are  the  total  head  and  static  pressures  measured  with  the 
reference  pitot  static  tube. 

The  error  In  Cp  attributable  to  variation  In  the  freestream 
static  pressure  with  position,  x,  was  calculated  to  be  less  than 
±0.01  when  the  boundary  layer  growth  on  the  wind  tunnel  walls  and 
tunnel  blockage  produced  by  the  wing  were  considered.  The  maximum 
error  in  Cp  produced  by  the  Baracel  measurement  system  was  an  addi¬ 
tional  ±  0.0007. 

The  wall  pressures  for  positions  x/t  <  0  upsteam  from  the 
wing  are  shown  In  Figure  4-11.  The  wall  pressure  Is  observed  to  be 
peaked  at  z/t  ■  0.0  along  the  z/t  axis,  with  a  cusp  in  the  direc¬ 
tion  along  the  x-axis.  The  Cp  surface  can  be  visualized  as  a 
pressure  ridge  extending  upstream  in  the  negative  direction  of  the 
x-axis.  The  Cp  contour  along  the  x-axls,  which  corresponds  to  the 
pressure  ridge,  is  shown  in  Figure  4-12.  It  is  seen  that  the  ridge 
extends  at  least  5  wing  thicknesses  upstream. 

The  pressure  coefficients  Cp  for  positions  along  the  side 
of  the  wing  (x/t  >  0)  are  shown  in  Figure  4-13.  The  accelerated 
flow  around  the  wing  leading  edge  produced  the  large  negative  coef- 
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flcients,  associated  with  low  pressures,  shown  In  the  figure,  the 
pressures  rapidly  relaxed  toward  Cp  -  0.0  aB  the  fluid  traveled 
downstream  In  the  direction  of  Increasing  x/t,  but  even  at  x/t  - 
2.0  the  pressures  were  still  negative. 

4.2.2  Three  Component  Velocities 

Measurements  were  made  of  the  three  velocity  components 
(u/UQ,  v/U0,  w/UQ)  In  the  plane  (x/t  *  2.00,  y/6,  z/t  >  0).  The 
measurement  probe  was  an  x-fllm  array,  and  successive  measurements 
were  made  with  the  "x"  In  two  planes  to  obtain  the  three  velocity 
components  (the  longitudinal  component  was  measured  In  both  cases). 
The  cross-stream  velocity  v  +  w  Is  shown  in  Figure  4-14,  and  iso- 
veloclty  contours  for  the  magnitude  of  the  longitudinal  velocity  T? 
are  shown  In  Figure  4-15.  The  significant  feature  for  both  figures 
is  the  evidence  that  the  transport  of  axial  momentum  by  the  cross¬ 
stream  velocities  associated  with  the  horseshoe  vortex  literally 
inverted  the  boundary  layer.  High  longitudinal  momentum  fluid  was 
transported  Into  the  boundary  layer  at  the  same  time  that  low 
longitudinal  momentum  was  transported  out.  The  effect  was  similar 
to  a  jetting  of  the  fluid  down  into  the  wing-body  intersection, 
then  out  along  the  body  (away  from  the  wing)  and  up  into  the 
wall  boundary  layer. 

It  should  be  noted  that  this  view  differs  in  perspective 
from  the  notion  that  "the  vortex"  turns  the  boundary  layer  over. 
The  difference  in  perspective  accounts  for  the  downstream 
deterioration  of  the  well-defined  vortex  line  at  the  wing  leading 
edge.  As  the  vortex  is  swept  around  the  leading  edge,  strong  shear 
and  dlffu8ivity  forces  rapidly  disperse  the  vortex  "core"  so  that 
the  notion  of  a  vortex,  in  contrast  to  distributed  vorticity,  is  no 


longer  appropriate.  Hence  Figures  4-14  and  4-15  are  not  reflective 
of  a  vortex  with  a  core. 

The  results  in  Figures  4-14  and  4-15  are  in  substantial 
agreement  with  the  features  observed  by  Shabaka  (1979)  and  by  Oguz 
(1981). 

4.2.3  Flow  Visualization  of  Separations 

Oil  dot  and  oil  film  flow  visualization  studies  were  per- 
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formed  to  evaluate  the  type  and  degree  of  three-dimensional  flow 
separation  associated  with  the  wing-body  junction  and  junction  wake 
flows,  and  to  make  comparisons  with  previously  published  results. 

In  addition,  visualization  studies  were  made  to  evaluate  the  size 
of  the  separation  bubble  produced  by  the  trip  wire  and  to 
demonstrate  that  the  flow  at  the  tip  of  the  wing  did  not  produce 
spanwise  distortion  of  the  flow  in  the  wing-body  junction. 

4. 2. 3.1  Wing-Body  Junction  Flow  Separations 

The  photographs  in  Figure  4-16  show  results  from  oil  dot 
applications  to  the  wing-body  junction  using  a  viscous  mixture  of 
motor  oil  and  lampblack.  The  oil  dots  were  applied  to  the  wing  and 
the  wall  in  the  leading  edge  region  of  the  a/b  ■  1.5  elliptic  nose 
infinite  chord  length  wing.  Because  the  wall  was  vertical  along 
the  direction  of  gravitation  forces,  the  oil  dot  mixture  was  more 
viscous  than  required  for  a  horizontal  surface  to  minimize  the  flow 
of  the  dots  in  response  to  gravitational  forces.  Of  course,  the 
sensitivity  of  the  oil  dot  flow  to  fluid  stress  was  diminished. 

The  variability  of  the  dot  sizes  in  Figure  4-16  is  accounted  for  by 
the  addition  of  ever  larger  dots  during  the  study.  The  tunnel  was 
stopped  to  allow  addition  of  oil.  This  technique  was  an  attempt  to 
force  the  oil  to  flow  in  regions  of  low  stress,  and  was  not  very 
successful. 

Figure  4-16  (a)  shows  two  important  features:  the  desired 
small  separation  produced  by  the  trip  wire  on  the  wing,  and  the 
equally  desired  absence  of  separation  at  the  wing  "shoulder"  where 
the  elliptic  nose  faired  into  the  parallel  middlebody.  The  trip 
wire  produced  a  small  separation  bubble  about  5  mm  long  behind  the 
wire.  It  is  concluded  that  the  trip  did  not  stimulate  large 
separation  on  the  wing  which  itself  could  be  the  source  for 
unsteady  fluid  motion. 

Figure  4-16(b)  shows  the  oil  dot  pattern  on  the  wall  in  the 
region  of  the  wing  leading  edge.  The  figure  is  actually  a  photo¬ 
graph  of  the  piece  of  plastic  tape  shown  in  Figure  4-16(a)  after 
having  been  removed  from  the  tunnel  and  mounted  on  a  sheet  of 
paper  for  closer  study.  An  outline  of  the  relative  position  of  the 
nose  is  also  shown  on  the  paper.  The  oil  dot  pattern  exhibits  the 


expected  features  associated  with  the  horseshoe  root  vortex.  The 
flow  on  the  wall  at  the  wing  leading  edge  was  initially  out  and 
away  from  the  wing  surface  before  being  swept  around  the  wing 
shoulder.  The  oil  dot  streaks  converge  in  the  region  of  flow 
separation  about  20  nm  from  the  wing  shoulder.  This  result  is 
expected.  However,  there  appears  to  be  a  second  convergence  about 
4  mm  from  the  wing  and  on  the  wall.  The  two  convergence  regions 
are  separated  by  a  region  where  the  dots  did  not  move;  hence  a 
region  of  relatively  low  wall  stress,  it  is  not  obvious  from 
Figure  4-16(b)  alone  whether  the  relatively  low  stresses  were  in  a 
region  of  flow  attachment  or  whether  the  second  set  of  convergence 
lines  were  perhaps  produced  by  flow  accelerations  around  the  wing 
trip  wire. 

A  second  flow  visualization  test  was  made  to  learn  more 
about  the  vortex  flow  at  the  leading  edge  of  the  wing  and  to  better 
define  the  wall  region  where  the  low  stresses  were  observed.  The 
results  are  shown  in  Figure  4-17. 

Figure  4-17 (a)  shows  photographed  resultB  for  the  oil  dots 
on  and  around  the  leading  edge  of  the  wing.  It  is  observed  that 
oil  dots  on  the  wing  near  the  extreme  leading  edge  traveled  toward 
the  wall  rather  than  in  a  downstream  direction.  On  the  other  hand, 
the  direction  of  travel  is  observed  to  have  been  mostly  in  a 
downstream  direction  for  oil  dots  to  the  Immediate  side  of  the 
extreme  leading  edge.  The  downward  force  of  the  fluid  in  the  root 
vortex  at  the  leading  edge  was  relatively  strong,  as  evidenced  by 
the  oil  dot  displacement,  for  at  least  a  distance  of  one-half  the 
boundary  layer  thickness. 

The  photograph  in  Figure  4-17(a)  does  not  clearly  show  the 
oil  dot  displacements  on  the  wall  at  the  side  of  the  wing  nose. 

In  this  case  white  tape  was  not  used,  and  reflections  from  the 
plexiglas  wall  were  superimposed  on  the  oil  dot  images  in  the  pho¬ 
tograph.  The  sketch  in  Figure  4-17(b)  shows  the  oil  dot  displace¬ 
ments,  with  the  arrow  origin  located  at  the  initial  position  of  the 
oil  dot,  the  shaft  the  path  of  the  displacement,  and  the  head  the 
final  position.  The  oil  dots  curved  toward  the  wing  at  x/t  ■  0.50, 
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Figure  4-17:  Additional  Surface  Flow  Visualization  Results 
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as  if  Indicating  an  attachment  line  between  the  wing  and  the  pri¬ 
mary  separation  line.  This  result  is  consistent  with  the  pattern 
shown  in  Figure  4-16(b),  and  discussed  above. 

Although  the  region  of  low  stress,  or  flow  attachment, 
suggests  the  existence  of  a  secondary  vortex  counterrotating  rela¬ 
tive  to  the  primary  vortex,  there  was  no  suggestion  of  a  secondary 
vortex  in  the  cross-stream  velocities  shown  in  Figure  4-14.  Those 
velocities  were  obtained  in  the  plane  (2.0,  y/6,  z/t),  which  was 
downstream  from  the  oil  dot  pattern  shown  in  Figure  4-17(b).  It  is 
possible  that  the  secondary  vortex  was  dispersed  sufficiently  while 
traveling  that  distance  to  not  be  immediately  evident  in  the  velo¬ 
city  measurements. 

On  the  other  hand,  the  curvature  of  the  oil  dots  toward  the 
wing  (at  x/t  ■  0.5)  may  in  part  be  produced  by  the  trip  wire  on  the 
wing.  The  high  velocities  in  the  junction  as  a  result  of  vortex 
transport  probably  fill  in  behind  the  small  separation  bubble  pro¬ 
duced  by  the  wire.  On  contact  with  the  wing,  if  this  is  the  case, 
the  high  velocity  flow  may  have  produced  the  small  spanwise  (toward 
the  tip)  displacement  of  oil  dots  on  the  wing,  and  near  the  wall, 
at  x/t  •  0.75. 

It  is  concluded  from  Figures  4-16  and  4-17  that  the  a/b  > 
1.5  elliptic  nose  wing  produced  a  primary  separation  vortex  of  the 
type  observed  by  previous  investigators,  and  that  there  may  be 
secondary  features  of  the  separated  flow  that  are  not  fully  evident 
from  the  present  study. 

4. 2. 3. 2  Wing-Body  Junction  Wake  Flow  Separations 

Figure  4-18  shows  oil  film  flow  visualizations  of  the 
separated  flow  at  the  trailing  edge  of  the  a/b  >1.5  elliptic  nose 
finite  chord  length  wing.  Figure  4-18(a)  shows  the  oil  film 
separation  line  on  the  wing.  Figure  4-18(b)  is  a  photograph  of  the 
oil  film,  which  was  actually  applied  to  a  piece  of  tape  which  was 
removed  from  the  wing  for  the  photograph.  Figure  4-19  is  a  drawing 
that  shows  the  location  of  the  separation  line  on  the  wing.  The 
significant  feature  of  the  separation  line  is  the  curvature  toward 
the  trailing  edge  inside  the  boundary  layer.  Whereas  the  separa- 


Cion  line  defining  che  separated  flow  region  was  parallel  Co  Che 
trailing  edge  outside  che  boundary  layer,  Che  mean  flow  was  unse¬ 
parated  within  the  inner  half  of  the  boundary  layer.  Apparently 
the  horseshoe  vortex  transports  fluid  with  sufficiently  high  energy 
into  the  junction  to  overcome  the  high  pressures  at  the  trailing 
edge  and  to  prevent  separation  in  this  case. 

4. 2. 3. 3  Wing  Tip  Flow  Separation 

Figure  4-20  shows  oil  dot  flow  visualizations  made  at  the 
wing  tip  to  show  that  it  was  not  associated  with  large  spanwise 
effects  distorting  the  wing-body  junction  flow.  The  measurements 
were  made  with  the  a/b  »  1.5  elliptic  nose  infinite  chord  length 
wing.  As  shown  in  the  figure,  the  oil  dots  represent  flow  up  and 
over  the  tip  that  produced  little  spanwise  displacement  of  the  oil 
dots  except  those  very  near  the  tip;  those  were  displaced  toward 
the  tip.  It  is  concluded  Chat  the  tip,  at  least  in  the  leading 
edge  region,  did  not  produce  flow  distortion  in  the  wing-body  junc- 


CHAPTER  5 


TECHNIQUE  FOR  MEASUREMENT  OF  LARGE  SCALE  TEMPORAL  FLOW 


5.1  Development  of  the  Technique 

It  was  anticipated  that  the  large  scale  temporal  flow  would 
be  a  subtle  feature  superimposed  on  the  mean  flow,  and  that 
sophisticated  techniques  would  be  required  to  measure  Its  presence 
and  extent.  A  tentative  hypothesis  was  made  that  the  temporal  flow 
would  be  associated  with  an  inherent  instability  or  natural  fre¬ 
quency  of  the  flow,  and  therefore  would  be  characterized  by  a 
single  frequency.  However,  the  excitation  or  forcing  function  for 
the  oscillation  was  hypothesized  to  be  associated  with  the  random 
flow  fluctuations  in  the  onset  boundary  layer.  It  was  anticipated 
that  the  temporal  flow  amplitude  might  vary  with  time,  and  that  the 
randomly  changing  flow  would  then  exhibit  a  range  cf  corresponding 
frequencies.  Therefore,  it  was  anticipated  that  the  temporal  flow 
would  be  a  narrow  bandwidth  disturbance  superimposed  on  the  general 
turbulence  spectrum. 

The  measurement  objective  was  first  to  detect  the  presence 
of  the  temporal  flow  structure,  and  second  to  determine  the  spatial 
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and  temporal  extent  of  the  structure.  The  detection  was  first 
attempted  using  a  method  developed  by  Strickland  and  Simpson 
(1975).  It  was  concluded  that  the  method  was  not  applicable  to  the 
present  problem,  at  least  In  Its  current  state  of  development.  A 
new  detection  method  was  developed  using  two-point  velocity  corre¬ 
lations.  This  new  application  of  a  standard  statistical  quantity 
solved  many  measurement  problems  simultaneously.  Before  describing 
the  correlation  technique,  the  application  of  the  Strickland- 
Slmpson  method  will  be  reviewed  both  because  it  is  instructive 
about  the  flow,  and  because  it  merits  further  development  for  use 
with  the  wing-body  junction  flow. 

Strickland  and  Simpson  (1975)  experimentally  related  the 
boundary  layer  characteristic  bursting  frequency  to  the  peak  fre¬ 
quency  in  the  log  normal  distribution  of  the  first  moment  of  the 
wall  shearing  stress  power  spectral  density.  The  reasoning  behind 
this  association  follows  from  the  observation  of  Rao,  et  al. 

(1971)  that  mean  burst  periods  are  distributed  according  to  the  log 
normal  law.  Strickland  and  Simpson  reasoned  that  a  long  time 
averaged  spectrum  of  the  shear  stress  would  be  almost  identical  in 
frequency  distribution  to  a  histogram  compiled  from  burst  frequen¬ 
cies  for  the  same  data.  However,  an  adjustment  had  to  be  made 
to  the  spectrum,  namely  multiplying  it  by  the  frequency,  to  obtain 
the  log  normal  distribution.  The  calculations  are  as  follows: 


PSD  df 


f  PSD 


d(ln  f) 


(/  PSD  df)  (/  PSD  df) 

where  PSD  is  the  shearing  stress  power  spectral  density  as  a  func¬ 
tion  of  the  frequency  f.  The  first  moment  of  the  power  spectral 
density,  f  PSD,  when  plotted  as  a  function  of  In  f  was  almost  iden¬ 
tical  to  a  normal  distribution  curve.  The  curve  was  nearly  iden¬ 
tical  to  the  histogram  for  the  burst  periods,  causing  Strickland  and 
Simpson  to  conclude  that  measured  long  time  spectral  averages  were 
equivalent  to  measured  periods  for  individual  bursts.  This  was  a 
significant  advancement  in  frequency  measurement  technology  because 
it  greatly  simplified  the  process. 
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Simpson  et  al  (1981)  later  applied  the  technique  to  separating 
flows  where  the  characteristic  velocity  fluctuation  as  a  function 
of  distance  from  the  wall  was  related  to  the  outer  flow  variables  U 
and  6,  the  boundary  layer  outer  velocity  and  thickness.  Again  the 
justification  for  this  technique  is  reasoned  by  Rao,  et  al  (1971) 
who  observed  an  inner  and  outer  boundary  layer  coupling  such  that 
burst  frequency  phenomena  extend  across  the  entire  boundary  layer. 
Upstream  from  the  separating  flow  Simpson's  results  produced  values 
for  U/f6  nearly  Identical  to  those  compiled  by  Rao  using  histograms. 

In  the  separating  flow  the  values  for  U/f£  were  up  to  five  times 
larger,  presumably  because  the  flow  speed  was  smaller  than  U 
although  the  flow  structure  was  unchanged  in  length  scale. 

The  Strickland-Simpson  method  was  applied  to  the  present 
flow  to  detect  the  presence  of  large-scale  structure.  However,  the 
peak  frequency  was  at  least  an  order  of  magnitude  larger  than 
the  frequency  that  was  eventually  measured  for  the  large  scale  low 
frequency  organized  motion  in  the  wing-body  junction. 

This  initial  analysis  also  demonstrated  that  the  amplitude 
of  the  power  spectral  density  associated  with  the  large  scale 
motion  would  usually  be  so  small  that  it  merely  produced  a  bulge  in 
the  general  background  spectrum  rather  than  a  peak.  Resolution  of 
the  bulge  was  very  difficult  because  the  general  spectrum  in  the 
wing-body  junction  flow  was  grossly  distorted  as  the  root  vortices 
literally  churned  the  boundary  layer.  There  was,  therefore,  no  com¬ 
mon  spectral  distribution  for  the  junction  flow,  and  the  presence 
of  a  bulge  in  the  spectrum  could  not  always  be  detected  or  associated 
with  large-scale  motion.  Conversely,  the  motion  did  not  always 
produce  a  bulge.  Sometimes  it  was  a  filling  in  of  an  inflected 
region.  These  observations  were  not  obvious  at  first,  but  were 
realized  as  the  preliminary  measurement  efforts  were  tediously  and 
laborously  carried  out. 

It  was  eventually  concluded  that  the  large-scale  motion  was 
not  the  turbulent  burst  pheonomenon  as  evidenced  by  the  peak  in  the 
log  normal  distribution  of  the  velocity  and  that  its  frequency 
could  not  be  surely  detected,  or  measured,  with  peaks  in  either  the 
spectrum  or  the  first  moment  of  the  spectrum. 


An  Informal  reexamination  of  the  measurement  objectives 
relative  to  the  nature  of  the  wing-body  junction  flow  was  then 
made.  It  was  considered  that  the  alternate  approach  of  discerning 
the  spatial  extent  of  the  organized  motion  first  might  prove  more 
fruitful  than  the  above  attempt  to  discern  frequency  first.  The 
definition  of  large  scale  structure  Includes  the  notion  of  a  phase 
related  distribution  of  fluid  velocities  in  an  organized  motion 
confined  in  space.  In  theory  the  fluctuating  velocities  in  the 
structure  are  correlated.  Therefore,  in  practice  the  spatial 
extent  of  the  structure  may  be  determined  by  mapping  contours  of 
correlation  values.  Because  the  organized  motion  is  quasl- 
deterministic,  and  therefore  has  characteristics  of  random  motion, 
the  evaluation  of  the  correlation  between  positions  in  the  flow 
must  be  made  using  practical  application  of  statistical  theory. 

A  second  stage  of  preliminary  measurements  produced  fruitful 
results.  Correlations  were  obtained  between  velocities  at  posi¬ 
tions  far  enough  apart  that  turbulence  Integral  length  scales  were 
greatly  exceeded.  This  Implied  that  the  correlations  were  asso¬ 
ciated  with  large  scale  structure  not  inherent  to  the  otherwise 
flat  plate  boundary  layer.  Most  rewarding  was  the  observation  that 
the  cross-spectral  power  spectral  density  produced  a  narrow  band¬ 
width  spectrum.  This  provided  additional  support  that  the  large- 
scale  structure  was  coherent. 

The  remarkable  result  from  this  two-point  measurement  tech¬ 
nique  was  that  the  large-scale  structure  could  be  discriminated 
from  the  background  turbulence  even  when  there  was  no  bulge  in 
either  of  the  power  spectral  density  distributions  at  the  two 
measurement  positions. 

An  added  benefit  from  this  measurement  technique  was  that 
the  fluctuations  in  the  individual  horseshoe  root  vortex  legs  were 
separable.  This  occured  because  the  large-scale  structure  origi¬ 
nated  in  the  wing-body  junction  flow  produced  antisymmetric 
fluctuations  in  the  vortex  legs  so  that  the  relative  phase  was 
nearly  180  degrees.  The  meaning  of  this  observation  will  now  be 
discussed  in  detail. 
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Consider  fluctuations  on  one  side  of  the  nominal  plane  of 
symmetry  of  the  flow,  as  defined  by  the  plane  of  symmetry  of  the 
wing,  relative  to  fluctuations  on  the  other  side.  It  would  be 
expected,  for  example,  that  temporal  oscillations  of  the  separating 
flow  at  the  horseshoe  vortex  origin  would  be  either  fore  and  aft  or 
side  to  side.  The  former  would  produce  fluctuations  in-phase  for 
downstream  positions  symmetric  to  the  plane  of  symmetry.  The 
meaning  of  symmetric  here  is  in  the  sense  of  a  mirror  image.  In- 
phase  fluctuations  imply  invariance  of  the  flow  to  reflection 
through  the  plane  of  symmetry.  Slde-to-side  oscillations  of  the 
flow  at  the  origin  would  also  produce  symmetric  flow,  but  the  phase 
relationship  between  the  velocities  would  be  180  degrees. 

If  the  frequency  of  the  motion  is  sufficiently  low,  then 
moderate  displacements  of  the  nominally  symmetric  measurement  posi¬ 
tions,  or  moderate  asymmetry  of  the  flow,  would  not  be  important  to 
the  relative  phase,  changing  it  by  a  small  value  from  its  basically 
in-phase  or  out-of -phase  value. 

In  this  study  the  fluctuations  were  found  to  be  out  of 
phase.  The  added  benefit  from  this  is  that  one  is  able  to 
distinguish  fluctuations  traveling  down  one  side  of  the  wing  from 
fluctuations  traveling  down  the  other  side,  even  after  the  flows 
merge  downstream  from  the  wing,  if  one  fluctuating  structure  domi¬ 
nates  the  other.  This  happens,  for  example,  for  the  wing  at  inci¬ 
dence  when  the  higher  pressure  vortex  temporally  dominates  the 
lower  pressure  vortex,  as  discussed  later. 

The  step-by-step  measurement  process  that  was  developed  is  as 
follows : 

1.  Analog  output  from  measurement  at  position  and  measurement 
at  position  Pg  is  simultaneously  sampled  at  a  predetermined  rate 
for  a  predetermined  length  of  time. 

2.  Harmonic  analysis  producing  magnitude  and  phase  using  a  Fast 

Fourier  Transform  is  performed  on  the  time  record  of  data. 

SA  "  I  C  .  i  C08  *A  +  J  I  SA  i  *ln  ♦a*  where  SA  "  SA  <f > 

S  .  - 1  SB  |  cos  ♦g  ♦  j  |  SB  |  sin  *B,  where  Sg  -  Sg  (f) 


3.  The  power  spectral  density  and  the  cross-power  spectral  density  are 
computed : 

IGaaI  “  SASA*  ’I  Sa!  2 
tGBBl  *  SBSB  “I  SB  I  2 

IGAB^  “  SBSA  “|  SB  (!  Sa|  cos  ^B  “  +  SA  ||  SB  I  sin  ($B  “ 

4.  Steps  1-3  are  repeated  until  predetermined  precision  criteria 
have  been  met.  The  results  are  ensemble  averaged  to  produce  the 
following  quantities: 

Power  spectral  density  gaa-  I  saI  2 

GBB  *  I  SB  I  2 

Cross  power  spectral  density ,  Sg  ||  S*|  amplitude 

(♦B  “  Phase 

Coherence 

tab2.  Gab  gab* 
gaa  gbb 

5.2  Statistical  Description  of  the  Temporal  Flow 

This  section  presents  the  reasoning  underlying  the  develop¬ 
ment  of  the  measurement  process  described  above. 

5.2.1  Historical  Approach:  Correlation  Function 

The  historical  approach  to  measuring  the  correlation  bet¬ 
ween  stochastic  processes  relied  on  the  cross-correlation  function 
R.«  (t)  defined  as  follows 
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I  ,P 

-  /  A(t)  B(t+x)  dt 

P  © 


This  measurement  was  appropriate  for  the  analog  data  analysis 
Instrumentation  available  before  the  invention  of  electronic  data 
processors  and  the  discovery  of  the  Fast  Fourier  Transform.  Never¬ 
theless,  it  is  Instructive  to  proceed  with  this  historical  approach 
in  developing  a  statistical  description  of  the  temporal  flow,  and 
the  transition  to  the  modern  spectral  approach  can  be  easily  mads. 

The  signals  A(t)  and  B(t)  are  the  velocities  at  two  posi¬ 
tions  in  the  flow.  Of  course,  the  velocities  can  be  for  the  same 
position  in  which  case  the  integrand  is  A(t)  A(t+x)  which  gives 
rise  to  the  auto-correlation  function  w  t).  The  velocities  are 
stochastic  processes,  and  are  considered  stationary  for  this  study. 
Although  the  large  scale  temporal  flow  is  best  described  by 
measuring  the  Instantaneous  velocity  at  all  positions  in  the  flow 
simultaneously,  practicality  does  not  allow  this,  and  in  this  study 
the  flow  was  measured  at  only  two  positions  at  one  instant.  Those 
positions  were  variable  throughout  the  flow. 

The  velocities  are  classed  as  a  bivariate  time  series  pair. 
Because  they  are  Interactive  they  are  said  to  arise  on  a  similar 
footing  contrasted  to  a  casually  related  pair  for  which  one 
variable  is  a  function  of  the  other  in  a  true  cause  and  effect 
relationship.  The  objective  for  this  study  was  to  study  the 
interaction  between  the  velocities  and  to  describe  it  in  terms  use¬ 
ful  for  further  analysis. 

Many  assumptions  were  made  about  the  bivariate  process 
which  were  necessary  to  make  headway  with  a  statistical  descrip¬ 
tion.  The  assumptions  made  are  described  in  detail  by  Jenkins  and 
Watts  (1968).  Principally  it  was  assumed  that  the  stochastic  sta¬ 
tionary  process  represented  by  the  bivariate  time  series  A(t)  and 
B(t)  may  be  adequately  described  by  the  lower  moments  of  their  pro¬ 
bability  distributions.  Specifically,  the  mean,  variance,  and 
covariance  were  assumed  sufficient  to  describe  the  process.  In 
practice  the  description  Involves  measured  variations  of  these 
fundamental  probability  distribution  quantities.  For  instance. 


time  lags  are  Introduced  to  measure  the  nonrandom  features  of  the 
processes,  and  Fourier  transforms  are  employed  to  allow  spectral 
analysis  techniques. 

The  present  study  of  large  scale  temporal  flow  relies  on 
the  theoretical  validity  of  the  cross-correlation  function, 
to  describe  the  predictive  association  between  the  series  A(t)  and 
B(t).  The  definition  of  the  cross-correlation  function  is 

rAB 

PAB(t)  -  (RAA(0)  Rbb(O))1/2 


where  and  Rgg( 0)  are  employed  to  normalize  the  correlation 

for  the  different  scales  at  A  and  B  and  where  the  means  for  A  and  B 
are  set  to  zero.  To  prevent  confusion  here  it  is  noted  that  other 
definitions  of  this  function  are  in  use.  For  example,  Bendat  and 
Piersol  (1980)  refer  to  this  function  as  the  correlation  coefficient 
function. 

If  (t)  Is  Identically  zero  for  all  values  of  t,  the 
bivariate  process  is  uncorrelated.  Furthermore,  if  A(t)  and  B(t) 
are  normal  processes,  they  are  mutually  independent.  If  they  are 
not  normal  processes,  they  may  or  may  not  be  Independent.  Fluid 
mechanics  phenomena  are  often  only  approximated  by  a  normal  pro¬ 
cess.  Therefore  it  Is  possible  that  A(t)  and  B(t)  would  be 
mutually  dependent  but  exhibit  zero  correlation  according  to  the 
definition  of  "correlation." 

If  p£g  (t)  is  not  identically  zero  for  all  values  of  t, 
the  bivariate  process  is  correlated.  There  is  at  least  a  linear 
dependence  of  one  variable  on  the  other.  The  meaning  of 
"dependence"  is  that  there  is  an  interaction  between  the  variables, 
not  necessarily  that  there  is  a  causal  relationship. 

Historically  analog  techniques  were  used  to  measure  the 
cross  correlation  function.  These  techniques  were  notoriously 
inaccurate  for  low  frequency  phenomena  because  of  the  long  measur- 
ment  times  required. 


5.2.2  Modern  Approach:  Spectral  Density 


The  introduction  in  modern  times  of  the  electronic  pro¬ 
cessor  and  the  Fast  Fourier  Transform  led  to  the  spectral  analysis 
approach  to  correlation  assessments,  which  was  used  in  this 
investigation.  The  analysis  is  digital  and  nearly  real  time.  The 
accuracy  of  low  frequency  processes,  such  as  those  encountered  in 
this  study,  is  vastly  improved  because  the  spectral  resolution  is 
not  dependent  on  technically  difficult  analog  bandpass  filtering. 

Mathematically  the  cross-correlation  function  and 

cross-spectral  density  G^g(f)  are  Fourier  transform  pairs. 
Therefore,  the  methods  are  in  principle  equivalent. 

5.3  Statistical  Measurements 

A  Nlcolet  Model  660B  Dual  Channel  FFT  Analyzer  was  used  to 
perform  the  statistical  measurements.  The  anlyzer  was  a  complete 
processor.  It  was  used  to  digitally  sample  the  analog  velocity 
signals,  A(t)  and  B(t),  and  to  compute  the  auto-spectral  density 
functions,  the  cross-spectral  density  function,  the  coherence,  and 
the  transfer  function.  These  quantities  provided  the  necessary 
information  to  statistically  describe  the  large-scale  temporal  flow 
and  to  permit  use  of  the  data  in  further  applications. 

The  basic  functions  of  the  Nlcolet  signal  processor  were  to 
digitize  a  predetermined  time  length  of  data,  called  the  record,  at 
a  predetermined  rate,  and  to  perform  a  Fourier  series  analysis 
using  the  nearly  real  time  Fast  Fourier  Transform  (FFT).  The 
results  from  successive  records  were  averaged  to  produce  statisti¬ 
cally  precise  results.  Data  acquisition  was  continuous  with  the 
data  being  simultaneously  processed  because  of  the  speed  of  the 
calculation  process. 

The  Fast  Fourier  Transform  converted  the  signals  A(t)  and  B(t) 
to  the  spectral  information  SA,  Sg,  $A,  and  where  the  S's  and 
$'s  are  the  amplitudes  and  phases  of  the  Fourier  series  terms: 


|  SA|  cos  ♦a  +  j  |  SA)  sin  $A  where,  eg:  A(t) 
I  SR  I  cos  ♦b  +  J  I  SB  I  sin  $B 


The  auto-spectral  density  was  computed  as  an  ensemble  average 
of  the  record  time  histories: 


CAA  ’  SA  S.* 


GBB  "  SB  SB 


The  cross-spectrum  density  was  similarly  calculated;  there  is 
however,  relative  phase  information  in  addition  to  amplitude 
information: 

GAB  *  SB  SA 

"  I  SA  |  SB  I  cos  ^B  “  $A^  +  -J  I  SA  !  SB  I  sln  ^B  “  ^A^ 

The  coherence,  about  which  more  will  be  said  in  the  next  sec¬ 
tion  was  calculated  as  follows: 


w*  GAB  GAB 

AB  gaa  ^b 


and  is  a  function  of  frequency,  both  in  the  numerator  and  the  deno¬ 
minator. 

The  transfer  function,  useful  for  further  application  of 
the  measurements,  was  calculated  as  follows: 


Hab  Gm 

5.4  The  Coherence 

n 

The  coherence,  y ^B  *s  a  funct*on  frequency.  A  discussion 
of  the  meaning  and  usefulness  of  the  coherence  is  in  order  because 
It  is  the  single  most  critical  measurement  result  on  which  this 
study  is  based.  It  is  used  to  assess  the  interaction  between  velo¬ 
cities,  and  subsequently  the  space-time  extent  of  the  large  scale 
temporal  flow. 


mr- 


93 


If  the  signals  A(t)  and  B(t)  are  correlated,  then  G^g  is 
not  identically  zero  because,  by  definition,  R^g(x)  is  not  iden- 
tically  zero.  Therefore,  y  aB  *8  nonzero  for  correlated  signals. 
An  Interpretative  statement  about  the  coherency  will  be  made  after 
demonstration  of  its  resemblance  to  the  squared  correlation  coef¬ 


ficient. 


The  cross-correlation  function  was  defined  previously 


PAB(t) 


RAB<t> _ 

(RAA  (o)  RBB(0))l/2 


For  each  value  of  the  lag  t,  p^g(x)  expresses  the  degree  of  linear 
dependence  between  A  and  B.  A  Fourier  transform  of  the  correlation 
function  produces 


KAB(f)  “  _AB(fj _ 

(RAA(O)  RBB(0))l/2 


which  provides  a  measure  of  the  correlation  between  the  spectral 
components  of  A  and  B.  Because  G^g  contains  phase  information,  the 
value  for  KAB  (f)  is  dependent  on  the  relative  phase  between  the 
signals.  A  measure  of  the  correlation  between  the  signal  magnitu¬ 
des  alone  can  be  made  by  forming  the  complex  conjugate  of  (f), 
and  taking  the  product 


KAB2  (f)  -  KAB(f)  KAB  <f)* 


GAB  gab* 
RAA(O)  RBB(O) 


Now,  R^(0)  and  R^g(0)  are  simply  the  mean  square  amplitudes  of  the 
signals  A  and  B.  If  the  signal  consisted  of  only  the  frequency  f, 
then  RAA<°>  ■  °AA  and  R^g(O)  *  Ggg.  Therefore 


GAB  gAB* 


K2AB(f)  -  -  =  TAB- 

gaa  gbb 


Hence,  y^ab*  "the  coherency  (sic)  plays  the  role  of  a  correlation 


1 1 1  »■: 
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coefficient  defined  at  each  frequency  f"  (Jenkins  and  Watts,  1968, 
p.353).  Note  that  no  statement  regarding  linearity  is  necessary 
here;  the  definition  of  correlation  coefficient  does  not  depend  on 
linearity;  only  certain  interpretations  of  the  coefficient  are 
based  on  linearity  assumptions. 

The  coherence  then  is  a  normalized  measure  of  the  correla¬ 
tion  between  the  spectral  components  of  A  and  B  where  the  biasing 
effects  of  phase  have  been  removed  by  the  product  of  and  its 
conjugate. 

It  will  be  shown  in  the  next  section  that  the  coherence  can 
be  used,  with  appropriate  assumptions,  to  predict  the  spatial 
extent  of  the  large  scale  temporal  flow. 

5.5  Further  Applications  of  the  Statistical  Description 

All  of  the  measurements  performed  in  this  study  are  necessary, 
but  not  rlgourously  sufficient  to  completely  describe  the  interac¬ 
tion  between  A  and  B.  If  A  and  B  are  not  normally  distributed  then 
higher  order  moments  are  necessary  to  provide  the  full  description. 
There  are  two  possible  sources  of  non-normal  fluctuations  in  the 
present  study:  the  nonlinear  relation  between  the  field  point 
velocity  of  a  vortex  and  the  distance  to  the  center  of  the  vortex, 
and  the  anisotropy  of  the  turbulent  boundary  layer.  The  effects  of 
these  non-normal  fluctuations  on  the  measurements  reported  in  this 
study  is  unknown.  But  it  is  likely  that  the  ultimate  description 
of  the  interaction  between  the  signals  A  and  B  can  be  expressed  as 
an  infinite  sum  of  probability  density  function  moments  (Jenkins 
and  Watts,  1968,  p.75),  of  which  the  linear  term  is  the  dominant. 

Therefore,  it  is  proposed  that  the  further  application  of 
the  results  of  this  study  can  safely  assume  linear  dependence  bet¬ 
ween  A  and  B  for  a  first  order  approximation.  In  this  context  the 
results  fall  neatly  into  the  established  theory  and  practice  of 
bivariate  linear  theory,  which  is  further  discussed  in  Jenkins  and 
Watts  (1968). 

5.6  Low  Frequency  Limitations 
5.6.1  Discussion  of  the  Problem 


Experimental  investigation  of  large  scale  temporal  flow  in 
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a  boundary  layer  generally  implies  the  measurement  of  low  frequency 
phenomena.  An  upper  bound  for  the  frequencies  of  interest  can  be 
obtained  using  the  boundary  layer  thickness  and  the  mean  flow  speed 
as  scales  for  length  and  time.  Although  the  upper  bound  was  600 
Hz  for  the  present  investigation,  the  measured  frequencies  of 
interest  were  as  low  as  0.5  Hz.  Until  recently,  practical  tech¬ 
niques  were  not  available  to  make  accurate  correlation  measurements 
for  signals  at  such  low  frequency  (see  Section  5.2). 

Frequency  measurements  in  the  past  were  made  by  Fourier 
transform  of  analog  generated  correlation  functions,  and  were  known 
to  be  in  error  for  low  frequency  signals.  For  at  least  this  reason 
lateral  correlations  (two-point  velocity  correlations  for  separa¬ 
tions  in  the  z-direction)  were  measured  only  to  the  first  zero 
crossing;  see  Favre,  Gaviglio,  and  Dumas  (1958),  Champagne,  Harris, 
and  Corrsin  (1970),  and  Comte-Bellot  and  Corrsin  (1971).  The  nega¬ 
tive  correlations  obtained  for  distances  beyond  the  first  zero 
crossing  in  these  earlier  studies  were  produced  by  instrumentation 
error  as  well  as  by  fluid  dynamic  phenomena.  Figure  5-1  schemati¬ 
cally  demonstrates  this  point.  The  dashed  line  represents  results 
of  questionable  accuracy  for  f  >  fQ.  Also  shown  in  Figure  5-1 
are  the  improved  results  obtained  with  the  state-of-the-art  signal 
processing  used  in  the  present  investigation.  This  processing 
incorporated  the  inventions  of  the  FFT  and  the  semi-conductor  com¬ 
puter  chip,  which  were  not  available  for  the  previous  investigations. 
As  a  result,  nearly  real-time  direct  digital  spectral  decomposition 
of  the  velocity  correlations  were  available  for  the  present 
investigation.  The  low  frequency  cut-off  for  the  instrumentation 
was  substantially  lowered  as  was  the  lower  frequency  bound  for  the 
measurements.  It  was  possible  to  make  accurate  measurements  in  the 
region  f  >  fQ. 

For  the  present  measurements  to  be  credible  it  is  necessary 
to  show  that  the  instrumentation  low  frequency  cut-off  was  smaller 
than  the  low  frequency  end  of  the  narrow  band  disturbance  charac¬ 
terizing  the  large-scale  temporal  flow  structure.  Evidence  that 
this  objective  was  achieved  is,  firstly,  that  the  lowest  frequency 
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of  interest  with  significant  non-zero  value  of  the  coherence  is 
larger  than  the  instrumentation  low  frequency  cut-off,  as  shown  in 
Figure  5-2,  and,  secondly,  the  second  zero  crossing  for  the  lateral 
correlation  function  was  measured,  as  shown  in  Figure  5-3. 

The  second  zero  crossing  was  a  significant  feature  shown  in 
the  schematic  Figure  5-1.  It  is  also  significant  for  another 
reason.  It  indicates  that  the  measured  flow  disturbances  were  not 
physically  bounded  by  the  wind  tunnel  end  walls,  and  hence  repre¬ 
sented  wall  boundary  phenomena  and  not  duct  flow,  as  indicated  by 
Comt e-Bel lot  and  Corrsin  (1971). 

5.7  Exloratory  Results  (6:1  Elliptic  Nose) 

Preliminary  two-point  velocity  correlations  were  made  to 
explore  the  large  scale  time  dependent  flow  in  the  wing-body  junc¬ 
tion.  The  measurements  were  made  to  test  the  following  hypotheses: 

1.  large  scale  time  dependent  flow  exists  in  the  wing-body 
junction  flow 

2.  the  frequency  bandwidth  is  relatively  narrow  and  consists 
of  low  frequencies  (i.e.,  f  <  0.1  U/6) 

3.  the  flow  structure  is  convected  downstream 

4.  two-point  velocity  correlation  techniques  using  the 
spectrally  distributed  coherence  can  Isolate  the  structure 
from  the  ambient  turbulent  background  even  in  cases  where 
the  structure  local  velocities  are  on  the  order  of  the 
ambient  velocities 

All  of  these  statements  proved  to  be  true  in  the  preliminary 
measurements.  Those  measurements  were  made  with  the  a/b  -  6.0 
elliptic  nose  profile  finite  chord  length  wing  (see  section  3.2). 
However,  the  detailed  correlation  measurements  that  were  made  sub¬ 
sequent  to  the  preliminary  measurements  were  made  with  a  blunter 
nose  profile  (a/b  -  1.5)  both  to  separate  the  disturbance  band- 
widths  and  to  provide  a  larger,  and  hence  more  easily  measured, 
spatial  extent  for  the  secondary  flow  (the  root  vortex). 

The  separation  of  bandwidths  was  desirable  for  the 
following  reasons.  The  downstream  disturbance  bandwidth  produced 


by  the  a/b  ■  6.0  nose  profile  was  numerically  a  subset  frequency 
bandwidth  within  the  larger  bandwidth  for  the  onset  boundary  layer 
structure.  On  the  other  hand,  the  a/b  ■  1.5  nose  profile  produced 
lower  frequency  bandwidth  disturbance  in  a  range  of  frequencies 
outside  the  bandwidth  of  the  onset  boundary  layer  structure. 
Although  the  structure  in  both  cases  must  be  considered  the  result 
of  the  wing-body  junction  flow,  the  structure  in  the  latter  case 
appears  to  be  newly  generated  by  the  junction  flow  rather  than  the 
result  of  a  re-distribution  of  energy  in  the  onset  boundary  layer 
structure.  The  newly  generated  structure  was  a  more  appealing 
feature  for  study  in  this  investigation  where  the  existence  of 
wing-body  junction  flow  structure  was  to  be  proven. 

The  spatial  extent  of  the  horseshoe  root  vortex  secondary 
flow  was  larger  for  the  a/b  ■  1.5  profile  wing  than  for  the  a/b  - 
6.0  profile.  Therefore  it  was  presumably  easier  to  obtain  accurate 
spatial  definition  of  the  correlation  contours  for  the  blunter  pro¬ 
file.  In  view  of  the  purpose  of  the  investigation  to  assess  the 
existence  of  flow  structure,  the  decision  to  use  the  blunter  pro¬ 
file  was  proper  because  it  minimized  risk. 

The  exploratory  results  from  the  preliminary  measurements 
are  reported  both  to  show  that  the  orginal  hypotheses  were  true  an- 
to  show,  In  comparison  with  the  results  for  the  detailed  measure¬ 
ments,  that  the  downstream  bandwidth  depends  on  the  wing  geometry. 
Two  of  the  specific  conclusions  from  the  exploratory  results,  that 
the  flow  structure  is  convected  downstream  and  that  it  grows  In 
spatial  extent  far  downstream  from  the  wing,  were  intentionally  not 
re-confirmed  in  the  time-consuming  detailed  measurements  because 
the  data  were  limited  to  measurements  in  planes  near  the  wing.  It 
is  believed  that  the  general  features  individually  observed  for  the 
preliminary  and  detailed  measurements  are  applicable  to  both  nose 
profiles. 

The  Information  presented  in  this  section  provides  a  global 
experimental  description  of  the  temporal  flow  upstream  of  the  wing, 
along  the  wing-body  junction,  and  downstream  of  the  wing.  The 
measurements  were  made  in  the  nominal  boundary  layer  and  were  not 
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the  result  of  wing  (Strouhal)  vortex  shedding  observed  outside  the 
boundary  layer.  Furthermore,  frequency  bandwidth  and  spatial 
extent  of  the  flow  fluctuations  were  different  than  those  of  the 
undisturbed  boundary  layer;  the  time  dependent  flow  structure 
downstream  from  the  trailing  edge  was  absent  when  the  wing  was 
removed  and  measurements  made  at  the  same  positions  as  with  the 
wing  present. 

Wind  tunnel  measurements  of  the  lateral  (z)  two-point  velo¬ 
city  correlations  made  without  the  wing  present,  l.e.  in  the  onset 
boundary  layer,  indicate  regions  of  negative  correlation  that  are 
associated  with  large-scale  structure  inherent  to  the  boundary 
layer.  That  data  was  obtained  after  the  preliminary  study  and  is 
presented  in  section  6.2. 

The  velocity  correlations  presented  in  this  section  as 
caused  by  the  wing-body  junction  flow  are  clearly  distinguishable 
from  the  onset  boundary  layer  because  1)  the  frequency  bandwidth  is 
much  smaller,  although  possibly  a  subset  of  the  onset  boundary 
layer  bandwidth  and  2)  the  spatial  extent  of  the  structure  in  the 
lateral  direction  is  larger. 

The  data  presented  in  this  section  were  obtained  for  a  wing 
incidence  of  2.4  deg.  The  air  speed  in  the  tunnel  was  approxi¬ 
mately  36.6  m/s. 

Voltages  nearly  proportional  to  velocities  were  measured 
simultaneously  at  two  positions  in  the  flow  with  TSI  hot  wire  ane- 
mometry.  A  single  element  was  used  at  each  position  to  measure  the 
magnitude  of  velocity  fluctuations.  The  analog  anemometer  output 
voltages  were  digitized,  correlated,  and  analyzed  by  the  Nlcolet 
660B  Dual  Channel  FFT  Analyzer. 

The  spectra  were  usually  obtained  for  a  bandwidth  from  0  to 
200  Hz.  This  bandwidth  was  chosen  as  a  compromise  between  fre¬ 
quency  resolution  and  sample  time  after  examination  of  preliminary 
measurements  which  showed  no  pertinent  Information  at  higher  fre¬ 
quencies. 

The  data  in  this  section  are  presented  as  displayed  by  the 
Nlcolet  FFT  Analyzer.  That  is,  the  presentation  is  in  the  units  of 
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volts  and  Hz.  The  conversion  from  volts  to  velocity  cannot  be 
readily  made  because  of  the  manner  in  which  the  span  for  the  anemo¬ 
meter  linearized  output  voltage  was  set;  the  probe  was  Inside  the 
boundary  layer  where  the  velocity  was  only  roughly  known.  A 
ballpark  approximation  is  a  sensitivity  of  46  m/volt-sec. 

The  voltages  associated  with  the  large-scale  structure 
indicate  velocity  fluctuations  on  the  order  of  one  percent  of  the 
local  mean  velocity.  Although  this  information  provides  conceptual 
notions  as  to  the  size  of  the  structure,  the  actual  existence  of 
the  structure  is  shown  by  relative  spectra  and  by  the  nondimen- 
slonal  coherence  rather  than  by  absolute  values  of  the  velocity. 

For  the  present,  the  direct  output  voltages  are  sufficient  to 
establish  these  relative  values  in  the  spectra.  In  either  case 
the  coherence  is  independent  of  linear  calibration  coefficients 
converting  the  voltages  to  velocities  because  they  cancel  in  the 
computation.  Therefore,  the  voltages  alone  without  conversion  to 
velocities  are  sufficient  to  determine  the  global  features  of  the 
large  scale  organized  temporal  flow  in  the  wing-wall  junction 
region. 

The  measurements  provide  strong  evidence  for  large  scale 
organized  motion  in  the  junction  flow.  The  evidence  is  developed 
by  presentation  of  leading  edge  flow  structure  followed  by  up¬ 
stream  correlations  indicating  the  convection  of  the  structure 
along  the  pressure  ridge.  Then  the  structure  will  be  followed  as 
it  traveled  around  the  wing,  interacted  with  the  trailing  edge  con¬ 
ditions,  and  finally  dominated  the  downstream  flow.  During  this 
transit  the  bandwidth,  and  perhaps  the  physical  make-up  of  the  flow 
structure,  changed  substantially. 

Initial  description  of  the  large  scale  temporal  flow 
features  was  obtained  by  two-point  velocity  measurements  in  the 
leading  edge  region  of  the  6:1  elliptic  nose.  Figure  3-4  shows  the 
measurement  planes  for  this  region  as  well  as  for  other  regions  to 
be  presented  later.  The  reference  position  (x/t,  y/6,  z/t)  - 
(-0.12,  0.25,  -0.25)  was  to  one  side  of  the  wing  leading  edge  posi¬ 
tion  and  Inside  the  upstream  separation  line  indicated  by  the  scour 
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Figure  5-4 :  Streamvlse  Locations  of  Measurement  Planes  for  Global 
Experimental  Description  of  Temporal  Flow 


line  found  during  oil  film  flow  visualizations.  On  this  basis,  it 
is  believed  that  the  reference  probe  was  inside  the  root  vortex.  A 
reference  velocity  measured  at  PA  (-0.12,  0.25,  -0.25)  was  corre¬ 
lated  with  velocity  measurements  from  a  second  probe  sequentially 
positioned  on  the  other  side  of  the  plane  of  symmetry  for  the  flow 
(at  positive  values  of  z/t),  then  further  from  the  wall  (at  larger 
values  of  y/5),  and  finally  upstream  from  the  reference  position 
(at  more  negative  values  of  x/t). 

5.7.1  Leading  Edge  Correlations 

Figure  5-5  shows  auto-spectral  densities,  GAA  and  Gg» , 
o 

coherence,  Y^B"  *  relative  phase  *or  a  correlati°n  posi¬ 

tion  Pjj"  near  the  reference  position,  PA  (see  Figure  5-5(a));  a 
correlation  position  PB"  in  the  symmetry  plane  for  the  wing-body 
configuration,  (see  Figure  5-5(b));  and  a  position  Pg  symmetrical 
(about  the  plane  of  symmetry)  with  the  reference  position  PA  (see 
Figure  5-5(c)).  This  series  of  measurements  was  intended  to  show  the 
effects  of  moving  the  probe  away  from  the  reference  position  PA  and 
across  the  plane  of  symmetry  for  the  flow. 

It  is  clearly  seen  in  Figure  5-5(c)  that  the  velocities  at 
the  symmetrical  positions  PA  and  Pg  were  coherent  from  2  to  200  Hz 
and  were  out-of-phase.  That  is,  the  temporal  flow  was  antisym¬ 
metric  although  the  wing-body  configuration  geometry  was  symmetric. 
The  two  other  graphs  in  Figures  5-5(a)  and  5-5(b)  show  that  the 
velocities  became  in-phase  as  the  probe  was  moved  across  the  plane 
of  symmetry  and  near  to  the  reference  position  PA.  This  was 
expected.  Note  that  in  no  case  were  the  power  spectral  densities 
peaked  for  the  coherent  bandwidth;  the  coherent  velocities  were 
apparently  small  in  amplitude. 

When  the  probes  were  very  near  each  other  the  coherence  had  a 
peak  at  a  frequency  of  approximately  50  Hz  (Figure  5-5(a)).  For 
correlations  on  opposite  sides  of  the  plane  of  symmetry  the  peak 
coherence  occurred  at  a  frequency  of  about  12  Hz  (Figure  5-5(c)). 

It  appears,  then,  that  coherence  at  higher  frequencies  decreased 
faster  with  increasing  probe  separation  distance  than  did  the 
coherence  for  the  lower  frequencies.  This  was  not  unexpected 
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because  lower  frequency  implies  larger  spatial  structure.  It  will 
be  shown  later  that  the  lower  frequency  structure  was  charac¬ 
teristically  of  large  spatial  extent  relative  to  the  high  frequency 
structure  throughout  the  wing-body  junction  flow. 

The  extent  of  the  large-scale  structure  across  the  boundary 
layer  normal  to  the  wall  was  evaluated  by  placing  a  probe  at  two 
positions  and  in  the  y-direction  and  correlating  the 
measurements  with  the  velocity  at  P^.  The  arrangement  is  shown  in 
Figure  5-6.  The  reference  probe  was  positioned  at  P^  (-0.12,  0.25, 
-0.25),  and  the  second  probe  was  placed  consecutively  at  Pc  (-0.12, 
0.25,  -0.20)  and  ?c'  (-0.12,  0.50,  -0.20).  All  three  probes  were 
inside  the  boundary  layer. 

Figure  5-6  (a)  shows  the  results  for  PC  atu*  ^A*  ver^  strong 
coherence  and  in-phase  signals  are  indicated.  This  is  expected; 

the  measurement  positions  were  only  slightly  separated  in  the  z- 
directlon.  The  2:1  ratio  of  the  auto-spectral  densities  may  be 
attributed  partly  to  the  proximity  of  the  reference  probe  to  the 
flow  separation,  whereas  the  second  probe  is  near  the  quiescent 
vortex  core,  and  partly  to  the  fact  that  no  attempt  was  made  to 
electronically  force  the  wire  sensitivities  to  be  identical.  . 

Figure  5-6(b)  shows  that  a  relocation  of  the  second  probe 
to  P*  toward  the  outer  boundary  layer  region  greatly  reduced,  but 
did  not  remove,  the  coherence,  and  produced  an  increase  in  the 
phase  relative  to  the  reference  probe  output,  The  phase 

shift  is  indicative  of  larger  velocities  at  the  second  probe. 

The  conclusion  is  that  at  PA  (-0.12,  0.25,  -0.25)  the  ver¬ 
tical  extent  of  the  structure  was  on  the  order  of  the  boundary 
layer  thickness.  It  was  already  shown  that  the  transverse  extent 
was  at  least  one-half  the  boundary  layer  thickness.  Therefore,  the 
large-scale  structure  near  the  wing  leading  edge  was  apparently  as 
wide  as  it  was  thick,  and  it  permeated  the  entire  boundary  layer. 

5.7.2  Upstream  Correlations 

Velocity  measurements  were  made  at  positions  upstream  of 
the  reference  probe  to  obtain  correlations  with  velocities  measured 
with  the  reference  probe  maintained  at  (-0.12,  0.25,  -0.25). 
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The  data  were  used  to  assess  the  streamwise  extent  of  the  temporal 
flow  structure  and  to  determine  whether  the  fluctuations  originated 
at  the  wing  leading  edge  separation  region,  producing  subsequent 
upstream  flow  reactions,  or  If  they  originated  upstream  and  were 
convected  downstream  as  organized  motion. 

The  convection  directions  of  flow  events  were  determined  by 
observing  the  measured  change  in  the  relative  phase  as  the  probe 
was  moved  upstream.  A  convected  disturbance  can  be  described  at  a 
fixed  point  "x"  by  the  function  cos  (2  v  f  (t-tQ)  -2  it  fx/Uc)  where 
Uc  is  the  convection  velocity  In  the  direction  of  increasing  x  and 
$“(-2irfx/U0)  is  the  measured  phase.  The  correlation  between  the 
moving  probe  and  the  reference  probe,  as  measured  in  the  actual 
experiment,  produced  a  relative  phase  ($ref  -  4>m0ving^  8hould 

have  decreased  as  the  probe  was  moved  upstream  for  convected 
disturbances:  (*ref  -  <|>moving)  -  (*A  "  <j>B)  -  (“2  it  fAx/Uc)  where 
Ax  -  Kfgf  “*moving  an^  increaae£*  88  probe  was  moved  upstream. 

In  some  cases  small  values  of  Ax  made  the  relative  phase 
too  small  to  be  accurately  measured.  However  the  effect  of 
Increasing  f  within  the  disturbance  bandwidth  in  those  cases 
revealed  the  sign  of  the  phase  (-2  ir  fA  x/Uc),  and  was  useful  for 
determining  the  convection  direction  of  the  event.  The  quantity 
(2irAx/Uc)  is  the  slope  of  ($ref  -  tm0ving^  aa  a  function  of  f.  In 
the  case  of  Ax>o,  with  the  moving  probe  upstream  from  the  fixed 
probe,  a  negative  slope  is  interpreted  to  mean  a  downstream  convec¬ 
tion. 

If  the  upsteam  fluctuations  are  actually  reactions  to  the 
changing  structure  at  the  leading  edge,  then  the  upsteam  event 
should  have  lagged  the  change  at  the  leading  edge.  In  this  case  no 
convection  of  organized  motion,  in  the  Lagrangian  sense,  would  be 
implied. 

Figure  5-7  shows  the  effects  of  moving  a  second  probe  up¬ 
stream  relative  to  the  reference  probe  fixed  at  (-0.12,  0.25, 
-0.25).  The  figure  shows  a  decrease  in  both  the  coherence  and  the 
relative  phase  as  the  second  probe  was  moved  from  x  »  -1.0  to  x  ■ 
-6.0.  This  result  indicates  that  disturbances  were  convected 
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downstream  (as  explained  above  $ref  ~  ♦moving  “  ♦a  “♦b  decreased  as 
Ax  increased).  The  correlation  distance  was  equivalent  to  at  least 
six  boundary  layer  thicknesses  or  six  wing  thicknesses,  and  is  a 
distance  on  the  order  of  the  distance  over  which  the  pressure 
ridge  produced  by  the  stagnation  pressure  at  the  wing  leading  edge 
extended  upstream.  The  autospectral  density  fell  to  one-third  of 
its  value  at  x/t  -  -0.25  over  the  sane  distance. 

Although  the  pressure  ridge  extended  well  upstream,  it  was 
not  the  source  for  the  out-of-phase  flow  structure  upstream  from 
the  wing  leading  edge  in  this  case.  Measurements  made  subsequent 
to  this  preliminary  study  indicated  that  the  far  upstream  velocity 
correlations  measured  here  were  actually  inherent  to  the  boundary 
layer  (see  Section  6.2).  The  pressure  ridge  provided  a  physical 
reference  separating  out-of-phase  structures,  however.  This  will 
be  discussed  later. 

Figure  5-8  shows  results  for  the  moving  probe  it  z  «  +0.25, 
similar  to  those  obtained  at  z  -  -0.25,  after  an  adjustment  is  made 
to  account  for  the  180  deg  phase  shift  that  occurred  as  the  probe 
was  moved  from  one  side  of  the  plane  of  symmetry  to  the  other. 

(The  phase  shift  occurred  because  of  the  flow  antisymmetry  described 
above.) 

It  is  concluded  that  the  large  scale  temporal  flow  originated 
well  upstream  of  the  leading  edge  region  and  was  perhaps  associated 
with  a  wave  that  traveled  along  the  pressure  ridge.  The  distur¬ 
bance  bandwidth  was  2-200  Hz  and  was  antisymmetric  about  the  plane 
of  symmetry  for  the  specific  correlation  positions  investigated. 

5.7.3  Mid-Wing  Correlations 

The  disturbance's  downstream  progress  was  tracked  by  placing 
one  probe  at  x/t  -  -0.12  and  the  other  at  x/t  -  3.0,  and  then 
correlating  the  measured  outputs.  The  probe  at  x/t  -  3.0  was 
located  where  the  wing  profile  changed  from  elliptic  to  circular 
arc;  it  was  also  the  location  of  the  wing  maximum  thickness. 
Therefore,  the  correlations  were  roughly  between  leading  edge  and 
mid-wing  flow  velocities.  The  results  are  shown  in  Figure  5-9, 
which  shows  that  velocities  at  (-0.12,  0.25,  -0.25)  were 
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strongly  correlated  with  those  downstream  at  Pg  (3.0,  0.25,-0.12). 
Furthermore,  the  phase  for  the  velocities  in  the  mid-wing  region 
relative  to  the  phase  for  the  velocities  in  the  leading  edge  region 
decreased  as  the  frequency  increased.  As  discussed  earlier,  this 
negative  relative  phase  implies  that  events  occurred  at  the  leading 
edge  before  they  occurred  at  mid-wing.  (Note  that  in  this  case  the 
relative  phase  is  ($moving  ”  $ref^*  *s  conclu<*ed  that  there  is 
a  disturbance  that  is  traveling  downstream  from  the  leading  edge  to 
the  mid-wing.  Furthermore,  that  distrubance  has  most  of  its  energy 
in  a  bandwidth  at  frequencies  less  than  200  Hz. 

To  prevent  confusion  it  should  be  mentioned  here  that  the 
power  spectral  density  in  Figure  5-9  is  smaller  than  that  for  Figure 
5-2  for  the  same  position  PA  (-.12,  0.25,  -0.25).  This  is  because 
different  probes  with  slightly  different  flow  orientations  and 
substantially  different  electronic  sensitivities  were  used  for  the 
two  measurements.  This  should  not  have  affected  the  phase  and 
coherence  measurements,  nor  the  relative  distribution  of  the 
autospectral  density  as  a  function  of  frequency. 

Additional  correlation  measurements  were  made  in  the  plane 
x/t  ■  3.0  (at  the  mid-wing)  to  assess  the  spatial  extent  of  the  tem¬ 
poral  organized  flow.  Specifically,  a  systematic  plane  of  symmetry 
correlation  study  was  conducted,  and  the  results  indicated  that 
two-point  mirror-image  positions  P±(3.0,  y/5,  ±j  z/t  |  )  produced 
correlated  velocities  only  for  |  z/t  |  <  0.5.  Figure  5-10  shows  the 
measurement  results  for  P±  (3.0,  0.12,  ±  0.5)  and  P±  (3.0,  0.25,  ± 
0.5).  For  these  positions,  and  for  all  positions  with  |  z/t  |  >  0.5 

the  coherence  between  the  probe  outputs  was  essentially  zero. 
However,  these  observations  apply  strictly  to  mirror-image  posi¬ 
tions. 

Figure  5-11  shows  correlation  results  for  |  z/t  |  <  0.5  and 

I  7/*  |  S.  0»5  for  mirror-image  positions  across  the  plane  of  sym¬ 
metry.  Figure  5-ll(a)  shows  moderate  coherence  for  the  positions 
P±  (3.0,  0.25,  ±  0.25),  which  were  mirror  image  positions  closer  to 
the  wing-body  junction  than  the  positions  (z/t  -  ±  0.5)  described 
above,  for  which  the  velocities  were  essentially  uncorrelated. 
Subsequent  repositioning  of  the  probes  even  closer  to  the  wing  pro- 
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Figure  5-11:  Mid-Wing  Plane-of -Symmetry  Correlations  for  |z/t|  <0.5 
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duced  slgnf leant  changes  In  the  correlation  features  as  shown  In 
Figure  5-ll(b).  The  figure  shows  large  coherence  for  the  positions 
P±  (3.0,  0.25,  ±  0.12),  which  were  very  near  to  the  wing.  The 
relative  phase  between  the  velocity  on  the  two  sides  of  the  wing 
was  180  deg. 

Figure  5-ll(c)  shows  significant  changes  in  the  correlation 
as  the  probes  were  maintained  close  to  the  wing  (z/t  *  ±0.12)  but 
were  positioned  closer  to  the  wall.  The  coherence  was  greatly 
diminished,  but  there  appear  to  have  been  two  bandwidths  of  corre¬ 
lated  flow:  an  out-of-phase  range  from  2  to  4  Hz,  and  an  in-phase 
range  from  20  to  100  Hz.  Either  range  is  considered  "low 
frequency"  in  the  context  of  this  study.  But  for  identification 
purposes  they  will  be  referred  to  in  this  preliminary  study  as  the 
low  and  high  frequency  bandwidths. 

The  correlated  velocities  had  different  spatial  extent  for 
the  two  bandwidths,  at  least  for  positions  near  the  wing  but  rela¬ 
tively  far  from  the  wall.  The  correlations  further  from  the  wall, 
but  at  positions  maintained  close  to  the  wing,  were  significant  for 
the  low  frequency  bandwidth  (2-7  Hz),  as  shown  in  Figure  5-11 (d) 
for  P±  (3.00,  0.5,  ±0.12).  On  the  other  hand,  the  high  frequency 
bandwidth  (20-100  Hz)  was  uncorrelated. 

It  appears,  then,  that  the  correlated  flow  for  mirror-image 
positions  at  x  -  3  is  bound  by  |  z/t  |  <0.5  and  |  y/6  |  0.12,  and  for 
the  high  frequency  structure  (f  >  20  Hz),  the  additional  limitation 
|  y/6  |  <0.5  is  imposed. 

It  is  significant  that  the  disturbance  bandwidth  for  the 
wing-body  junction  flow  was  different  than  that  for  the  far 
upstream  flow.  Either  new  disturbances  were  created  or  only  a  sub¬ 
band  of  onset  boundary  layer  frequencies  passed  through  the  junc¬ 
tion  while  the  remainder  were  suppressed. 

The  spatial  extent  of  the  flow  temporal  structure  was  further 
investigated  by  placing  a  reference  probe  at  PA  (3.0,  0.25,  -0.12) 
and  moving  a  second  probe  to  positions  Pg  (3.0,  y/6,  0.12)  on  the 
opposite  side  of  the  wing.  This  correlation  was  different  than  that 
described  above  because  the  probes  in  this  case  (except  for  the 


obvious  exception)  were  not  at  mirror-image  positions. 

Figure  5-12  shows  the  results  for  positions  varying  from 
y/6  ■  0.3  to  y/6  ■  2.5.  As  y/6  increased,  the  high  frequency  (f  >  20 
Hz)  coherence  rapidly  approached  zero.  The  low  frequency  coherence 
centered  on  f  ■  4  Hz  as  y/6  increased,  and,  although  diminished,  was 
not  zero  for  y/6  -  2.5.  This  position  was  greater  than  two  nominal 
boundary  layer  thicknesses  from  the  wall.  It  appears,  then,  that 
the  low  frequency  temporal  flow,  that  traveled  along  the  pressure 
ridge  upstream  from  the  wing,  had  effects  on  the  wing  flow  well 
outside  the  nominal  boundary  layer.  It  is  not  clear  if  this  effect 
was  produced  by  intermediate  fluctuations  of  the  wing  boundary 
layer,  or  directly  by  potential  flow  reactions  to  the  temporal 
wing-body  junction  flow. 

It  is  reasonable  to  conjecture,  because  of  the  mean  root 
vortex  transport  that  the  flow  in  the  outer  boundary  layer  would  be 
correlated  with  the  inner  boundary  layer  flow  on  the  other  side  of 
the  wing,  although  not  correlated  with  the  outer  boundary  layer 
flow  on  the  other  side.  Figure  5-12(c)  shows  this  strong  "inverted 
boundary  layer"  correlation.  The  coherence  between  velocities  at 
y/6  -  0.5  and  y/6  «  0.25  on  opposite  sides  of  the  wing  was  larger 
than  for  the  mirror-image  positions  y/6  »  ±  0.5  (Figure  5-ll(d))  and 
it  was  comparable  to  the  coherence  for  y/6  -  ±  0.25  (Figure  5- 11(b)) 
for  the  low  frequency  bandwidth. 

The  high  frequency  bandwidth  was  not  correlated  for  the 
"inverted  boundary  layer"  correlation,  although  it  was  for  the 
mirror  image  correlation  y/6  -  ±  0.25.  Therefore,  the  inverted 
boundary  layer  correlation  seems  to  apply  only  to  the  low  frequency 
bandwidth. 

The  cross-stream  spatial  extent  of  the  low  frequency  struc¬ 
ture  outside  the  nominal  boundary  layer  beginning  at  Pg  (3.0,  2.5, 
0.12)  and  extending  in  the  z-direction  is  shown  in  Figure  5-13. 

The  correlations  with  PA  (3,0,  0.25,  -0.12)  were  extended  to  PB 
(3.0,  2.5,  2.0)  in  the  z-dlrectlon,  where  the  correlation  seemed  to 
approach  zero.  It  is  concluded  that  the  temporal  flow  structure, 
initiated  upstream  wholly  within  the  boundary  layer,  had  velocity 
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Figure  5-13:  Plane-of-Symmetry  Correlations  Showing  Cross-Stream 
Spatial  Extent  Well  Outside  Nominal  Boundary  Layer 


perturbation  effects  well  outside  the  boundary  layer  in  the  flow 
around  the  wing.  It  is  not  clear  from  the  measurements  if  these 
effects  were  potential  flow  responses  to  the  temporal  flow  in  the 
boundary  layer.  It  is  clear  that  the  velocity  correlations  for 
these  positions  were  zero  when  the  wing  was  removed  and  that  the 
velocity  fluctuations  are  not  produced  by  mechanical  oscillation  of 
the  wing  as  evidenced  by  accelerometer  correlations.  The  certainty 
of  these  statements  is  argued  in  Section  5.8. 

Further  substantiation  of  the  structure's  extent  is  shown 
in  Figure  5-14  for  the  position  Pg  (3.0,  1.5,  2.0)  which  is  closer 
to  the  boundary  layer  edge  than  the  outermost  position  of  Figure 
5-13.  The  coherence  at  this  position  was  also  non-zero. 

The  outer  boundary  layer  correlations  were  further  studied 
by  moving  the  reference  probe  away  from  the  wall  to  the  position 
(3.0,  0.75,  -0.25),  which  is  in  the  outer  part  of  the  boundary  layer, 
and  correlating  its  output  with  that  from  a  probe  placed  at  the 
positions  P(3.0,  y/6,  0.25).  This  procedure  permitted  plane-of- 
symmetry  correlations  with  the  reference  probe  in  the  outer  boun¬ 
dary  layer  to  be  compared  with  those  described  above  with  the 
reference  probe  further  inside  the  inner  boundary  layer.  Figure  5- 
15(a)  shows  a  small  non-zero  correlation  at  low  frequency  for  the 
mirror  image  position.  Figures  5-15(b)  and  5-15(c)  show  an 
increase  in  the  coherence  as  the  probe  was  moved  toward  the  wall. 

This  is  consistent  with  the  above  description  of  the  "inverted 
boundary  layer"  correlation. 

It  appears  that  the  strongest  correlations  across  the  plane 
of  symmetry  for  the  low  frequency  bandwidth  disturbance  are  for 
positions  such  that  one  is  in  the  inner  boundary  layer  on  one  side 
of  the  wing  and  the  other  anywhere  in  the  (inner  or  outer  parts) 
boundary  layer  on  the  other  side  of  the  wing.  The  correlation  bet¬ 
ween  the  outer  boundary  layer  velocities  across  the  plane  of  sym¬ 
metry  are  small.  All  of  the  correlations  show  an  out-of-phase 
relationship  across  the  plane  of  symmetry. 

It  is  here  noted,  in  preparation  for  the  following  presen¬ 
tation  of  the  correlated  flow  downstream  of  the  wing,  that  the 
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Figure  5-14:  Plane-of-Symmetry  Velocity  Correlations  Showing  Non-Zerp 
Coherence  Outside  Nominal  Boundary  Layer  and  Relatively 
Far  From  Wing 
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Figure  5-15:  Plane-of-Symmetry  Correlations  Showing  Strong 
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angle  of  incidence  (a  *  2.4  deg)  of  the  wing  did  have  global 
measurable  effects  on  the  mean  velocity.  Mean  velocity  profiles 
shown  in  Figure  5-16  for  z/t  >  0  show  a  trough  associated  with  the 
vortex  driven  outflow  of  low  energy  boundary  layer  fluid,  while 
profiles  for  z/t  <  0  do  not  show  a  definite  trough.  The  reason  is 
that  the  vortex  leg  on  the  z/t  <  0  side  of  the  wing  was  spatially 
expanded  by  the  high  pressure  so  that  the  outflow  was  spread  over  a 
larger  area,  producing  a  shallower  trough.  An  investigation  of  the 
effects  of  incidence  is  discussed  in  Section  5.7.5. 

5.7.4  Trailing  Edge  and  Downstream  Correlations 

The  low  frequency  disturbance  was  tracked  in  its  downstream 
movement  by  positioning  the  reference  probe  at  the  mid-wing  position 
P^(3.0,  0.25,  -0.12)  and  the  second  probe  at  a  position  Immediately 
downstream  from  the  trailing  edge,  Pg  (5.12,  0.25,  -0.18).  Figure 
5-17  shows  significant  coherence,  centered  on  the  frequency  4  Hz, 
between  the  velocity  measurements  at  the  two  positions.  The  rela¬ 
tive  phase  is  a  small  positive  angle.  This  result  is  expected 
because  the  low  frequency  (4  Hz)  is  associated  with  large  wave¬ 
lengths  such  that  the  velocities  were  nearly  in-phase  except  that 
the  peak  should  have  occurred  in  time  at  the  upstream  mid-wing 
position  before  it  occurred  at  the  trailing  edge.  Using  a  nominal 
flow  speed,  U  »  30.5  m/s,  a  separation  distance  from  the  mid-wing 
position  to  the  trailing  edge  position,  t  ■  10  cm,  and  a  peak 
coherence  at  a  frequency,  f  -  4  Hz,  the  expected  relative  phase  is 
♦a  ”  ♦b  *  i/(U/f)  360  deg  -  5  deg.  This  is  consistent  with  the 
observed  measurements.  Furthermore,  the  positive  slope  of  the 
relative  phase  as  a  function  of  the  frequency  indicates  downstream 
convection  (in  this  case  Ax  <  o).  It  is  concluded  that  low  fre¬ 
quency  organized  motion  of  the  flow,  most  likely  as  convected  spa¬ 
tial  structure,  passed  from  the  mid-wing  to  the  trailing  edge. 

'i  j>)  point  velocity  correlations  were  then  obtained  to  form 
a  quantitative  measure  of  the  spatial  distribution  of  relative 
phase  within  one  leg  of  the  horseshoe  vortex.  Figure  5-18  is  a 
compilation  of  those  measured  two-point  velocity  correlations, 
which  were  obtained  Immediately  downstream  from  the  wing  trailing 
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Figure  5-16:  Mean  Velocity  Profiles  for  x/t  *  3.0  Showing  Asymmetry 
Produced  By  Wing  at  Small  Incidence  Angle 
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Figure  5-18:  Two-Point  Velocity  Correletione  in  One  Leg  of  Root  Vortex 
at  Trailing  Edge  (x/t  »  5*121 
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edge.  The  reference  probe  was  maintained  at  (5.12,  0.90,  0.50), 

while  the  survey  probe  was  positioned,  as  shown  in  the  figure, 
within  the  same  vortex  leg  as  the  reference  probe.  This  procedure 
was  different  than  that  for  most  of  the  correlations  presented  so 
far,  in  which  the  probes  were  on  opposite  sides  of  the  wing  and 
hence  in  separate  vortex  legs  (although  always  within  the  same  vor¬ 
tex  "tube”  wrapped  around  the  wing  leading  edge). 

Considerable  information  is  presented  in  Figure  5-18  con¬ 
cerning  coherence,  cross-spectral  density,  and  phase.  The  figure 
is  presented  in  two  parts  representing  the  two  frequency  bandwldths 
discussed  previously  for  the  disturbances. 

The  most  relevant  features  of  the  figure  are  the  phase 
boundaries  separating  in-phase  fluctuations  from  out-of-phase  fluc¬ 
tuations.  The  noticeable  difference  between  the  two  frequency 
bandwldths  is  for  the  relative  phase  for  z/t  <  0.5;  the  rela¬ 
tionship  is  out-of-phase  for  f  <  20  Hz  and  in-phase  for  f  >  20  Hz. 
The  determination  of  this  difference  in  fact  was  made  from  the  one 
measurement  shown  in  detail  in  Figure  5-19.  Clearly  the  lower  fre¬ 
quency  bandwidth  had  significant  coherence  with  an  out-of-phase 
relationship  between  the  velocities.  On  the  other  hand,  the 
coherence  for  the  higher  frequency  was  almost  zero  and  the  phase 
was  more  nearly  in-phase  than  out-of-phase. 

A  previous  hint  that  the  two  bandwidth  frequencies  are 
associated  with  different  structures  was  presented  in  Figure 
5-ll(c)  where  they  had  different  phase  relationships  for  plane-of- 
symmetry  correlations  at  mid-wing  (x/t  ■  3.0).  The  lower  frequency 
bandwidth  disturbance  was  out-of-phase  and  the  higher  frequency 
bandwidth  disturbance  was  nearly  in-phase.  Admittedly  these  obser¬ 
vations  Involve  correlations  with  low  values  of  coherence  and  hence 
additional  measurements  are  required  to  confirm  the  phase  boundary 
location. 

Figure  5-18  also  shows,  for  both  bandwidth  disturbances, 
that  the  relative  phase  (♦reference  "  ♦moving*  increased  with  fre¬ 
quency,  except  for  one  instance.  Assuming  that  the  fluctuations  at 
the  trailing  edge  result  from  a  common  upstream  disturbance  even¬ 
tually  passing  through  the  two  measurement  positions,  the  change  in 


relative  phase  with  frequency  must  be  attributable  to  different 
flow  velocities  for  the  two  pathlines  passing  through  the  measure¬ 
ment  positions.  Since,  [($ref  -  $moving)/f]  ~  ('x)/uref  " 
(-x)/Umov£ng)  where  x  is  the  distance  along  the  pathline  from  the 
origin  of  the  disturbance,  the  implication  is  that  the  disturbance 
inside  the  nominal  boundary  layer  had  a  smaller  average  velocity 
than  that  for  the  outer  part  of  the  nominal  boundary  layer,  where 
the  reference  probe  was  located. 

Figure  5-20  adds  credibility  to  the  phase  boundary  pre¬ 
sented  in  Figure  5-18  by  showing  measurement  results  for  two  addi¬ 
tional  positions.  One  position  is  inside  the  phase  ’'trough**  in  the 
z/t  -  0.5  neighborhood  for  f  <  20  Hz  shown  in  Figure  18,  and  the 
other  is  outside  the  trough  (larger  value  for  z/t).  Figure  5-20 
demonstrates  that  the  relative  phase  changed  as  the  phase  boundary 
was  crossed. 

5.7.5  Downstream  Correlations 

Two-point  velocity  correlations  were  measured  in  the  plane 
(18.0,  y/6,  z/t)  located  2.5  chordlengths  downstream  from  the  wing 

trailing  edge.  The  purpose  of  the  measurements  was  to  obtain 
information  about  the  growth  of  the  large  scale  low  frequency  orga¬ 
nized  motion  of  the  flow  originating  in  the  wing-body  junction. 
Specifically,  it  was  determined  that  there  was  dramatic  growth  of 
the  fluctuations  described  above  as  they  were  convected  downstream. 
Far  downstream  from  the  wing,  as  in  the  present  case,  the  flow  was 
a  wing  wake-vort ex-boundary  layer  interaction  problem. 

Figure  5-4  schematically  shows  the  measurement  plane  at 
x/t  ■  18  relative  to  the  wing  position  and  the  free  stream  velocity 
As  shown  in  the  figure,  the  onset  boundary  layer  was  at  small  inci¬ 
dence  to  the  wing.  The  small  incidence  produced  dramatic  changes  in 
the  relative  sizes  of  the  two  vortex  legs  far  downstream  from  the 
trailing  edge.  The  leg  on  the  high  pressure  side  of  the  wing  was 
spatially  enlarged  relative  to  the  leg  on  the  low  pressure  side. 
Some  evidence  for  this  feature  was  shown  in  Figure  5-16,  where  it 
is  shown  that  the  peaks  and  troughs  of  the  velocity  profile  were 
not  symmetric  about  the  wing  at  incidence.  The  accelerated  flow 


Figure  5-20:  Trailing  Edge  Correlations  Showing  Phase 
Increased  Values  of  z/t 
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on  the  low  pressure  side  of  the  wing  stretched  and  intensified  the 
root  vortex  leg  on  that  side,  while  the  root  vortex  leg  on  the  high 
pressure  side  of  the  wing  was  compressed  and  enlarged  in  spatial 
extent.  The  legs  therefore  transported  axial  momentum  Into  the  wing- 
body  junction  with  different  results.  Furthermore,  the  vorticles 
did  not  appear  to  recover,  or  relax,  to  equal  sizes  even  far  down¬ 
stream  from  the  wing  trailing  edge,  as  discussed  below. 

Figure  5-21  shows  the  mean  velocity  profiles  at  x/t  -  18  as 
a  function  of  z/t  for  various  positions  y/6.  The  profiles  show  the 
expected  trough  for  large  y/6  due  to  the  wing  wake  momentum  defi¬ 
cit,  and  the  peak  for  small  values  for  y/6  due  to  root  vortex 
transport  of  high  momentum  fluid  into  the  wing  wake.  The  asymmetry 
of  the  profiles  is  attributed  to  the  differences  in  size  of  the 
root  vortex  legs.  The  larger  vortex,  on  the  high  pressure  side, 
dominated  the  flow.  The  trough  on  the  low  pressure  side  is  pro¬ 
duced  by  the  transport  of  low  momentum  fluid  by  the  spatially 
smaller,  but  more  Intense,  other  vortex  leg.  This  feature  is  dif¬ 
ferent  than  the  trough  at  y  -  4.0,  which  is  produced  by  the  momen¬ 
tum  deficient  wing  wake.  Although  the  larger  vortex  also  displaced 
low  momentum  fluid,  it  was  spread  over  a  much  larger  area  to  pro¬ 
duce  a  shallow  dip  in  the  profile  rather  than  a  distinct  trough. 

Two-point  velocity  correlations  were  measured  in  this  asym¬ 
metric  wake.  The  phase  distributions  between  measurement  positions 
confirm  that  the  temporal  flow  associated  with  the  larger  vortex 
spatially  dominated  the  flow.  Figure  5-22  presents  the  measurement 
positions  for  two  placements  of  the  reference  probe.  The  figure 
also  shows  the  locus  of  positions  for  the  profile  troughs  shown  in 
Figure  5-21.  Because  the  trough  is  postulated  to  have  been  an  effect 
of  low  momentum  fluid  transport  by  the  smaller  vortex  leg,  it  is 
surmised  that  the  smaller  vortex  core  is  to  the  left  of  the  locus. 

The  dashed  line  is  a  possible  representation  of  the  division  between 
large  vortex  dominated  flow,  and  small  vortex  dominated  flow.  Speci¬ 
fically,  in-phase  fluctuations  occur  to  the  left,  and  out-of-phase 
fluctuations  occur  to  the  right,  of  the  dashed  line. 
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All  the  measurements  for  Figure  5-22  exhibited  correlation 
bandwidths  centered  on  approximately  4  Hz  with  one  exception,  those 
for  the  position  Pg  (18*0,  0.5,  0.50,  0.50)  in  Figure  22(b),  which 
was  centered  on  50  Hz.  Figure  5-23  shows  results  for  the  correla¬ 
tions  at  4  Hz  that  support  the  thesis  that  the  large  vortex  spa¬ 
tially  dominated  the  flow. 

Representative  measurements  of  the  coherence  bandwidths  are 
shown  in  Figure  5-23  for  the  reference  position  Pg  (18.0,  0.50, 
-1.50)  shown  in  Figure  5-22  (a).  Results  are  shown  for  the 
out-of-phase  correlations  Pg  (18.0,  0.50,  +  1.50)  and  the  in-phase 
correlations  at  Pg  (18.0,  0.25,  +  0.50). 

As  noted  above,  the  correlation  for  Pg  (18.0,  0.50,  0.50) 
shown  in  Figure  5-22  (b)  indicates  a  correlation  across  the  dashed 
line  for  a  frequency  bandwidth  centered  on  50  Hz.  Figure  5-24 
shows  the  measured  coherence  for  this  correlation.  This  result  is 
consistent  with  measurements  made  for  x/t  -  3,  which  Indicated  two 
non-zero  correlation  bandwidths  for  mirror-image  positions  on 
either  side  of  the  wing.  Figure  5-11  (c)  shows  the  two  bandwidths. 
The  higher  frequency  is  from  20  to  100  Hz  and  is  in-phase  across 
the  wing.  It  was  shown  that  this  higher  frequency  bandwidth  is  not 
spatially  correlated  for  a  very  large  distance.  It  is  likely  that 
this  higher  frequency  motion  is  also  evident  in  Figure  5-22  (b) , 
appearing  as  the  in-phase  correlation  across  the  dashed  line. 

For  contrast,  Figure  5-25  shows  the  correlation  bandwidths 
for  the  position  P^  (18.0,  1.00,  0.50)  which  is  centered  on  the 
lower  frequency  4  Hz  and  is  out-of-phase  with  the  reference 
(18.0,  0.25,  0.50).  The  higher  frequency  correlation  for  this  case 
is  zero,  consistent  with  the  small  spatial  extent  observed  at  x/t  » 
3  for  the  higher  frequency  disturbance. 

Of  special  significance  is  the  very  narrow  frequency  band¬ 
width  in  Figure  5-25.  The  velocity  correlations  were  non-zero 
within  a  2-11  Hz  bandwidth,  which  is  dramatically  smaller  than  the 
2-200  Hz  bandwidth  in  the  onset  boundary  layer.  However,  it  is 
difficult  to  distinguish  the  physics  governing  the  two  frequency 
bandwidths  because  they  overlap  for  this  particular  wing. 
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The  subsequent  detailed  correlation  studies  were  made  with 
the  blunter  a/b  •  1.5  nose  profile,  and  the  frequency  bandwidths 
were  observed  to  be  separated. 

5.7.6  Significant  Conclusions 

Several  significant  conclusions  were  made  from  the  explora¬ 
tory  results,  and  they  provided  the  motivation  for  the  detailed 
study  that  followed.  Those  conclusions  were: 

1.  Large  scale  time  dependent  flow  does  exist  in  the  wing-body 
junction. 

2.  There  are  at  least  two  relatively  narrow  low  frequency  band¬ 
width  disturbances  present,  one  from  the  onset  boundary  layer 
structure  and  another  produced  by  the  presence  of  the  wing. 

3.  The  flow  structure  is  convected  downstream,  and  appears  to  grow 
in  spatial  extent  in  the  wing  wake. 

4.  Two-point  velocity  correlation  techniqes  employing  the 
spectrally  distributed  coherence  appear  to  separate  the  flow  struc¬ 
ture  from  the  ambient  turbulence  even  where  the  structure  veloci¬ 
ties  are  on  the  order  of  the  ambient  values. 

5.8  Tests  to  Verify  Fluid  Dynamic  Origin  of  Temporal  Flow 

The  large  scale  time  dependent  flow  measured  in  the  wing- 
body  junction  was  of  low  level  and  had  subtle  characteristics.  It 
was  prudent  to  perform  several  tests  to  rule  out  instrumentation  or 
facility  related  noise  origins  for  the  measured  temporal  features. 
The  following  tests  were  conducted: 

1.  Sequential  tests  with  and  without  the  wing  present  to 
demonstrate  the  absolute  requirement  that  the  wing  be 
present  to  establish  the  flow. 

2.  Free  stream  velocity  and  acoustic  pressure  correlations 
with  the  boundary  layer  velocities  to  rule  out  possibility 
of  external  flow  driven  phenomena. 

3.  Accelerometer /anemometer  output  correlations  to  rule  out 
mechanical  vibrations  as  causes  for  the  flow. 

4.  Velocity  correlations  for  probes  in  dose  proximity 
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(single  point  two  probe  correlations)  to  rule  out  the 
possibility  of  instrument  hook-up  error,  such  as  an 
improperly  set  inversion  (sign)  switch. 

5.  Variations  in  the  velocity  to  test  if  the  velocity 

correlations  were  functions  of  velocity  to  demonstrate  the 
fundamental  fluid  dynamic  origin  of  the  time-dependent 
flow. 

The  above  tests  were  conducted  during  the  preliminary 
exploratory  measurements  of  the  wing-body  flow.  Consequently,  they 
were  performed  at  different  times  and  with  different  experimental 
set-ups,  such  as  different  nose  shapes  for  the  wing.  However,  it 
is  believed  that  the  results  in  each  case  are  applicable  to  the 
experiment  in  general.  None  of  the  tests  above  produced  evidence 
of  non-fluid  dynamic  sources  for  the  temporal  measurements. 

Example  results  from  the  tests  are  presented  below. 

The  tests  consistently  showed  that  the  large  scale  time 
dependent  flow  in  the  wing-body  Junction  is  an  inherent,  naturally 
occurlng  feature  of  the  flow.  It  is  not  a  product  of  the  par¬ 
ticular  experimental  arrangement. 

5.8.1  Sequential  Measurements  With  and  Without  Wing 

Sequential  velocity  correlations  were  made  with  and  without 
the  wing  to  show  the  presence  of  the  large  scale  time  dependent 
flow  only  in  the  presence  of  the  wing.  An  example  of  the  results 
is  shown  in  Figure  5-26,  which  compares  correlation  data  on  two 
sides  of  the  wing  at  midchord  with  data  at  the  same  measurement 
positions  but  with  the  wing  removed.  Clearly  the  relative  phase 
between  the  correlated  velocity  measurements  is  180  degrees  with 
the  wing  present.  With  the  wing  removed  the  data  are  uncorrelated. 
This  test  showed  that  the  wing  produces  the  large  scale  time- 
dependent  flow.  It  does  not,  however,  discern  if  the  wing  mechani¬ 
cally  amplifies  flow  structure  inherent  to  the  boundary  layer,  or 
if  it  "creates"  the  flow. 

5*8.2  Acoustic  Pressure  and  Free  Stream  Velocity  Fluctuations 

Acoustic  pressure  and  flow  velocity  measurements  in  the 
free  stream  flow  were  correlated  with  hot  wire  anemometer  outputs 
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from  the  undisturbed  boundary  layer.  The  correlations  were  zero, 
indicating  that  the  external  flow  was  not  temporally  forcing  the 
boundary  layer. 

An  example  of  the  measured  acoustic  correlations  is  shown 
in  Figure  5-27.  The  data  were  obtained  with  no  wing  present. 

There  is  no  correlation  between  the  boundary  layer  velocity  fluc¬ 
tuations  and  the  acoustic  pressure.  Two  of  the  peaks  in  the 
acoustic  pressure  spectrum  shown  in  the  figure  are  identified  with 
the  rotation  of  the  downstream  propeller  of  the  wind  tunnel.  The 
data  were  obtained  with  a  nominal  propeller  rotation  speed  equal  to 
25  rps;  the  propeller  had  four  blades.  The  shaft  frequency  of  25 
Hz  and  the  propeller  blade  frequency  of  100  Hz  are  both  clearly 
seen  as  peaks  in  the  pressure  spectrum.  Neither  of  those  peaks  is 
measured  as  driving  the  boundary  layer  velocity  fluctuations.  The 
tests,  therefore,  showed  that  the  acoustic  pressure  in  the  free 
stream  was  not  the  source  for  any  velocity  fluctuations  in  the 
boundary  layer,  and  it  is  reasoned  that  it  is  not  the  source  for 
the  wing-body  junction  time  dependent  flow. 

Additional  measurements  showed  that  the  time  dependent  flow 
is  not  correlated  with  the  very  small  essentially  unmeasurable  free 
stream  velocity  spectrum. 

It  is  concluded  that  the  time  dependent  flow  in  the  wing- 
body  junction  is  created  by  the  interaction  between  the  wing  and 
the  upstream  boundary  layer  and  not  by  free  stream  velocity  or 
acoustic  pressure  fluctuations. 

5.8.3  Mechanical  Vibration  Tests 

Mechanical  vibration  tests  were  performed  by  correlating 
accelerometer  output  with  the  hot  wire  anemometer  output.  The 
accelerometer  was  placed  at  many  locations  throughout  the  experi¬ 
mental  set-up.  In  no  case  was  there  a  correlation  between  the 
accelerometer  and  the  anemometer  outputs.  An  example  measurement 
is  shown  in  Figure  5-28.  The  accelerometer  output  shows  sharp 
spectral  peaks  at  49  Hz  and  at  101  Hz,  with  the  accelerometer 
placed  to  be  sensitive  to  wing  "flapping."  For  the  preliminary  two 
point  velocity  correlation  measurements  one  of  the  frequency  bands 


Figure  5-28:  Example  Mechanical  Vibration  Test  Results  for 
Correlations  of  Structural  Acceleration  With 
Flow  Velocity  in  Wing-Body  Junction 


produced  by  Che  presence  of  Che  wing  (20-100  Hz)  did  have  a  maximum 
coherence  at  approximately  SO  Hz.  However,  the  spatial  region  in 
which  this  frequency  band  was  evident  was  very  small  which  indi¬ 
cated  that  it  was  not  produced  solely  by  wing  oscillations.  For 
the  detailed  measurements,  neither  natural  frequency  was  in  the 
range  of  the  measured  time  dependent  flow  frequencies  produced  by 
the  presence  of  the  wing.  It  is  concluded  that  the  large  scale 
time  dependent  flow  is  not  produced  by  mechanical  oscillation  of 
the  test  arrangement  apparatus. 

5.8.4  Single-Point  Two  Probe  Correlations 

Several  times  throughout  the  course  of  the  investigation 
the  two  hot  wire  probes  used  for  the  velocity  correlations  were 
placed  next  to  each  other  to  verify  that  the  outputs  were  nearly 
identical,  correlated,  and  in-phase.  This  test  showed  that  the 
subtle  features  observed  for  the  large  scale  time  dependent  flow, 
such  as  the  out-of-phase  phenomenon,  were  not  products  of  error 
instrumentation  hook-ups.  Figure  5-29  is  an  example  comparison  for 
two  probes  first  placed  on  opposite  sides  of  then  wing  placed  in 
close  proximity.  The  spectra  are  in-phase  for  the  latter  case 
although  out-of-phase  for  the  former  case.  The  fall-off  of  the 
coherence  for  the  larger  frequencies  with  the  probes  in  close 
proximity  is  expected  as  the  wavelengths  become  less  than  the 
separation  distance  between  the  probes. 

5.8.5  Effects  of  Varying  Velocity 

If  the  flow  fluctuations  were  fluid  dynamic  in  origin  then 
they  should  have  been  functions  of  the  free  stream  velocity.  A 
simple  test  was  made  by  observing  the  peak  frequency  for  the 
coherence  as  a  function  of  speed.  An  example  of  the  results  is 
shown  in  Figure  5-30  for  a  selected  wing.  In  this  particular  case 
there  were  two  peaks  in  the  coherence.  As  shown  in  the  figure, 
both  peak  frequencies  are  nearly  proportional  to  the  free  stream 
velocity. 

5.9  Identification  and  Isolation  of  Trailing  Edge  Shedding  Frequencies 

Measurements  were  made  at  the  trailing  edge  of  the  wing  to 
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Figure  5-30:  Measurements  of  the  Peak  Coherence  Frequency 
as  a  Function  of  Velocity 
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Identify  the  shedding  frequency.  The  measured  frequencies  outside 
the  boundary  layer  were  much  larger  than  those  for  the  large-scale 
structure  investigated  inside  the  boundary  layer.  The  frequencies 
increased  with  increasing  distance  from  the  wall  as  would  be 
expected  for  constant  Strouhal  number  fr/U  as  the  velocity 
increases,  where  r  is  the  effective  radius  of  curvature.  However, 
the  frequency,  or  at  least  the  shedding  phenomena,  is  absent  from 
the  inner  half  of  the  boundary  layer.  Figure  5-31  shows  an  example 
set  of  measurements  obtained  at  the  wing  trailing  edge  for  the  a/b 
■  6  nose  profile.  As  the  distance  from  the  wall  Increased  the 
large  shedding  frequency  first  becomes  noticeable  for  y/{  ■  0.5. 
Outside  the  boundary  layer  the  prominent  shedding  frequency  is  400 
Hz.  Figure  5-32  shows  similar  results  for  the  a/b  ■  1.5  nose  pro¬ 
file. 

It  is  concluded  that  the  measurements  have  isolated  the 
shedding  frequency,  and  that  it  is  not  the  cause  for  the  low  fre¬ 
quency  flow  in  the  wing-body  junction. 


250  -  500  Hz 
(shedding  bandwidth) 


Figure  5-32:  Data  Showing  Shedding  Frequency  at  Trailing 
Edge  of  a/b  •  1.5  Elliptic  Nose  Wing 


CHAPTER  6 


LARGE  SCALE  TEMPORAL  FLOW  MEASUREMENT 
RESULTS  (a/b  -  1.5  ELLIPTIC  NOSE) 


6.1  Method  for  Presentation 

6.1.1  Correlation  Contours 

A  method  was  needed  to  present  the  data  in  a  way  that  the 
temporal  flow  structure  over  a  substantial  portion  of  the  flow 
could  be  depicted  in  a  single  graphical  display.  The  method  that 
evolved  produces  correlation  contours  in  the  (y/t,  z/t)  plane  nor¬ 
mal  to  the  free  stream  flow  direction.  The  contours  are  of  the 
contrived  quantity  St,  which  is  similar  to  the  linear  correlation 
coefficient.  Sy  numerically  decribes  the  correlations  between  the 
velocity  at  a  fixed  probe  P^  and  the  velocity  at  a  traversing  probe 
Pg.  The  quantity  Sy  is  defined  to  be  the  square  root  of  the 
nominal  peak  value  of  the  coherence  multiplied  by  the  cosine  of  the 
relative  phase  between  the  velocities  at  the  two  probe  positions; 

Sy  =  cos  ($b  “  ^  )  >/yab^*  For  the  low  frequency  motions  investi¬ 
gated  in  this  experiment,  cos  )  is  approximately  +1  or  -1 

for  probes  at  the  same  longitudinal  position  in  the  flow. 
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Therefore  Sy  5  iy^g  where  the  sign  is  determined  by  the  relative 
phase  of  the  velocities;  in-phase  implies  +1,  and  out-of-phase 
Implies  -1. 

The  justification  for  equating  y^g  to  the  linear  correla¬ 
tion  coefficient  was  described  in  Section  5-4.  In  that  section 

9 

was  described  as  the  square  of  the  linear  correlation  coef¬ 
ficient  for  each  frequency. 

The  magnitude  of  the  contour  value  for  y^g^  was  obtained  in 
the  following  manner  for  each  pair  of  correlation  measurement  posi¬ 
tions,  PA  (*0/t,  y0/t,  zQ/t)  (fixed  position)  and  Pg  (xQ/t,  y/t, 
z/t)  (variable  position).  The  nominal  peak  value  of  the  coherence 
was  observed  to  occur  within  a  band  of  frequencies  characterizing 
the  disturbance,  as  shown  in  Figure  6-1.  The  fluctuations  of  the 
coherence  with  frequency  about  the  nominal  peak  value  were  taken  as 
manifestations  of  the  statistical  uncertainty  in  the  measurement  of 
the  signal.  Therefore,  a  line  was  drawn  through  the  nominal  mean 
value  for  the  coherence,  and  projected  to  obtain  the  value  y  ^g  for 
calculation  of  S^.  The  sign  for  was  obtained  directly  from  the 
phase  relationship  at  the  nominal  peak  value  of  the  coherence. 

This  procedure  is  illustrated  in  Figure  6-1. 

The  band  of  frequencies  for  the  non-zero  coherence  was  not 
identical  for  all  the  measurements.  In  general  there  were  three 
frequency  bands,  and  a  notation  as  to  which  bandwidth  Sy  applies  is 
made  for  each  presentation  of  results.  The  three  frequency  band- 
widths,  and  the  flow  regions  in  which  they  appear  to  originate  are 
as  follows: 

Bandwidth  I:  2-65  Hz;  onset  boundary  layer  (characteristic 

center  frequency  f0) 

Bandwidth  II:  25-80  Hz;  wing-body  junction  flow  (characteristic 

center  frequency  fgj 

Bandwidth  III:  0.5-10  Hz;  wing-body  junction  wake  flow 

(characteristic  center  frequency  fw^) 

Although  the  numerical  values  of  the  frequency  bandwidths  overlap, 
they  usually  do  not  all  appear  simultaneously  in  the  measurements 
and  are  distinguishable  in  the  data. 


The  contours  for  Sy  were  obtained  by  plotting  the  measured 
data  at  the  y/t  and  z/t  positions  of  the  traversing  probe,  Pfi 
(xQ/t,  y/t,  z/t). 

Caution  should  be  exercised  when  interpreting  the  contours 
of  Sy:  they  apply  specifically  to  correlations  relative  to  veloci¬ 
ties  at  the  fixed  probe  position,  P^  (xQ/t,  yQ/t,  zQ/t),  and  are 
not  necessarily  appropriate  for  other  positions.  Specifically, 
correlations  between  velocities  at  two  arbitrary  positions  of  the 
traversed  probe  cannot  be  deduced  from  correlations  at  each  of 
those  positions  relative  to  the  fixed  probe.  This  will  be 
discussed  in  detail  in  Chapter  7. 

6.1.2  Scaling  Quantities 

This  section  describes  Che  reasoning  behind  a  rational 
choice  for  scaling  parameters.  The  dimensionless  frequency  is 
defined  as  f6/UQ  where  6  is  the  onset  boundary  layer  thickness  and 
U0  is  the  free  stream  velocity;  and  the  dimensionless  velocity 
fluctuation  is  u'/U0.  Correspondingly,  the  power  spectra  density 
♦(f)  -  u'2(f)/df  of  the  velocity  fluctuations  u'  is  made  dimen- 

♦(f)  u'2  (f) 

_  ,  —  _  A/ 

sionless  by  forming  ,  „  =  ♦(f).  It  is  reasonable  to 

SUo  At  A  Uo 

expect  that  the  frequency  band  of  the  large  scale  temporal  flow  in 

the  wake  region  is  physically  governed  by  the  non-dimensional  set 
of  numbers: 

Sfwb 

*  (^/t,  a/b,  c/ t ,  u,  0,  Cf ,  Rd ,  fo6/Uo»^z/6) 

Uo 

where  c  ■  wing  chord  length 

a  ■  wing  incidence  angle 
0  -  Included  angle  of  circular  arc  tail 
Rg  -  boundary  layer  thickness  Reynolds  number 
4_  -  lateral  extent  of  flow  structure  in  onset  boundary  layer 

The  independent  variables  from  the  onset  boundary  layer  are  Cf, 

R*.  f-6/U,.  and  4„/6.  Because  test  data  were  obtained  for  a  fixed 
o  *  O  o 9  z 

tunnel  geometry  and  free  stream  velocity  (30.5  m/s)  those  indepen¬ 
dent  variables  were  fixed  as  follows: 
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Cf  -  0.00225  («  de/dx) 

R6  -  104  (Rq  -  105) 

fQ6/U0  *  0.05  (where  the  center  frequency  fQ  was  chosen  from  band¬ 
width  I) 

A2/5  •  1.0  (where  the  characteristic  length  was  chosen  to  be  the 
width  of  the  positive  correlation  region  in  the  onset 
boundary  layer) 

Similarly,  the  wing  geometry  independent  variables  (6/t,  a/b,  c/t, 
a  and  8)  were 

6/t  »  1.0  (onset  boundary  layer  thickness /wing  thickness  ratio) 
a/b  ■  1.5  (ratio  of  major/minor  axis  of  semi-elliptical  nose) 
c/t  ■  5  and  19.5  (chord  length/wing  thickness  ratio) 
a  ■  0  radians  (wing  angle  of  incidence) 

3  “  0.49  radians  (28  deg)  (trailing  edge  angle) 

6.1.3  The  Investigated  Flow  Regions 

The  large  scale  temporal  flow  structure  was  Investigated  in 
Regions  I,  II,  III,  IV,  and  VI,  which  were  defined  in  Section  2.2. 
The  investigations  were  performed  with  two  geometrical  chord 
lengths  for  the  wing,  described  in  Section  3.2  and  shown  in  Figure 
6-2.  In  the  first  case,  the  chord  length  to  wing  thickness  ratio 
c/t  was  19.5.  This  configuration  was  designated  the  infinite  chord 
length  wing,  and  was  used  to  eliminate  trailing  edge  effects  on  the 
flow  in  the  wing-body  junction.  It  was  assumed  that  the  relatively 
large  distance  between  the  nose  and  the  tall  would  allow  this 
assumption.  In  the  second  case,  the  value  of  c/t  was  5,  and  the 
effects  of  the  trailing  edge  were  studied  both  upstream  along  the 
side  of  the  wing  and  in  the  downstream  wing-body  junction  wake 
flow.  This  configuration  is  designated  the  finite  chord  length 
wing.  The  streamwise  locations  xQ  of  the  planes  (y/t,  z/t)  in 
which  the  two  point  velocity  correlations  were  measured  for  each 
configuration  nre  identified  in  Figure  6-2  by  the  arrow  symbols. 

6.2  Spatial  Structure  of  Onset  Boundary  Layer  (Region  I) 

Measurements  of  the  two-point  velocity  correlations  for 
lateral  (z/t)  and  transverse  (y/5)  separations  of  the  probes  were 
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Figure  6-2:  Schematic  Showing  Locations  of  Measurement  Planes  for  th 
Infinite  and  Finite  Chord  Length  Wings 


made  with  the  wing  removed  to  characterize  the  inherent  spatial 
flow  structure  of  the  onset  boundary  layer  (Region  I). 

The  correlations  were  made  with  the  fixed  probe  at  a 
tranverse  location  within  the  onset  boundary  layer  (y/6  -  constant) 
but  at  two  different  locations  on  the  wall.  Measurements  were  made 
at  the  two  locations  in  order  to  show  that  the  results  were  inde¬ 
pendent  of  wall  position  and  therefore  representataive  of  classical 
flate  plate  data.  The  results  for  PA  (-1.00,  0.22,  -0.25)  are 
shown  in  Figure  6-3,  and  the  results  for  PA  (-0.45,  0.22,  0.30)  are 
shown  in  Figure  6-4.  The  correlation  contours  in  the  two  figures 
are  in  substantial  agreement;  most  importantly,  the  widths  of  the 
positive  and  negative  correlation  regions  are  similar.  The  data 
clearly  show  a  lateral  spatial  structure  in  the  flow  that  is  con¬ 
sistent  with  the  notion  of  adjacent  regions  of  flow  acceleration 
and  deceleration.  For  the  reasons  given  in  Section  5.6  the  region 
of  negative  correlation  is  believed  to  have  been  the  result  of 
fluid  dynamic  phenomena,  and  not  of  instrumentation  limitations  or 
endwall  boundaries  of  the  wind  tunnel. 

An  example  of  the  frequency  distribution  of  the  velocity 
power  spectra  *A  (f)  and  *B  (f),  the  coherence  and  relative 

phase  -  $A  are  shown  in  Figure  6-5  for  velocity  correlations 
at  PA  (-1.0,  0.22,  -0.25)  and  PB  (-1.0,  0.25,  +  0.25).  These 
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positions  were  negatively  correlated  with  Y  ^  >  0.10  over  the 
coherent  frequency  range  2-65  Hz.  Note  that  the  coherent  structure 
is  not  immediately  evident  with  a  large  peak  in  the  power  spectra 
♦A  (f)  and  *B  (f). 

Additional  frequency  distributions  of  velocity  correlation 
information  are  presented  in  Figure  6-6  where  similar  information 
to  that  described  above  is  shown  for  PA  (-0.45,  0.22,  0.30)  and  PB 
(-0.45,  0.35,  0.80)  at  two  speeds  UQ  -  30.5  m/s  and  UQ  -  15.2  m/s. 
The  noteworthy  feature  is  the  shift  in  the  frequency  band  of  the 
coherent  fluctuations  to  lower  values  in  proportion  to  the  lower 
velocity;  this  result  indicates  a  fluid  dynamic  source  for  the 
measured  fluctuations. 

The  effect  of  the  transverse  position  (y/6)  of  the  fixed 


Figure  6-3: 
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probe  P A  on  the  correlation  contours  is  shown  in  Figures  6-7  and 
6-8  with  PA  (-0.45,  0.40,  0.30)  and  PA  (-0.45,  0.60,  0.30).  These 
figures,  in  conjunction  with  the  above  presented  Figure  6-4  with  PA 
(-0.45,  0.22,  0.30),  show  two  features.  First,  the  lateral  corre¬ 
lations  were  strongest  for  position  PA  at  y/6  <  0.5.  And,  second, 
non-zero  correlations  were  obtained  outside  the  boundary  layer 
for  positions  PA  at  y/6  >  0.5.  This  result  most  likely  reflects 
the  interaittancy  of  the  edge  of  the  boundary  layer. 

6.3  Spatial  Structure  of  Wing-Body  Junction  Flow 

6.3.1  Infinite  Chord  Length  Wing 

The  spatial  structure  of  the  wing-body  junction  flow  was 
investigated  using  the  infinite  chord  length  wing.  Measurements 
were  made  in  the  planes  Indicated  in  Figure  6-2  for  Regions  II,  III, 
and  IV. 

6. 3. 1.1  Region  II  (Pressure  Driven  Boundary  Layer)(-10  <  x/t  <  -0.8) 

Figure  6-9  shows  the  correlation  contours  in  Region  II  for 
the  fixed  probe  at  PA  (-1.00,  0.25,  -0.25)  and  the  traversed  probe 
at  PB  (-1.00,  y/6,  z/t).  The  measurement  plane  was  immediately 
upstream  from  the  location  of  the  mean  flow  leading  edge  separa¬ 
tion.  The  correlation  contours  in  Figure  6-9  are  similar  to  the 
correlation  contours  for  the  onset  boundary  layer,  Figures  6-3  and 
6-4.  The  noticeable  difference  is  that  the  regions  of  positive  and 
negative  phase  are  extended  in  the  lateral  dimension  because  the 
flow  slowed  and  spread  out  to  pass  around  the  wing.  Figure  6-10 
shows  example  spectral  distributions  of  the  power  spectral  den¬ 
sities,  the  coherence,  and  the  relative  phase  for  the  travesed 
probe  at  Pg  (-1.00,  0.25,  0.25).  The  distributions  are  similar  to 
those  obtained  for  the  onset  boundary  layer  shown  in  Figure  6-5. 

6.3.1. 2  Region  III  (Leading  Edge  Separation  Region)(-0.8  T  x/t  <  4) 

Figure  6-11  presents  the  correlation  contours  in  Region  III 
with  the  fixed  probe  at  PA  (-0.12,  0.25,  -0.50).  This  position  was 
immediately  downstream  from  the  leading  edge  mean  flow  separation 
but  upstream  from  the  ving  leading  edge.  The  positive  and  nega- 
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tive  correlation  regions  were  extended  even  further  in  the  lateral 
direction  than  they  were  in  Region  II.  The  extension  is  about  50 
percent  of  the  lateral  distance  for  the  onset  boundary  layer  in 
Region  I  (Figures  6-3  and  6-4).  By  comparing  Figures  6-9  and  6-11 
It  is  seen  that  there  was  no  abrupt  distortion  of  the  correlation 
contours  as  the  flow  passed  through  the  leading  edge  separation 
where  the  horseshoe  vortex  was  formed.  It  is  tentatively  concluded 
that  the  horseshoe  vortex  introduced  no  new  frequencies  or 
magnified  amplitudes  at  its  origin  in  the  flow  separation 
immediately  upstream  from  the  wing  leading  edge. 

Although  the  leading  edge  separation  did  not  immediately 
produce  gross  distortion  of  the  correlation  contours,  the 
coexisting  horseshoe  vortex  did  eventually  distort  the  contours 
downstream  but  still  inside  Region  III.  As  the  relatively 
undisturbed  flow  passed  around  the  wing  leading  edge,  the  flow  was 
inverted  by  the  lateral  and  transverse  velocities  of  the  vortex  to 
produce  the  contours  shown  in  Figures  6-12  and  6-13.  In  these 
figures,  the  correlations  are  for  both  the  fixed  and  the  traversing 
probes  on  the  same  side  of  the  wing.  Figure  6-12  shows  the  results 
obtained  for  the  fixed  probe  near  the  wing,  PA  (2.00,  0.15,  0.15) , 
and  Figure  6-13  for  the  probe  further  from  the  wing,  PA  (2.00, 

0.20,  0.80).  Notable  In  the  figures  is  the  absence  of  the  nearly 
vertical  division  between  the  regions  of  positive  and  negative 
correlation  that  was  observed  for  the  onset  boundary  layer. 

Evident  In  both  figures  is  the  effect  of  the  "jetting"  flow  pro¬ 
duced  by  the  streamwlse  vorticity  as  It  circulated  axial  momentum. 
From  a  comparison  of  the  correlation  contours  in  Figures  6-12  and 
6-13  with  the  velocity  contours  shown  in  Figure  4-15,  it  is  observed 
that  there  was  a  "jetting"  of  both  the  correlation  and  velocity  con¬ 
tours  from  the  high  momentum  (velocity)  region  into  the  contours 
for  the  lower  momentum  (velocity).  It  is  tentatively  concluded 
that  the  distortions  of  the  correlation  contours  were  produced  by 
vortical  transport  of  the  fluid  by  the  horseshoe  root  vortex. 

A  new  feature  of  the  flow  was  observed  for  the  correlations 
shown  in  Figure  6-12:  the  "finger"  of  the  positively  correlated 
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region  jetting  up  into  the  negative  correlation  region  exhibited  a 
frequency  bandwidth  different  from  that  in  the  rest  of  the  measurement 
plane.  Figure  6-14  shows  spectral  measurement  results  for  the  two 
frequency  bandwidths.  The  spectrum  in  Figure  14(a)  for  the  positi¬ 
vely  correlated  region  with  Pg  (2.00,  0.25,  0.15)  is  similar  to 
that  for  the  onset  boundary  layer  with  Bandwidth  I  (2-65  Hz).  The 
spectrum  shown  in  Figure  14(b)  from  the  region  of  the  jetting 
"finger"  at  Pg  (2.00,  0.25,  0.75)  in  the  negatively  correlated 
region  exhibits  a  new  frequency  bandwidth  (25-80  Hz)  which  is 
denoted  as  Bandwidth  II.  It  is  not  clear  from  the  measurements  if 
the  origin  of  this  higher  frequency  bandwidth  was  in  fact  a  result 
of  the  presence  of  the  wing  in  the  wall  boundary  layer.  It  is 
possible  that  it  was  a  flow  structure  inherent  in  the  onset  boun¬ 
dary  layer  inner  wall  region  that  was  displaced  away  from  the  wall 
by  the  horseshoe  vortex  transport.  In  either  case,  it  was  a  newly 
created  frequency  band  in  the  flow  structure  observed  at  positions 
away  from  the  wall,  and,  is  a  consequence  of  the  Interaction  of  the 
flow  in  the  wing-body  junction  region. 

The  meaning  of  the  negative  relative  phase  less  than  -180 
degrees  for  the  correlation  contours  at  positions  z/t  >  0.8  in 
Figure  6-13  may  be  that  there  was  a  difference  in  the  mean  veloci¬ 
ties  or  convection  paths  for  parts  of  flow  structure  originating 
upstream  and  traveling  to  positions  P^  and  Pg.  The  structure  parts 
could  then  have  passed  Pg  after  first  passing  P^  (the  fixed  probe 
position).  This  explanation  is  reasonable  if  the  flow  at  Pg  was 
either  at  a  slower  speed  or  a  longer  distance  in  its  convection 
from  the  origin  of  the  disturbance. 

In  addition  to  the  above  correlation  measurements  made  with 
the  two  probes  on  the  same  side  of  the  wing,  a  series  of  correla¬ 
tion  contours  were  measured  at  positions  Pg  (2.00,  y/t,  z/t)  in 
Region  III  with  the  fixed  probe  on  the  other  side  of  the  wing 
at  six  positions  inside  the  spatial  zone  (-0.75  z/t  _<  -0.15)  and 
(0.15  _<  y/6  _<  0.60).  The  results  for  these  "across-wing"  correla¬ 
tions  are  presented  in  Figures  16-15  through  6-22.  The  figures 
show  that  the  flow  was  correlated  for  positions  of  the  traversing 
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probe  as  far  from  the  wing  as  z/t  ■  0.75.  The  correlation  contours 
show  the  distortion  effect  of  the  vortlclty  transport  described 
above,  where  the  correlations  were  swept  out  from  the  Intersection 
of  the  wing  and  wall  In  a  direction  along  the  wall,  followed  by  a 
lifting  away  from  the  wall.  The  correlations  were  greatest  in 
magnitude  and  spatial  extent  for  the  fixed  probe  near  the  wing  and 
wall  intersection,  and  decreased  as  the  fixed  probe  was  moved  in 
the  transverse  or  lateral  directions. 

Examples  of  the  "across-wing"  spectral  distributions  of  the 
power  spectral  density,  the  coherence,  and  the  relative  phase  are 
given  in  Figures  6-23  and  6-24.  The  fixed  and  traversed  probes 
were  at  symmetric  positions  relative  to  the  wing  for  the  measure¬ 
ments  presented  in  Figure  6-23,  PA  (2.0,  0.15,  -0.15)  and  PB  (2.0, 
0.15,  +0.15).  The  correlations  are  in-phase  and  within  the 
Bandwidth  I  (2-65  Hz).  In  Figures  6-24  the  probe  positions  were 
located  at  non-symaetrlc  positions  with  PB  further  away  from  the 
wing  and  near  the  boundary  between  the  positive  and  negative  corre¬ 
lations.  These  locations  were  PA  (2.00,  0.15,  -0.15)  for  the  fixed 
probe  and  both  PB  (2.00,  0.10,  0.55)  (Figure  6-24(a))  and  PB  (2.00, 
0.20,  0.55)  (Figure  6-24(b))  for  the  traversed  probe.  In  Figure 
6-24(a),  the  correlations  were  in  Bandwidth  I  (2-65  Hz)  similar 
to  the  symmetric  condition.  Figure  6-24(b)  presents  evidence  that 
Bandwidth  11  (25-80  Hz)  was  present  with  small  coherence.  This  was 
the  only  measurement  in  which  Bandwidth  II  was  observed  for  the 
"across-wing"  correlations  in  Region  III. 

6. 3. 1.3  Region  IV  (Developed  Secondary  Flow)  (4  <  x/t  <  •) 

Additional  two-point  velocity  correlations  were  measured  in 
the  plane  at  the  streamwlse  location  at  xQ  -  6.05  in  Region  IV. 

The  "across-wing"  correlation  contours  for  the  traversing  probe  are 
presented  in  Figure  6-25  for  PA  (6.05,0.15,-0.15),  in  Figure  6-26 
for  PA  (6.05,0.40,-0.40),  and  in  Figure  6-27  for  PA  (6.05,  0.15, 
-0.40).  The  correlations  were  substantially  diminished  in  magni¬ 
tude,  but  the  effective  displacement  of  the  correlation  regions  by 
the  vortlclty  transport  was  greater,  than  for  the  contours  in  the 
upstream  Region  III  for  PA  (2.00,  y/5,  z/t).  This  observation  is 
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Figure  6-24:  Example  Across-Wlng  Spectral  Distributions  of  Correlation 

Measurement  Data  for  Non-Symmetrlc  Probe  Positions  In  Region  III 


made  by  comparing  Figures  6-25,  6-26,  and  6-27  with  Figures  6-15, 
6-17,  6-19,  respectively.  These  results  are  not  surprising;  the 
diminished  correlation  magnitudes  were  a  consequence  of  the  shear 
and  the  diffusive  forces  produced  by  the  vorticity  as  it  jetted 
high  axial  momentum  flow  into  the  low  momentum  region  of  the  flow. 

6. 3. 1.4  Summary  Comments 

Summary  comments  drawn  from  the  observations  of  the  flow  in 
the  wing-body  junction  are  first,  that  the  vortical  motion  caused 
by  the  horseshoe  vortex  effectively  displaced  the  correlation  con¬ 
tours  in  the  lateral  and  transverse  directions;  second,  that  the 
jetting  of  high  axial  momentum  fluid  into  the  low  momentum  boundary 
layer  resulted  in  partial  destruction  of  the  flow  structure;  and 
finally,  that  the  onset  boundary  layer  frequency  Bandwidth  I  was 
preserved  in  the  junction  flow  although  at  positions  away  from  the 
wall  a  new  higher  frequency  Bandwidth  II  evolved. 

6.3.2  Finite  Chord  Length  Wing 

Measurements  were  made  with  the  finite  chord  length  wing  to 
establish  the  effect  of  the  trailing  edge  condition  on  the  wing- 
body  junction  flow  and  to  observe  the  temporal  flow  structure  in 
the  wing-body  junction  wake  flow. 

6. 3. 2.1  Region  III  (Leading  Edge  Separation  Region)(-0.8  <  x/t  T  4) 

Two-point  velocity  correlations  were  repeated  at  the 
downstream  end  of  Region  III  at  the  plane  location  xQ/t  *  2.0  with 
the  finite  chord  length  wing.  The  purpose  was  to  examine  the 
effect  of  the  trailing  edge  condition  on  the  upstream  flow  struc¬ 
ture.  "Across-wing"  correlation  contours  are  presented  in  Figure 
6-28  for  the  fixed  probe  at  PA  (2.00,  0.15,  -0.15)  and  in  Figure 
6-29  at  PA  (2.00,  0.25,  -0.25).  The  fixed  probe  PA  was  located  on 
one  side  of  the  wing  while  the  traversing  probe  Pg  was  located  on 
the  other  side.  A  direct  comparison  can  be  made  between  the  con¬ 
tours  in  Figures  6-28  and  6-29  with  those  in  Figures  6-15  and  6-16 
for  the  infinite  chord  length  wing.  The  comparison  shows  that  the 
positive  and  negative  correlation  regions  for  the  finite  chord 
length  wing  were  smaller  in  spatial  extent  and  magnitude  than  those 
for  the  infinite  chord  length  wing.  This  result  was  expected 
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because  the  finite  chord  length  wing  experienced  flow  acceleration 
along  its  side  that  was  not  present  for  the  infinite  chord  length 
wing.  This  aspect  of  the  flow  will  be  discussed  in  greater  detail 
in  Chapter  7.  Another  observation  from  the  comparison  is  that,  in 
general,  the  effect  of  the  finite  chord  length  was  to  suppress 
positive  correlations  or  to  enhance  negative  correlations  perhaps 
because  of  a  change  in  the  mean  flow  for  the  horseshoe  vortex. 

Although  no  measurements  were  obtained  showing  a  distur¬ 
bance  in  Bandwidth  II  (25-80  Hz),  it  is  believed  that  this  result 
was  because  of  the  limited  number  and  locations  for  the  probe  posi¬ 
tions.  It  is  likely  that  a  more  complete  mapping  would  show 
Bandwidth  II,  which  did  reappear  in  measurements  downstream  to  be 
discussed  in  Section  6-4. 

No  other  new  features  were  apparent  in  the  spectral  distri¬ 
bution  of  the  correlation  data.  It  is  tentatively  concluded  that 
the  flow  in  Region  III  was  slightly  modified  by  the  change  in  the 
mean  flow,  and  that  the  temporal  structure  in  this  region  was  basi¬ 
cally  unchanged  by  the  trailing  edge  conditions  of  the  wing. 

6.4  Spatial  Structure  of  Wing-Body  Junction  Wake  Flow  (Region  VI) 

Measurements  were  continued  with  the  finite  chord  length 
wing  at  positions  in  a  plane  downstream  from  the  wing  trailing 
edge.  Those  measurements  were  made  in  Region  VI,  which  is  part  of 
the  wing-body  junction  wake  flow  and  clearly  the  most  complex 
region.  The  results  obtained  in  this  region  must  be  considered 
exploratory  since  more  systematic  measurements  are  required  to 
clarify  the  complex  nature  of  the  flow.  It  is  shown  that  the  wing- 
body  junction  wake  flow  is  distinguished  from  the  upstream  wing- 
body  junction  flow  both  in  the  spatial  distribution  and  frequency 
content  of  the  large  scale  temporal  flow.  The  active  Ingredient 
producing  the  new  features  in  the  flow  was  the  adverse  pressure 
gradient  in  the  trailing  edge  region  of  the  wing.  The  spatial 
structure  and  frequency  of  the  wake  flow  were  distinctly  different 
Inside  the  boundary  layer  than  outside  the  boundary  layer,  where 
for  the  latter  condition  ordinary  Strouhal-type  shedding  occurred 
with  a  much  higher  frequency  (See  Section  5.9).  The  spatial  struc¬ 
ture  of  this  complex  flow  is  shown  by  the  correlation  contours  in 
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Figure  6-30  for  the  fixed  probe  at  P^  (6.10,  0.25,  -0.38)  and  the 
traversing  probe  at  PB  (6.10,  y/6,  z/t). 

The  correlation  contours  show  that  the  flow  fluctuations 
were  antisymmetric  about  the  wing  plane-of-symmetry.  That  is,  the 
velocity  increased  on  one  side  of  the  wing  as  it  decreased  at  the 
reflected  position  on  the  other  side.  This  result  is  In  contrast 
to  the  correlation  contours  obtained  for  Regions  I-IV,  where  the 
flow  was  apparently  symmetric  although  positively  and  negatively 
correlated. 

Example  spectral  distributions  of  the  power  spectral  den¬ 
sities,  coherences,  and  relative  phases  are  shown  in  Figures  6-31 
(a),  (b),  and  (c)  for  three  positions  of  the  traversing  probe:  PB 

(6.10,  1.25,  -0.25),  PB  (6.10,  1.50,  0.25),  and  PB  (6.10,  0.25, 
0.25).  These  positions  as  well  as  the  fixed  probe  position  P^ 
(6.10,  0.25,  -0.38),  were  in  the  wing  "shadow"  in  the  downstream 
direction.  The  results  for  these  positions  PB  represent  measure¬ 
ments  in  three  correlation  regions.  As  shown  in  Figure  6-30,  there 
were  at  least  two  additional  correlation  regions  distinguished  by 
relative  phase,  one  of  which  was  outside  the  wing  "shadow",  another 
of  which  extended  from  inside  the  "shadow"  to  outside  the  "shadow." 
A  third  Bandwidth  III  (0.5-10  Hz)  was  observed  in  these  measure¬ 
ments  and  was  lower  than  either  of  the  two  bandwldths  observed  in 
Regions  I-IV.  Those  two  were  the  onset  boundary  layer  Bandwidth  I 
(2-65  Hz),  and  the  high  frequency  Bandwidth  II  (25-80  Hz)  that 
appeared  in  the  wing-body  junction  flow.  Apparently  the  trailing 
edge  condition  introduced  a  "natural  frequency"  into  the  spectral 
distribution  of  large-scale  temporal  flow  structure  that  gave  rise 
to  Bandwidth  III.  It  is  interesting  that  neither  the  onset  boun¬ 
dary  layer  frequency  Bandwidth  1  nor  the  wing-body  Bandwidth  II 
appeared  in  the  correlations  for  the  fixed  probe  at  P^  (6.10,  0.25, 
-0.38)  shown  in  Figure  6-30;  only  the  frequency  Bandwidth  III 
(0.5-10  Hz)  was  observed.  Except  in  the  vicinity  of  the  wake  cen¬ 
terline  where  values  of  the  power  spectral  density  were  large,  the 
power  spectral  density  for  Bandwidth  III  (0.5-10  Hz)  was  comparable 
in  magnitude  to  the  previous  measurements  for  Bandwidth  I  (2-65  Hz) 
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and  Bandwidth  II  (25-80  Hz).  Apparently  the  new  Bandwidth  III 
(0.5-10  Hz)  did  not  simply  drown  out  the  two  other  bands,  which  are 
apparently  missing  for  correlations  with  (6.10,  0.25,  -0.38). 

On  the  other  hand,  positioning  the  probe  PA  further  from 
the  wing  wake  in  the  lateral  direction  produced  correlations  for 
Bandwidth  II  (25-80  Hz)  as  well  as  Bandwidth  III  (0.5-10  Hz). 

Figures  6-32  and  6-33  show  the  correlation  contours  for  the  fixed 
probe  at  P^  (6.1,  0.25,  -0.75),  which  was  twice  as  far  from  the 
nominal  wake  centerplane  (z/t  *  0.0)  as  was  the  position  for  the 
contours  presented  in  Figure  6-30.  In  fact,  the  position  P^  (6.1, 
0.25,  -0.75)  was  outside  the  projected  profile  of  the  wing 
(thickness)  onto  the  measurement  plane. 

The  correlation  contours  shown  in  Figure  6-32  are  for  the 
frequency  Bandwidth  III  (0.5-10  Hz)  and  can  be  compared  to  the  con¬ 
tours  shown  in  Figure  6-30.  The  comparison  shows  that  the  correla¬ 
tions  are  quite  different.  With  the  fixed  probe  closer  to  the  wake 
centerline  the  correlation  regions,  shown  in  Figure  6-30,  were 
"pockets"  of  positive  and  negative  values  which,  for  one  thing, 
did  not  extend  completely  through  the  boundary  layer.  On  the  other 
hand,  the  correlations  obtained  with  the  fixed  probe  further  from 
the  wake  centerline,  shown  In  Figure  6-32,  clearly  show  that  both 
the  regions  of  positive  and  negative  correlation  extend  from  the 
wall  to  the  edge  of  the  boundary  layer. 

The  correlation  contours  shown  in  Figure  6-32  are  similar 
to  the  contours  shown  in  Figure  6-4  for  the  onset  boundary  layer. 
There  are  two  distinctions  however.  First,  the  frequency  Bandwidth 
III  (0.5-10  Hz)  for  the  wing-body  junction  wake  flow  (Figure  6-32) 
is  different  than  that  for  frequncy  Bandwidth  I  (2-65  Hz)  for  the 
onset  boundary  layer  (Figure  6-4).  Second,  there  were  three  correla¬ 
tion  regions  measured  in  the  direction  of  increasing  z/t  relative 
to  the  position  of  the  fixed  probe  in  the  wing-body  junction  wake 
flow,  while  the  onset  boundary  layer  measurements  showed  only  two 
regions . 

Figure  6-33  shows  the  correlation  contours  for  Bandwidth  II 
(25-80  Hz).  That  frequency  band  was  first  observed  in  the  wing- 
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Figure  6-32:  Contours  of  Constant  Sy  for  Region  VI  for  Bandwidth  III 
(P .  Outside  Wing  Shadow) 
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body  junction  flow.  Apparently  the  disturbance  associated  with 
that  frequency  band  persisted  into  the  wake  flow  and  grew  in  spa¬ 
tial  extent.  The  contours  in  Figure  6-33  for  Bandwidth  II  (25-80 
Hz)  are  not  similar  to  those  obtained  at  the  same  location  for  the 
Bandwidth  III  (0.5-10  Hz).  For  one  thing  the  width  of  the  positive 
region  surrounding  the  fixed  probe  is  twice  as  wide.  Furthermore, 
there  is  a  region  of  gradual  spatial  change  in  the  relative  phase 
between  the  large  region  of  positive  phase  and  the  small  region  of 
negative  phase.  There  are  only  two,  rather  than  three  correlation 
regions  in  the  direction  of  increasing  z/t. 

Figure  6-34  is  an  example  of  the  measurement  spectra  that 
clearly  shows  that  the  two  frequency  bands  described  above  were 
experimentally  distinguishable. 

The  measurements  indicate  that  there  were  two  kinds  of  flow 
structure  in  the  wing-body  junction  wake:  a  vertically  arranged 
structure  outside  the  wing  wake  proper  with  regions  of  positive  and 
negative  phase  each  extending  from  the  wall  to  the  edge  of  the 
boundary  layer,  and  an  antisymmetrically  arranged  structure  in  the 
wake  proper  that  was  characterized  by  pockets  of  positive  and  nega¬ 
tive  correlation.  The  vertically  arranged  structure  exhibited  fre¬ 
quencies  in  Bandwidth  II  (25-80  Hz)  and  Bandwlth  III  (0.5-10  Hz). 
The  spatial  patterns  of  the  correlation  regions  were  different  for 
the  two  bandwldths.  On  the  other  hand,  antisymmetrically  arranged 
structure  exhibited  frequencies  only  in  Bandwidth  III. 


Bandwidth  III  (0.5-10  ■*)  Bandwidth  II  (25-BO  B*) 


6-34 (a):  Pfi  (6.10,0.50,0.00) 


Figure  6-34:  Example  Spectral  Distributions  for  Correlation  Measurement 
Data  for  Region  VI  Showing  Bandwidth  Separation 
P.  (6.10,0.25,-0.75) 
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CHAPTER  7 


THE  MEASUREMENT  RESULTS:  INTERPRETATION, 
CONJECTURE,  CONCLUSIONS 


In  view  of  the  complexity  of  the  wing-body  junction  flow, 
and  the  limited  data  available  in  the  open  literature,  it  was 
expected  that  this  experimental  Investigation  would  raise  new 
questions  about  the  governing  physics  of  the  flow.  This  clearly 
has  been  the  case.  The  investigation  accomplished  the  objective  of 
showing  that  the  wing-body  junction  flow  is  characterized  by  large 
scale  time  dependent  flow  structure.  The  investigation  established 
an  orderly  framework  of  discussion  and  observation  around  which 
additional  research  can  be  performed  to  uncover  the  quantitative 
governing  physics.  As  shown  in  the  results,  distinct  features  of 
the  flow  structure  are  inherent  in  the  onset  boundary  layer,  and 
others  are  Introduced  by  the  wing-body  junction  flow,  and  still 
other  unique  features  are  created  in  the  wing-body  junction  wake 
flow  in  the  trailing  edge  region. 

The  results  show  that  the  frequency  and  spatial  extent  of 
the  structure  originated  at  the  trailing  edge  is  governed  in  part 
by  the  mean  flow  quantities  around  the  wing  nose  geometry.  The 


mechanics  of  the  structure  originated  in  the  trailing  edge  region, 
as  a  consequence  of  the  wing  geometry,  were  only  broadly  defined  by 
the  present  investigation.  Furthermore,  the  apparent  stability 
observed  in  this  investigation  of  the  flow  in  the  leading  edge 
region  should  not  be  construed  as  a  universal  feature  of  the  wing- 
body  junction  flow.  Rather,  the  conclusion  should  be  that  the  par¬ 
ticular  fixed  parameters  in  this  investigation  constrained  the 
motion  of  the  leading  edge  vortex  separated  flow. 

The  author  believes  that  variation  of  the  defining  parame¬ 
ters  will,  in  some  cases,  produce  an  unstable  leading  edge  vortex 
flow.  There  is  no  _a  priori  reason  that  flow  separation  in  that 
region  should  be  steady;  all  the  information  at  hand  suggests  that 
it  should  be  unsteady;  see,  e.g.  Simpson,  et  al.  (1981)  in  regard 
to  unsteady  two-dimensional  separation  in  an  adverse  pressure  gra¬ 
dient.  In  fact,  support  for  this  statement  is  found  in  the  present 
investigation  at  the  trailing  edge  of  the  wing.  An  observed 
separation  pattern  indicative  of  adverse  pressure  gradients  was 
associated  with  the  production  of  large  scale  time  dependent  flow 
structures  in  the  wing-body  junction  wake  flow. 

Inasmuch  as  the  measurements  require  additional  detail  to 
determine  the  quantitative  governing  physics,  the  discussion  of 
results  will  address  only  the  phenomenological  description  of  the 
flow.  The  notable  contribution  of  the  interpretation  presented  in 
this  chapter  is  to  provide  a  guide  for  future  research,  and  to  pro¬ 
vide  credibility  for  the  present  investigation  by  presenting  a  phy¬ 
sical  framework  for  discussion  in  which  the  measurements  are 
mutually  consistent. 

The  discussion  is  presented  in  three  sections.  The  first 
section  examines  the  meaning  and  limitations  of  the  two-point  velo¬ 
city  correlations  and  the  correlation  contours;  the  second  section 
discusses  the  physical  interpretation  of  the  flow  phenomena  that 
can  be  made  from  the  measurements;  and  finally,  the  third  section 
presents  a  list  of  conclusions  drawn  from  this  investigation. 

7,1  Interpretation  of  the  Two-Point  Correlation  Contours 

The  measurements  of  velocity  correlation  contours  presented 


in  this  investigation  require  careful  interpretation  because  of  the 
averaging  process  involved,  and  because  of  the  extrapolation  of 
two-point  statistics  to  multipoint  statistics  that  is  implied  by 
the  correlation  contours.  The  coherence  obtained  between  two 
measurement  positions  was  the  result  of  an  averaging  process  that 
necessarily  looses  Information.  For  example,  it  is  not  apparent 
whether  the  measured  frequency  bandwidth  was  produced  by  multi¬ 
frequency  structures  of  constant  amplitude  or  if  it  was  produced  by 
constant  frequency  structures  with  variable  amplitudes.  Neverthe¬ 
less,  the  Implication  of  non-zero  coherence  (the  quantity  from 
which  the  correlation  contours  were  derived)  is  that  there  was 
mutual  fluid  motion  at  the  two  measurement  positions.  It  is 
expected  that  the  averaging  process  understated  the  true  spatial 
correlation  of  Individual  flow  structures. 

The  correlation  coefficient  obtained  from  the  square  root 
of  the  coherence  is  a  measure  of  joint  probability  density  under 
Gaussian  conditions.  The  correlation  coefficient  is  very  simply 
interpreted  as  a  measure  of  the  fractional  part  of  the  velocity 
fluctuation  at  position  B  that  is  linearly  related  to  the  fluc¬ 
tuation  at  position  A.  But  this  interpretation  must  be  broadened 
by  a  reference  to  joint  probability  density  statistics  to  have  phy¬ 
sical  meaning. 

An  understanding  of  the  physical  meaning  of  the  concept  of 
a  correlation  coefficient  for  each  frequency  can  be  obtained  by 
examining  the  probability  density  associated  with  the  two-point 
velocity  correlations.  It  is  developed  in  this  section  that  the 
correlation  coefficient  specifies  both  the  slope  of  the  linear 
curve  fit  between  the  two  velocity  fluctuations  and  the  statistical 
distribution  of  the  velocities  about  the  linear  fit.  The  reference 
for  the  following  discussion  is  Jenkins  and  Watts  (1968,  p.67,  74, 
109).  Information  in  that  reference  has  been  specially  applied  by 
this  author  to  an  analysis  of  the  present  measurements. 

The  joint  probability  density  function  for  a  given  fre¬ 
quency  of  the  velocity  fluctuations  u^'and  uB '  at  positions  and 
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where  cr^,  ag  denote  the  root  mean  square  velocity  fluctuations  at 

uA'  uB' 

PA  and  Pb  for  given  frequency  f  and  pAB-  — —  is  the  corre- 
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correlation  coefficient.  In  the  case  where  the  fluctuations  are 

either  in-phase  or  out-of-phase,  p^g  -  Sy, 

The  argument  of  the  exponential  is  the  equation  for  an 

ellipse  oblique  to  the  u^'  and  Ug1  axis.  For  p^g  *  0,  the  equation 

describes  a  circle,  and  the  lines  of  constant  probability  density 

*AB  (“A**  UB*^  are  concentr^-c  circles.  For  non-zero  values  for 

PAB*  ttle  l*nes  constant  probability  density  form  ellipses 

rotated  in  the  u^'  and  Ug*  plane.  Figure  7-1  shows  the  contours 

2*°AaB^AB  ub* )  *  0.5  for  two  values  of  p^g,  0.0  and  0.5.  The 
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least  squares  linear  fit  of  uB  as  a  function  of  A  is  also 
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shown  in  the  figure;  the  slope  of  the  line  is  equal  to  the  value  of 
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Figure  7-2  shows  additional  significance  of  the  linear  fit 
determined  by  the  value  of  p^g.  The  several  shaded  bars  in  the 
figure  represent  ranges  of  values  for  Ug'  corresponding  to  a  cer¬ 
tain  probability  for  occurrence.  There  are  two  ranges,  or  bars, 
presented  for  each  of  five  values  for  u^'/o^,  (-1.0,  -0.5,  0.0, 
+0.5,  and  +1.0).  One  of  the  two  ranges  is  for  p^g  ■  0.0,  and  the 
other  is  for  p^  -  0.5;  the  two  ranges  are  noted  in  the  figure  and 
distinguished  by  the  shading  patterns.  The  ranges  represent  the 
interval,  centered  on  the  mean  value  of  Ug'/og',  in  which  70  per¬ 
cent  of  the  values  of  Ug'/o^  will  lie  when 
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Figure  7-2  shows  that  increasing  p^g  decreases  the  sta¬ 
tistical  spread  as  it  Increases  the  expected  mean  value  of  Ug'  for  a 
given  value  of  u^' .  Note,  however,  that  the  figure  cannot  be  used 
to  make  probability  predictions  at  the  different  values  for  u^';  to 
make  that  comparison  the  joint  probability  along  the  u^'  axis  would 
have  to  be  examined. 

A  physical  interpretation  of  the  meaning  of  Figure  7-2  is 
presented  in  Jenkins  and  Watts  (1974,  p.74),  where  typical  scatter 
diagrams  for  several  values  of  p^g  are  presented.  The  scatter 
diagrams  show  sets  of  discrete  measured  values  of  (u^’,  uB' )  that 
could  be  obtained  in  an  actual  situation.  Jenkins  and  Watts  state 
that  data  for  p^g  »  0.5  is  weakly  correlated. 

To  apply  this  statement  to  the  present  measurements  will 
require  additional  experimental  information.  The  statement  is 
surely  correct  considering  all  the  measurement  data  with  equal 
weighting.  However,  if  the  flow  structure  is  only  present  part  of 
the  measurement  time,  and  absent  for  the  remainder,  the  type  of 
measurement  presented  in  this  investigation  would  lead  to  weak 
correlations.  It  would  be  more  exact  and  revealing  to  measure 
correlation  only  while  the  structure  was  present  and  to  separately 
measure  the  fraction  of  time  that  it  is  present.  This  technique  of 
conditional  sampling  has  yet  to  be  developed  for  the  present 
problem. 

The  correlation  contours  seem  to  imply  that  the  flow 
through  the  measurement  plane  is,  on  the  average,  mutually  related 
at  arbitrary  positions  according  to  values  of  the  contours  at  the 
specific  positions.  This  is,  in  general,  not  correct  unless  cer¬ 
tain  other  conditions  are  present  because  the  contours  are  func¬ 
tions  of  the  positions  of  the  reference  probe.  For  example,  the 
correlations  AB^  and  the  AB2  do  not  imply  anything  about  the  corre¬ 
lation  B-^.  Hence,  the  contours  cannot  be  used  to  describe  mutual 
association  between  velocities  at  two  measurement  positions  neither 
of  which  is  at  the  reference  position.  In  fact,  study  of  the 
correlation  contours  shows  quite  explicitly  that  the  numerical  con¬ 
tours  can  be  quite  different  for  different  reference  positions. 
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Physical  meaning  in  terms  of  large  scale  movement  of  the 
flow  can  be  obtained  from  the  correlation  contours.  To  obtain  this 
meaning  one  must  have  certain  perceptions  of  the  flow.  These  per¬ 
ceptions  are  developed  below  in  two  steps.  First,  the  physical 
meaning  of  the  contours  in  the  onset  boundary  layer  is  presented. 

7.1.1  Contours  in  the  Onset  Boundary  Layer 

The  onset  boundary  layer  is  a  classical  equilibrium  boun¬ 
dary  layer  and  the  mean  flow  results  can  be  compared  with  those  in 
the  open  literature  (Section  4.1). 

The  correlation  contours  obtained  in  this  investigation  are 
also  consistent  with  results  obtained  previously  by  Favre,  et  al 
(1958).  The  measured  positive  correlation  contours  surrounding  the 
reference  probe  position  in  this  investigation  are  across-stream 
planar  cuts  from  the  streamwise  elongated  correlation  surfaces 
reported  by  Favre,  et  al.  These  elongated  surfaces  can  be  envi¬ 
sioned  as  representing  "cigar-shaped"  structures. 

Although  three-dimensional  correlation  contours  were  not 
obtained  in  the  present  measurements,  other  information  that  was 
obtained  can  be  used  to  demonstrate  that  the  flow  structure  was 
streamwise  elongated.  The  low  frequency  disturbance  bandwidth  for 
the  present  data  implies  a  long  correlation  distance  in  the  longi¬ 
tudinal  direction.  Example  auto-correlation  functions  at 
(-1.00,  0.25,  -0.25)  and  Pg  (-1.00,  0.10,  0.25)  are  shown  in  Figure 
7-3.  It  is  estimated  from  the  figure  that  flow  structure  existed 
at  each  measurement  position  for  0.75  seconds.  Invoking  Taylor's 
frozen  eddy  hypothesis  and  estimating  the  local  mean  flow  speed,  27 
m/s,  one  arrives  at  a  correlation  distance  equal  to  20m.  This 
length  is  two  orders  of  magnitude  larger  than  the  lateral  or  trans¬ 
verse  dimensions  of  the  structure;  hence  the  structure  is 
elongated. 

The  present  investigation  elaborated  on  the  phenomena 
reported  by  Favre  by  measuring  the  adjacent  negatively  correlated 
finite  flow  region  adjacent  to  the  elongated  positively  corre¬ 
lated  region.  The  appearance  of  negatively  correlated  regions  of 
the  flow  makes  sense;  a  thrust  of  fluid,  as  in  the  form  of  a  finger 
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of  jetting  fluid,  would  be  expected  to  be  accompanied  by  a  deceler¬ 
ation  of  the  fluid  on  either  side  to  satisfy  continuity  of  flow. 

This  result  is  not  new;  Champagne,  et  al  (1970)  measured  negative 
correlations  outside  the  positive  region,  i.e.,  at  distances  beyond 
the  first  zero  crossing  of  the  correlation  function,  but  the  corre¬ 
lation  function  measurement  technique  employed  at  that  time  was  not 
refined  enough  to  produce  accurate  results.  Furthermore,  the  nega¬ 
tive  correlations  were  observed  to  extend  to  the  wall  bounding  the 
fluid.  In  that  case  it  would  be  expected  that  the  large-scale  flow 
structure  was  representative  of  duct  flow  rather  than  boundary 
layer  flow  because  the  end  wall  conditions  contributed  to  the 
measured  flow. 

In  the  present  experiment  the  technique  employed  Fast 
Fourier  Transforms  and  solid  state  computer  micro-chip  electronics  to 
provide  both  nearly  real  time  low  frequency  analysis  and  a  significant 
improvement  in  accuracy.  Furthermore,  the  observed  negative  corre¬ 
lation  region  was  an  order  of  magnitude  smaller  than  the  wind  tun¬ 
nel  width.  The  small  size  of  the  negatively  correlated  region, 
relative  to  the  tunnel  width,  satisfied  one  requirement  of 
Comte-Bellot  and  Corrsin  (1971)  for  separating  fluid  dynamic  nega¬ 
tive  correlations  from  measurement  instrumentation,  or  network, 
negative  correlations  (see  Section  5.6  for  a  discussion). 

Furthermore,  the  finite  lateral  negative  correlation  region 
decoupled  the  measured  flow  from  the  tunnel  end  wall  boundary 
layers.  Therefore,  it  is  believed  that  the  measured  low  frequency 
"eddies"  were  not  products  of  duct  flow  but  rather  of  the  flat  plate 
boundary  layer. 

The  specific  correlation  contours  for  these  "eddies" 
generally  depend  on  the  location  of  the  fixed  probe.  It  is  dif¬ 
ficult  to  interpret  the  meaning  of  "superimposed"  two-point  cor¬ 
relations  from  two  reference  positions.  However,  certain 
conclusions  are  plausible  from  two-point  correlations  if  the  domain 
of  influence  of  the  velocity  at  one  reference  position  is 
coexistent  with  the  domain  at  another  position,  providing  the  two 
reference  positions  are  correlated.  Consider  the  correlation  con- 
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tours  for  the  onset  boundary  layer  in  Figures  6-4,  6-7,  and  6-8. 

The  figures  show  that  in  general  the  regions  of  positive  and  nega¬ 
tive  correlation  are  independent  of  the  position  y  for  reference 
positions  (xQ,  y,  zQ)  where  xQ  and  zQ  are  points  inside  the  boun¬ 
dary  layer. 

Under  these  conditions  it  is  plausible  that  the  flow  moved 
as  occasional  (quasi-random)  wholesale  movements  in  fairly  large 
regions  defined  by  the  signs  of  the  correlation  regions,  and  that 
this  movement  occurred  at  random  positions  in  the  boundary  layer. 
Furthermore,  the  movements  may  have  had  overlapping  structure. 

7.1.2  Contours  in  the  Wing-Body  Junction  Flow 

The  large  scale  time  dependent  flow  structure  in  the  onset 
boundary  layer  has  phenomenological  description  in  the  notion  of 
elongated  correlation  surfaces.  However,  that  description  falls  in 
the  wing-body  junction  flow  for  two  reasons.  First,  the  mean  vor¬ 
tex  flow  inverts  the  boundary  layer,  thereby  deforming  the 
elongated  correlation  surfaces  in  the  lateral  and  transverse  direc¬ 
tions.  Second,  the  Inversion  transports  high  energy  fluid  into  the 
lower  energy  region  of  the  boundary  layer.  There  is  a  longitudinal 
"jetting**  of  high  momentum  fluid  through  the  low  momentum  fluid. 
Hence,  the  correlation  contours  become  distorted  in  the  longitudi¬ 
nal  direction. 

The  result  is  that  a  material  line  of  particles  coinciding 
with  a  contour  of  constant  correlation  in  the  measurement  plane  is 
distorted  as  it  moves  downstream.  Original  material  planes  become 
distorted  surfaces.  Hence  it  would  be  expected  that  correlations 
in  (y,  z)  planes  would  decrease  as  x0  Increases.  This  apparent 
dissipation  is  in  addition  to  the  shear  and  diffusive  forces  which 
also  decrease  the  correlations. 

In  this  investigation  only  correlation  contour  measurements 
in  planes  were  made.  However,  no  apologies  are  made  for  not 
mapping  correlations  in  three-dimensional  space  because  the  precise 
manner  in  which  the  data  will  be  used  is  open  to  speculation. 

A  suggested  first  order  treatment  of  the  data  for  large 
scale  force  calculations  is  presented  in  the  next  section. 


7.1.3  Integrated  Flow  Structure 


A  numerical  measure  of  the  size  of  the  large  scale  flow 
structure  is  needed  to  visualize  its  importance  relative  to  the 
rest  of  the  flow  and  in  order  to  estimate  its  significance  in  engi¬ 
neering  calculations.  An  example  of  engineering  type  calculations 
is  the  vibratory  impact  of  the  flow  structure  on  rotating  machinery 
downstream  from  the  wing. 

The  reasoning  for  the  numerical  measure  proposed  below 
begins  with  the  notion  that  the  structure  is  elongated  in  the 
streamwise  direction.  The  present  investigation  shows  that  in 
almost  all  cases  the  correlations  were  either  almost  in-phase  or 
out-of-phase.  Because  the  large  scale  structure  had  a  wavelength 
much  larger  than  either  the  wing  chord  length  or  the  boundary  layer 
thickness,  it  is  estimated  that  the  longitudinal  distortion  of  the 
correlation  surface  produced  insignificant  relative  phase  changes 
in  planes  within  a  few  wing  chord  lengths  downstream  from  the 
trailing  edge.  This  estimate  may  not  be  accurate  further 
downstream. 

It  is  proposed  that  the  numerical  measure  of  the  size  of 
the  large  scale  flow  structure  be  the  area  bounded  by  the  Sy  ■  0.3 
correlation  contour,  where  the  reference  probe  is  placed  at  a  posi¬ 
tion  that  maximizes  the  area  within  the  0.3  correlation  contour 
surrounding  the  reference  probe.  The  wholesale  coherent  movement 
of  fluid  within  that  boundary  is  defined  to  be  the  integrated  flow 
structure  applicable  for  engineering  type  calculations.  The  choice 
of  the  0.3  correlation  contour  was  made  from  observations  that  the 
data  for  Sy  <  0.3  (coherence  ~  0.1)  were  noticeably  less  reliable 
than  data  for  Sy  >  0.3. 

An  example  application  of  the  Integrated  flow  structure  to 
the  data  in  the  measurement  plane  (6.05,  y/6,  z/t)  is  shown  in 
Figure  7-4.  In  this  figure,  which  was  prepared  from  Figure  6-30, 
it  is  assumed  that  the  reference  probe  was  at  the  maximal  position 
described  above. 

7.2  Conjecture  About  The  Large  Scale  Time  Dependent  Flow  Structure 

This  section  contains  speculative  arguments  which  are 
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intended  to  provide  a  framework  for  further  discussion  of  the  time 
dependent  flow  in  the  wing-body  junction.  Clearly  the  flow  is 
complex  and  a  complete  mapping  of  detailed  measurements  sufficient 
to  define  the  entire  flow  over  a  range  of  parameters  would  be  an 
expensive  undertaking.  Hence,  it  is  desirable  to  establish  a  model 
flow  in  terms  of  a  visual  construct  that  can  suggest  additional 
measurements  to  test  the  validity  of  the  model.  The  model  must  be 
continuously  modified  to  be  consistent  with  the  new  evolving  infor¬ 
mation.  At  any  time,  a  model  must  be  considered  consistent  if  it 
explains  the  data  in  hand  and  if  it  Is  otherwise  based  on  fundamen¬ 
tal  fluid  mechanics  theory  or  observation. 

The  phenomenological  description  presented  below  is  based 
on  the  measured  velocity  correlations,  mean  flow  velocity  and 
pressure,  and  flow  separation  patterns  obtained  during  the  present 
study.  However,  the  results  from  previous  studies  are  readily 
Introduced  where  appropriate  to  elaborate  on  the  flow  description. 

In  some  respects  the  description  to  follow  is  a  conjecture  by  the 
author  gleaned  from  all  the  bits  and  pieces  of  information,  both  in 
the  literature  and  from  actual  experiment,  that  were  obtained 
throughout  this  investigation. 

7.2.1  The  Onset  Boundary  Layer  (Region  I) 

The  observations  Indicate  that  the  onset  boundary  layer. 
Independent  of  the  presence  of  the  wing,  contains  large  scale  time 
dependent  flow  structure,  or  organized  motion  of  the  flow  that  is 
confined  both  in  time  and  in  space.  The  structure  is  convected  by 
the  mean  flow  and  produces  features  similar  to  those  expected  for  a 
longitudinal  acceleration,  or  jetting,  of  fluid  with  flow  decelera¬ 
tion,  or  reversal,  on  either  side  of  the  jet;  see  Figure  7-5.  The 
amplitude  of  the  fluctuating  velocity  in  the  structure  is  a  small 
fraction  of  the  mean  flow  speed.  Hence  the  structure  is  manifested 
as  a  small  perturbation  of  the  mean  flow.  The  structure  has  spa¬ 
tial  extent  with  height  and  width  on  the  order  of  the  boundary  lay¬ 
er  thickness,  and  length  an  order  of  magnitude  larger.  The  frequen¬ 
cy  for  the  structure  is  low  (f6/U0  *  0.04),  on  the  order  of  one- 
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fifth  the  typical  energy  carrying  frequencies  (f6/UQ-0.2);  see  Rao 
(1971). 

The  structure  is  similar  to  an  elongated  body  of  fluid, 
with  parts  further  from  the  wall  moving  faster  and  with  the 
velocity  at  the  center  of  the  body  fluctuating  in  the  opposite 
direction  to  that  at  the  outer  parts  in  the  lateral  directions. 

The  concept  of  a  jet  of  fluid  is  possibly  a  too  restrictive 
description  of  the  structure.  A  more  versatile  model  that  is  also 
consistent  with  the  measured  velocity  correlations  is  a  large  eddy 
motion  in  which  the  structure  consists  of  convected  counterrotating 
eddy  pairs,  or  even  higher  multiple  order  arrangements;  see  Figure 
7-6. 

The  readily  apparent  difference  between  the  two  models  is 
the  degree  of  cross-velocity  in  the  lateral  (z)  direction.  The 
larger  the  cross-velocity,  the  more  eddy-like  the  structure.  To 
assess  the  appropriateness  of  either  of  the  jet  and  eddy  models, 
detailed  measurement  of  the  spatial  and  temporal  features  of  the 
lateral  velocity  component  would  have  to  be  made  under  carefully 
controlled  test  conditions. 

For  both  models,  the  fluid  motion  is  vertically  structured. 

That  is,  it  has  nearly  constant  phase  in  the  transverse  (y)  direc¬ 
tion  in  the  boundary  layer,  i.e.,  from  the  wall  to  the  free  stream. 
The  phase  is,  however,  a  function  of  lateral  (z)  position  con¬ 
sistent  with  the  observed  flow  reversals. 

The  actual  motion  of  the  fluid  may  be  a  superposition  of  a 
large  number  of  the  flow  structures  described  above,  or  it  may  con¬ 
sist  of  a  small  number  of  structures  of  larger  amplitude.  The  sub¬ 
sequent  flow  produced  by  the  wing  may  depend  on  the  ’'density'*  of 
the  structures  in  the  onset  boundary  layer. 

7.2.2  The  Effect  of  the  Wing 

7. 2. 2.1  Region  II  (Pressure  Driven  Boundary  Layer) 

In  the  present  investigation  the  upstream  onset  boundary 
layer  structure  was  largely  unaffected  by  the  presence  of  the  pres¬ 
sure  gradient  produced  by  the  presence  of  the  wing.  In  fact  the 
structure  was  observed  to  maintain  its  spectral  distribution  even  as 
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It  passed  downstream  through  the  leading  edge  separation  region. 
This  observation  can  result  from  at  least  two  possible  descriptions 
of  the  flow  dynamics  depending  on  whether  the  flow  structures  are 
sparce  or  superimposed. 

In  the  first  description,  the  flow  structure  is  a  single 
organized  motion  that,  in  general,  passes  to  one  side  or  the 
other  of  the  pressure  ridge.  The  pressure  differential  across  the 
ridge  produces  an  oscillation  of  the  ridge  which  if  present  in 
this  study  was  not  sufficient  to  result  in  a  resonant  motion.  The 
second  description  envisions  a  superposition  of  many  structures 
randomly  positioned  in  space.  The  overlapping  structures,  passing 
to  either  side  of  the  pressure  ridge,  are  of  sufficient  number  that 
the  net  pressure  differential  across  the  pressure  ridge  is  almost 
zero. 

The  present  measurements  are  not  in  sufficient  detail  to 
resolve  which,  if  either,  of  the  descriptions  is  appropriate.  The 
reason  is  that  the  measurements  are  ensemble  averages  of  the  time 
dependent  characteristics  of  the  flow.  It  is  not  possible  to 
determine  if  the  measurements  were  of  successive  structures  passing 
the  probes,  or  of  superimposed  structures  continuously  passing  the 
probes.  Additional  Information  about  the  structure  will  have  to 
be  obtained  from  conditional  sampling  techniques.  It  is  conjecture 
by  this  investigator  that  the  true  answer  is  that  the  structure  is 
a  superposition  of  many  structures,  but  not  very  dense  or  else 
the  correlations  would  be  nearly  zero  for  even  very  small  distances 

7. 2. 2. 2  Region  III  (Leading  Edge  Separation  Region) 

These  slightly  overlapping  structures  pass  almost  intact 
through  the  separation  region  near  and  upstream  from  the  leading 
edge  of  the  wing.  However,  the  structure  is  stretched  in  the  lat¬ 
eral  direction  as  the  mean  flow  passes  to  either  side  of  the  wing. 
As  the  structure  initially  passes  to  the  side  of  the  wing  nose,  it 
is  "turned  over"  by  the  mean  cross-flow,  or  circulation,  in  the 
root  vortex.  The  mean  flow  circulation  transports  momentum,  first 
from  outside  the  boundary  layer  down  into  the  wing-wall  junction; 
secondly,  out  from  the  junction  along  the  wall;  and  finally,  up 
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into  Che  boundary  layer.  Because  of  this  transport  a  second  struc¬ 
ture  with  a  bandwidth  of  larger  frequencies  than  those  measured 
In  the  onset  boundary  layer  Is  literally  pumped  from  the  wall 
region  and  up  Into  the  middle  regions  of  the  boundary  layer.  It  Is 
not  clear  from  the  measurements  in  this  study  If  the  second  distur¬ 
bance  Is  produced  by  the  presence  of  the  wing  such  as  from  an 
oscillation  at  the  separation  line  produced  on  the  wall  by  the  root 
vortex.  This  author  suspects  that  It  is  the  structure  present  in 
the  wall  region  of  the  onset  boundary  layer  that  Is  transported 
away  from  the  wall  by  the  vortical  circulation.  Evidence  sup¬ 
porting  this  statement  is  the  small  peak  in  the  power  spectral  den¬ 
sity  in  this  higher  frequency  bandwidth  observed  for  probe 
positions  very  near  the  wall. 

The  velocity  correlations  are  positive  for  positions  both 
near  to  the  wing  and  positioned  symmetric  to  the  plane-of -symmetry 
of  the  wing.  In  effect  the  wing  simply  divides  the  onset  boundary 
layer  structure  in  a  way  that  does  not  destroy  its  basic  correla¬ 
tions,  except  very  near  the  wing  where  the  velocity  and  hence  the 
correlation  tend  to  zero  as  a  consequence  of  the  no-slip  condition. 
In  the  presence  of  the  wing,  two  fluid  elements  maintain  their  ori¬ 
ginal  correlation  though  the  wing  forces  them  further  apart. 

The  correlations  on  either  side  of  the  wing  are  arranged  in 
positive  and  negative  regions  that  are  generally  independent  of  the 
location  of  the  fixed  probe  within  a  fairly  large  latitude.  Hence 
the  flow  around  the  wing  has  large-scale  structure  that  is  in-phase 
along  the  sides  of  the  wing  and  out-of-phase  further  from  the  wing. 
The  significance  of  this  feature  is  that  the  wing  provides  an  abso¬ 
lute  reference  for  the  large-scale  structure,  whereas  there  was  no 
reference  in  the  onset  boundary  layer.  The  importance  of  this 
reference  is  that  there  will  be  definite  phase  relationships  bet¬ 
ween  the  fluctuations  on  one  side  of  the  wing  and  on  the  other 
side. 

In  the  case  of  an  infinite  wing  the  regions  of  positive 
correlation  remain  quite  large  on  either  side  of  the  wing;  they  are 
about  equal  to  the  boundary  layer  thickness  in  width.  However,  the 


finite  chord  length  wing  exhibits  touch  smaller  widths  for  the  posi¬ 
tive  regions,  which  introduces  the  possibility  that  the  flow  at  the 
two  sides  of  the  trailing  edge  can  be  out-of-phase. 

If  the  wing  chord  length  is  finite,  the  flow  is  at  higher 
velocity  along  the  sides  of  the  wing  than  it  is  for  an  infinite 
chord  length  wing.  The  velocities  at  the  sides  of  finite  chord 
length  wings  can  be  sufficiently  large  to  maintain  a  relatively 
large  vortlcity  intensity.  This  effect  is  that  the  vortical 
transport  can  be  sufficient  to  produce  out-of-phase  correlations  at 
symmetric  positions  about  the  wing  plane  of  symmetry;  whereas, 
in-phase  correlations  are  obtained  for  the  infinite  chord  length 
wing.  The  explanation  for  the  negative  correlations  is  not  that 
flow  originally  in-phase  suddenly  becomes  out-of-phase  flow  dose 
to  the  wing.  The  Intense  vortlcity  associated  with  the  finite 
chord  length  wing  transports  flow  toward  the  wing  that  was,  in  the 
onset  boundary  layer,  separated  by  a  larger  lateral  distance. 

That  larger  distance  means  that  the  flow  was,  in  the  onset  boundary 
layer  out-of-phase  (or  negatively  correlated).  Therefore  the  flow 
along  the  sides  of  the  wing  exhibits  out-of-phase  fluctuations  and 
appears  to  have  antisymmetric  fluctuations  about  the  wing  plane  of 
symmetry. 

7. 2. 2. 3  Region  IV  (Developed  Secondary  Flow) 

The  root  vortex  produces  large  variations  in  the  axial 
momentum  distribution  that  decrease  the  structural  features  in  two 
ways.  First,  the  vortex  motions  literally  produce  a  high  speed 
fluid  jet  in  the  boundary  layer  near  the  wing-body  intersection, 
and  a  low  speed  jet  in  the  outer  part  of  the  boundary  layer  away 
from  the  wing-body  intersection.  The  resulting  velocity  gradients 
produce  shearing  forces  that  dissipate  the  structure.  Second,  the 
inversion  of  the  boundary  layer  produces  distortion  of  the  struc¬ 
ture  in  the  lateral,  transverse,  and  longitudinal  directions.  The 
longitudinal  distortion  would  be  expected  to  decrease  two-point 
velocity  correlations  in  the  transverse  and  lateral  directions. 


7.2.3  The  Special  Effects  of  the  Trailing  Edge 

7. 2. 3.1  Region  V  (Trailing  Edge  Separation  Region) 

The  flow  apparently  undergoes  dramatic  changes  in  its 
characteristics  as  it  passes  through  the  trailing  edge  region  of 
the  wing.  In  this  region  the  mean  separation  line  on  the  wing, 
parallel  to  the  trailing  edge  outside  the  boundary  layer,  is  curved 
toward  the  trailing  edge  for  a  relatively  thick  wing  as  a  result  of 
the  root  vortex  transport  of  high  momentum  into  the  momentum  defi¬ 
cient  boundary  layer  to  delay  separation. 

It  is  believed  that  the  trailing  edge  separation  in  Region 
V  is  the  cause  for  the  radical  changes  in  both  the  frequency  band- 
widths  and  the  correlation  patterns  observed  in  Region  VI 
downstream  from  the  trailing  edge.  These  changes  may  be  the  result 
of  a  selection  process  in  which  there  is  enhancement  of  the  corre¬ 
lations  of  fluctuations  at  certain  frequencies  of  the  onset  distur¬ 
bance,  and  suppression  of  others.  The  enhancement  is  governed  by 
the  "natural  frequencies"  of  the  separations  on  the  wing  and  on  the 
wall. 

A  schematic  of  the  trailing  edge  separated  flow  is  shown  in 
Figure  7-7.  The  schematic  shows  the  following  major  features:  (1) 
the  ordinary  Strouhal-type  shedding  phenomenon  outside  the  boundary 
layer,  (2)  a  lower  frequency  shedding  Inside  the  boundary  layer 
associated  with  the  wing  separation  line  curvature  inside  the  boun¬ 
dary  layer  for  Region  V,  (3)  a  divergent  pair  of  separation  lines  on 
the  wall  for  Region  VI,  separating  the  wake  proper  from  the  high 
energy  outer  parts  of  the  wake,  and  (4)  the  rolled  over  vortex  sheet 
in  Region  VI  associated  with  the  wall  separation. 

7. 2. 3. 2  Region  VI  (Wake-Boundary  Layer  Region) 

The  divergent  separation  lines  on  the  wall  in  the  wing  wake 
are  dependent  on  the  vortical  transport  and  the  wing  separation 
lines  for  both  their  Included  angle  and  their  virtual  origin  (see 
Figure  7-7).  In  the  present  investigation  the  separation  line 
on  the  wing  curved  to  the  trailing  edge,  delaying  separation  on  the 
wall  so  that  the  separation  lines  on  the  wall  were  separated  by  an 
Included  angle  of  70  deg  and  had  a  virtual  origin,  from  which  the 
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Figure  7-7:  Schematic  of  Flow  Structure  and  Separations  at 
Wing  Trailing  Edge 
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lines  appeared  to  radiate,  a  distance  equal  to  one-half  the  wing 
thickness  upstream  from  the  wing  geometric  trailing  edge. 

Figure  7-8  shows  detail  of  the  conjectured  flow  in  the  vici 
nity  of  the  rolled  up  vortex.  The  detail  is  of  the  instantaneous 
flow  in  a  plane  normal  to  the  mean  flow;  l.e.  in  the  lateral  and 
transverse  directions.  The  significant  feature  of  the  detail  is 
the  saddle  point  produced  by  the  counterrotating  vortices  from 
the  horseshoe  vortex  and  the  vortex  sheet  on  the  wall.  Topological 
considerations  of  continuous  fluids  require  a  third  vortex  near 
the  wall,  as  depicted  in  the  figure.  The  third  vortex  contains 
fluid  with  low  longitudinal  momentum. 

The  above  features  are  conjectured  to  produce  the  antisym¬ 
metric  pockets  of  correlated  fluid  shown  in  Figures  7-9.  The 
figure  corresponds  to  measured  data  given  previously  in  Figure 
6-30.  Figure  7-9  contains  a  description  of  the  flow  accelerations 
and  decelerations,  and  the  corresponding  displacement  of  the  third 
(low  momentum)  vortex  producing  the  various  "pockets"  of  flow 
structure. 

7.3  Conclusions 

Consideration  of  the  evidence  obtained  in  this  investiga¬ 
tion  of  the  large  scale  temporal  flow  in  the  wing-body  junction 
supports  the  following  conclusions: 

(1)  The  onset  boundary  layer  temporal  flow  structure  passes  Intact 
although  greatly  deformed,  through  the  wing-body  junction  flow  to 
the  trailing  edge  region  of  the  wing,  where  it  encounters  adverse 
pressure  gradients  and  is  modified  or  replaced  to  form  new  flow 
structures  downstream  from  the  wing. 

(2)  For  the  particular  configuration  studied  the  wing-body  junc¬ 
tion  flow  is  characterised  by  three  disturbance  structures  distin¬ 
guished  by  frequency  bandwidth.  Bandwidth  I  is  initiated  in  the 
onset  boundary  layer  and  not  in  the  wing-body  junction  wake  flow. 
Bandwidth  II  may  also  be  present  in  the  onset  boundary  layer  at 
positions  close  to  the  wall  but,  in  any  event,  appears  away 


from  the  wall  In  the  boundary  layer  Interaction  flow  around  and 
downstream  from  the  wing.  Bandwidth  III  appears  downstream  in  the 
wing-body  junction  wake  flow  and  is  clearly  the  result  of  the  pre¬ 
sence  of  the  wing.  The  three  bandwldths  contain  frequencies  at 
least  an  order  of  magnitude  smaller  than  the  ordinary  Strouhal-type 
shedding  from  the  wing  trailing  edge  outside  the  boundary  layer. 

(3)  The  large-scale  structure  downstream  from  the  wing  (Bandwidth 
III)  is  dependent  for  its  frequency  on  the  nose  shape  of  the  wing, 
probably  because  of  its  effect  on  root  vortex  Intensity. 

(4)  The  horseshoe  vortex  apex  (origin  at  the  wall  separation 
upstream  from  the  wing  leading  edge)  does  not  exhibit  a  charac¬ 
teristic  or  resonant  frequency  at  least  for  the  set  of  parameters 
studied  in  this  Investigation. 

These  conclusions  on  the  temporal  flow  structure  should  be 
further  tested  to  ascertain  their  range  of  validity. 

In  the  course  of  this  investigation  limited  studies  were 
made  to  verify  the  general  features  of  the  mean  flow.  From  those 
studies,  the  following  additional  conclusions  were  made: 

(5)  The  mean  flow  in  the  wing-body  junction  is  characterized  by  a 
root  vortex  originating  with  the  wall  separation  Immediately 
upstream  from  the  wing  leading  edge,  along  with  various  other  forms 
of  three-dimensional  separations  originating  in  the  trailing  edge 
region. 

(6)  The  pressure  effects  of  the  wing  leading  edge  extend  upstream  a 
distance  on  the  order  of  10  wing  thicknesses. 
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